Ordering and Defect Dynamics in Passive and Active Nematopolars

Fabio Aprile; Giuseppe Gonnella; Massimiliano Semeraro

Ordering and Defect Dynamics in Passive and Active Nematopolars

T0 review · 3 major / 6 minor · reviewed 2026-07-09 · glm-5.2

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read One field, two symmetries: a minimal model for nematopolar ordering

desk verdict Single-field nematopolar model reproduces depolarization strings and coarsening; active case shows MICS and arrested coarsening but lacks finite-size scaling check read the letter →

arxiv 2607.07656 v2 pith:PS5PMOAX submitted 2026-07-08 cond-mat.soft

Fabio Aprile , Massimiliano Semeraro , Giuseppe Gonnella This is my paper · ORCID

classification cond-mat.soft

keywords defectsystemsdefectsdynamicsnematicorderingpolaractive

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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

What would settle it

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.

Watch

Extended reading notes

Core claim

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

Load-bearing premise

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

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

read the 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.

Referee Report

3 major / 6 minor

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.

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 (3)

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 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 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.

minor comments (6)

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.
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 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 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.
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.
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.

Simulated Author's Rebuttal

3 responses · 1 unresolved

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.

read point-by-point responses

Referee: 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.

Authors: 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: yes
Referee: 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.

Authors: 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: partial
Referee: 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.

Authors: 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: yes

standing simulated objections not resolved

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.

Circularity Check

1 steps flagged · score 2.0 of 10

Minor: unspecified monotonic function g(k̄) assumed increasing, then trend 'predicted' from that assumption; central claims rest on simulations, not this argument.

fitted input called prediction [Section IV.A, Equations 8 and 10, and surrounding text]
"E_string ~ σ ℓ g(k̄) ... g(k̄) is defined for k̄≥0, vanishes at k̄=0 (for which no string emerges) and increases monotonically. ... ℓ_eq ~ (k_p + 2k_n)/(σ g(k̄)) = ℓ_nat / h(k̄), where we used σ ~ α_p ℓ_nat f(k̄) and defined h(k̄) ≡ f(k̄)g(k̄), which is an increasing function of k̄. ... Equation 10 naturally implies that, when expressed in natural units ... the dimensionless separation ℓ_eq_nat ≡ ℓ_eq/ℓ_nat = 1/h(k̄) decreases with k̄, as shown in the inset of Figure 2(d)."

The function g(k̄) is not computed from the free energy; it is defined solely by the property that it 'increases monotonically' with k̄. Since f(k̄) = √(k̄/(k̄+2)) is also increasing, h(k̄) = f(k̄)g(k̄) is increasing by construction. The 'prediction' that ℓ_eq_nat = 1/h(k̄) decreases with k̄ is then a direct algebraic consequence of the assumed monotonicity of g. The observed trend in Figure 2(d) is not independently predicted; it is the assumption restated. The same h(k̄) reappears in Eq. 14 for the loop critical radius R_c, where the same logic applies: the predicted trend follows from the assumed property of g. However, the paper is transparent that g is an effective screening function, not a first-principles result, and frames these as qualitative rationalizations rather than定量预测. The纸

full rationale

The paper's central claims — (t/ln t)^{1/2} coarsening, MICS, arrested coarsening — are supported by numerical simulations (Figs. 4–6), not by the energetic arguments involving g(k̄) and h(k̄). The unspecified function g(k̄), defined only by its monotonicity, is used to 'predict' the observed decreasing trend of d_eq with k̄ (Eq. 10) and the phase boundary in the loop diagram (Eq. 14). These trends follow tautologically from the assumption that g is increasing. This is a minor circular element in a qualitative rationalization, not in a load-bearing derivation. The model (Eqs. 1, 3) is self-contained and newly introduced; no self-citation chain underpins the central results. The coarsening law is derived from standard KT arguments applied to the system and verified numerically. The MICS and arrest claims rest on simulation data, not on the g(k̄) argument. Score 2 reflects the minor self-referential character of the energetic 'predictions' without elevating to a structural problem.

Assumptions & free parameters 7 free parameters · 4 assumptions · 2 invented entities

The model introduces standard physical parameters (α_p, k_p, k_n, Λ, Γ) fixed to specific values. The heuristic functions g(k̄) and h(k̄) are free parameters in the analytical arguments, defined only by qualitative properties to fit observed trends. The axioms are standard domain assumptions for dry active matter. No fundamentally new physical entities are postulated beyond named phenomenological structures (strings, MICS) that are observed in simulation.

free parameters (7)

α_p = 0.1
Bulk polarization parameter controlling the double-well potential depth.
k_p = varied (e.g., 0.03)
Polar elastic constant; varied via k̄ while keeping k_p fixed in some sweeps.
k_n = varied (e.g., 0.1)
Nematic elastic constant; varied to tune k̄.
Λ = 0.1 (active case)
Self-advection strength controlling activity.
Γ = 1
Rotational viscosity.
c_p = c_n = 0.2
Threshold for extracting correlation lengths L_p(t) and L_n(t).
g(k̄), h(k̄) = unspecified
Heuristic functions introduced in Eqs. 8-10 to parameterize string energy and equilibrium separation; defined only by monotonicity, used to fit qualitative trends.

assumptions (4)

domain assumption Model A (non-conserved) relaxational dynamics for the polar field p (Eq. 3).
Section II: assumes purely relaxational dynamics for the order parameter.
domain assumption Single elastic constant approximation for polar and nematic gradients.
Section II, Eq. 1: assumes k_p and k_n are constants independent of direction.
domain assumption Dry active matter limit: activity enters only via self-advection Λ(p·∇)p, neglecting hydrodynamic velocity fields.
Section II: the model omits solvent hydrodynamics, assuming the dry limit is sufficient.
domain assumption Two-dimensional system (d=2).
Section II: all simulations and derivations restricted to 2D.

invented entities (2)

Motility-induced charge symmetry breaking (MICS) independent evidence
purpose: Describes the activity-driven selection of positive integer and negative half-integer defects.
Observed in simulations (Fig. 6) as a persistent steady-state phenomenon; falsifiable in active matter experiments.
Depolarization strings independent evidence
purpose: Lines of vanishing polarization connecting half-integer defects.
Directly visualized in simulations (Fig. 1, 2) and consistent with prior experimental/theoretical reports in ferroelectric liquids.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Ordering and Defect Dynamics in Passive and Active Nematopolars." pith.science (2026). https://pith.science/paper/PS5PMOAX

@misc{pith2026260707656,
  author       = {Pith},
  title        = {Pith review of: Ordering and Defect Dynamics in Passive and Active Nematopolars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PS5PMOAX}},
  note         = {Machine review of arXiv:2607.07656}
}

abstract

The coexistence of polar and nematic interactions, observed in a broad range of biological and synthetic active systems, gives rise to a rich phenomenology that continues to challenge our theoretical understanding of non-equilibrium collective behaviour. In this paper, we numerically investigate phase ordering and defect dynamics in a newly introduced minimal single-field model for dry nematopolar systems, where competing polar and nematic contributions enter the free energy, and activity is implemented through self-advection. At optimal balance, the system develops depolarization strings connecting half-integer defects and separating domains with opposite polarization, together with closed depolarization loops. We first characterize the elementary relaxation mechanisms of defect pairs and loops, showing that the interplay between polar and nematic alignment gives rise to non-monotonic string-mediated interactions, finite equilibrium separations and distinct loop-collapse pathways. Large-scale simulations from disordered states instead show dynamic scaling with a characteristic length growing as $\sim(t/\ln t)^{1/2}$, consistent with coarsening in systems with non-conserved order parameters and point-like defects. Upon introducing self-advection, sufficiently strong activity leads to the coexistence of positive integer and negative half-integer defects, which we term motility-induced charge symmetry breaking, and to saturation of the characteristic length scales, ultimately resulting in arrested coarsening. Overall, our results provide a simple unified framework for understanding the ordering and defect dynamics in biological and synthetic nematopolar systems.