Pith. sign in

REVIEW 3 major objections 6 minor 63 references

One field, two symmetries: a minimal model for nematopolar ordering

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · glm-5.2

2026-07-09 03:32 UTC pith:PS5PMOAX

load-bearing objection Single-field nematopolar model reproduces depolarization strings and coarsening; active case shows MICS and arrested coarsening but lacks finite-size scaling check the 3 major comments →

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

Ordering and Defect Dynamics in Passive and Active Nematopolars

Fabio Aprile , Massimiliano Semeraro , Giuseppe Gonnella This is my paper
classification cond-mat.soft
keywords defectsystemsdefectsdynamicsnematicorderingpolaractive
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

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.

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

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

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

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 this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

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

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)
  1. 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.
  2. 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.
  3. 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)
  1. 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.
  2. 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.
  3. 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.
  4. 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.
  5. 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.
  6. 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

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

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

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

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

specific steps
  1. 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.

Axiom & Free-Parameter Ledger

7 free parameters · 4 axioms · 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.
axioms (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.

pith-pipeline@v1.1.0-glm · 23342 in / 2688 out tokens · 558841 ms · 2026-07-09T03:32:06.662379+00:00 · methodology

0 comments
read the original 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 a self-advection contribution. At optimal balance between the two alignments, 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.

Figures

Figures reproduced from arXiv: 2607.07656 by Fabio Aprile, Giuseppe Gonnella, Massimiliano Semeraro.

Figure 1
Figure 1. Figure 1: A typical nematopolar configuration with string and loop structures. (a) Magnified view of a representative nematopolar configuration at t = 103 for the reference case k¯ = 0.3, started from a disordered initial state. (b)–(d) Defect structures emerging during evolution: a string connecting two like-charged half-integer defects (b), a string connecting oppositely charged half-integer defects (c) and a clos… view at source ↗
Figure 2
Figure 2. Figure 2: Pairs of half-integer defects in controlled settings. (a) Magnified view of a pair of oppositely charged defects evolved in the reference case k¯ = 0.3 and connected by a string. The corresponding far-field configuration displays an ordered arrangement, with arrows all pointing upwards. Inset: same defect pair evolved in the case k¯ = 0.05. (b) Annihilation velocity v of oppositely charged defect pairs as … view at source ↗
Figure 3
Figure 3. Figure 3: Loop phenomenology: shrinking, rotational relaxation and rupture. (a) to (c) Evolution of a depolarization loop at (k¯ = 0.2, ℓnat = 3.0). Panels report successive snapshots of the gradually shrinking loop. (d) Relaxation of a depolarization loop through a continuous rotation of the polarization field at (k¯ = 0.6, ℓnat = 3.5). (e) to (f) Rotational relaxation dynamics, followed by loop rupturing into two … view at source ↗
Figure 4
Figure 4. Figure 4: Correlation functions and typical lengths. (a) Nematic correlation function Cn(r, t) for the pure nematic case. (b) and (c) Polar and nematic correlation functions Cp(r, t) and Cn(r, t) for the reference case. Insets report the curve collapse obtained by rescaling distances according to the relevant length scale Ln(t) or Lp(t) from (d). (d) Typical nematic Ln(t) (main) and polar Lp(t) (inset) lengths for d… view at source ↗
Figure 5
Figure 5. Figure 5: Defect elimination dynamics in the ordering regime. (a)Time trend of the ratio between the difference and the sum of the number of integer and half-integer defects of both signs in nematopolar configurations for different k¯. (b) Total number of defects n tot def, given by the sum of integer and half-integer ones, as a function of time for increasing values of k¯. As reference, the black line reports the c… view at source ↗
Figure 6
Figure 6. Figure 6: Nematopolar ordering with advection. (a) Magnified view of a representative non-equilibrium nematopolar config￾uration at t = 5 × 103 with positive advection strength, started from a disordered initial state and depicted with the same representation style as in [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages · 1 internal anchor

  1. [1]

    Asymptotic structure factor and power-law tails for phase ordering in systems with continuous symmetry,

    A. J. Bray and S. Puri, “Asymptotic structure factor and power-law tails for phase ordering in systems with continuous symmetry,” Phys. Rev. Lett.67, 2670–2673 (1991)

    work page 1991

  2. [2]

    Theory of phase-ordering kinetics,

    A. J. Bray, “Theory of phase-ordering kinetics,” Adv. Phys.43, 357–459 (1994). 16

    work page 1994

  3. [3]

    Theory of dynamic critical phenomena,

    P. C. Hohenberg and B. I. Halperin, “Theory of dynamic critical phenomena,” Rev. Mod. Phys.49, 435–479 (1977)

    work page 1977

  4. [4]

    Coarsening phenomena,

    L. F. Cugliandolo, “Coarsening phenomena,” C. R. Phys.16, 257–266 (2015)

    work page 2015

  5. [5]

    Critical phenomena and renormalization-group theory,

    A. Pelissetto and E. Vicari, “Critical phenomena and renormalization-group theory,” Phys. Rep.368, 549–727 (2002)

    work page 2002

  6. [6]

    The critical properties of the two-dimensional xy model,

    J. M. Kosterlitz, “The critical properties of the two-dimensional xy model,” J. Phys. C7, 1046 (1974)

    work page 1974

  7. [7]

    Breakdown of scaling in the nonequilibrium critical dynamics of the two- dimensional XY model,

    A. J. Bray, A. J. Briant, and D. K. Jervis, “Breakdown of scaling in the nonequilibrium critical dynamics of the two- dimensional XY model,” Phys. Rev. Lett.84, 1503–1506 (2000)

    work page 2000

  8. [8]

    Persistence exponents and scaling in two-dimensionalXYmodel and a nematic model,

    Subhrajit Dutta and Soumen Kumar Roy, “Persistence exponents and scaling in two-dimensionalXYmodel and a nematic model,” J. Phys. A Math. Gen.38, 5859–5868 (2005)

    work page 2005

  9. [9]

    Kinetics of phase ordering in uniaxial and biaxial nematic films,

    M. Zapotocky, P. M. Goldbart, and N. Goldenfeld, “Kinetics of phase ordering in uniaxial and biaxial nematic films,” Phys. Rev. E51, 1216–1235 (1995)

    work page 1995

  10. [10]

    Phase ordering in nematic liquid crystals,

    C. Denniston, E. Orlandini, and J. M. Yeomans, “Phase ordering in nematic liquid crystals,” Phys. Rev. E64, 021701 (2001)

    work page 2001

  11. [11]

    Theory of hybrid defects, with coupled orientational order parameters, on flat and curved surfaces,

    L. Paik and J. V. Selinger, “Theory of hybrid defects, with coupled orientational order parameters, on flat and curved surfaces,” Soft Matter22, 4151–4160 (2026)

    work page 2026

  12. [12]

    First-principles experimental demonstration of ferroelectricity in a thermotropic nematic liquid crystal: Polar domains and striking electro-optics,

    X. Chen, E. Korblova, D. Dong, X. Wei, R. Shao, L. Radzihovsky, M. A. Glaser, J. E Maclennan, D. Bedrov, D. M. Walba, and N. A. Clark, “First-principles experimental demonstration of ferroelectricity in a thermotropic nematic liquid crystal: Polar domains and striking electro-optics,” Proc. Natl. Acad. Sci.117, 14021–14031 (2020)

    work page 2020

  13. [13]

    Ferroelectric nematic liquid crystal, a century in waiting,

    O. D. Lavrentovich, “Ferroelectric nematic liquid crystal, a century in waiting,” Proc. Natl. Acad. Sci117, 14629–14631 (2020)

    work page 2020

  14. [14]

    Soliton walls paired by polar surface interactions in a ferroelectric nematic liquid crystal,

    B. Basnet, M. Rajabi, H. Wang, P. Kumari, K. Thapa, S. Paul, M. O. Lavrentovich, and O. D. Lavrentovich, “Soliton walls paired by polar surface interactions in a ferroelectric nematic liquid crystal,” Nat. Commun.13, 3932 (2022)

    work page 2022

  15. [15]

    Ferroelectric nematic liquids with conics,

    P. Kumari, B. Basnet, H. Wang, and O. D. Lavrentovich, “Ferroelectric nematic liquids with conics,” Nat. Commun.14, 748 (2023)

    work page 2023

  16. [16]

    Half-integer topological defects paired via string micelles in polar liquids,

    Z. Ma, M. Jiang, Y. Song, A. Sun, S. Yi, C. Zhou, X. Huang, M. Huang, S. Aya, and Q.-H. Wei, “Half-integer topological defects paired via string micelles in polar liquids,” Proc. Natl. Acad. Sci. Nexus3(2024)

    work page 2024

  17. [17]

    Living liquid crystals,

    S. Zhou, A. Sokolov, O. D. Lavrentovich, and I. S. Aranson, “Living liquid crystals,” Proc. Natl. Acad. Sci. U. S. A.111, 1265–1270 (2014)

    work page 2014

  18. [18]

    Topological defects in a living nematic ensnare swimming bacteria,

    M. M. Genkin, A. Sokolov, O. D. Lavrentovich, and I. S. Aranson, “Topological defects in a living nematic ensnare swimming bacteria,” Phys. Rev. X.7(2017)

    work page 2017

  19. [19]

    Individual behavior and pairwise interactions between microswimmers in anisotropic liquid,

    A. Sokolov, S. Zhou, O. D. Lavrentovich, and I. S. Aranson, “Individual behavior and pairwise interactions between microswimmers in anisotropic liquid,” Phys. Rev. E91, 013009 (2015)

    work page 2015

  20. [20]

    Polar jets of swimming bacteria condensed by a patterned liquid crystal,

    T. Turiv, R. Koizumi, K. Thijssen, M. M. Genkin, H. Yu, C. Peng, Q.-H. Wei, J. M. Yeomans, I. S. Aranson, A. Doost- mohammadi, and O. D. Lavrentovich, “Polar jets of swimming bacteria condensed by a patterned liquid crystal,” Nat. Phys.16, 481–487 (2020)

    work page 2020

  21. [21]

    Biomechanical ordering of dense cell populations,

    D. Volfson, S. Cookson, J. Hasty, and L. S. Tsimring, “Biomechanical ordering of dense cell populations,” Proc. Natl. Acad. Sci.105, 15346–15351 (2008)

    work page 2008

  22. [22]

    Defect-mediated morphologies in growing cell colonies,

    A. Doostmohammadi, S. P. Thampi, and J. M. Yeomans, “Defect-mediated morphologies in growing cell colonies,” Phys. Rev. Lett.117, 048102 (2016)

    work page 2016

  23. [23]

    Bacteria solve the problem of crowding by moving slowly,

    O. J. Meacock, A. Doostmohammadi, K. R. Foster, J. M. Yeomans, and W. M. Durham, “Bacteria solve the problem of crowding by moving slowly,” Nat. Phys.17, 205–210 (2021)

    work page 2021

  24. [24]

    Individual bacterial cells can use spatial sensing of chemical gradients to direct chemotaxis on surfaces,

    J. H. R. Wheeler, K. R. Foster, and W. M. Durham, “Individual bacterial cells can use spatial sensing of chemical gradients to direct chemotaxis on surfaces,” Nat. Microbiol.9, 2308–2322 (2024)

    work page 2024

  25. [25]

    Local polar order controls mechanical stress and triggers layer formation in myxococcus xanthus colonies,

    E. Han, C. Fei, R. Alert, K. Copenhagen, M. D Koch, N. S. Wingreen, and J. W. Shaevitz, “Local polar order controls mechanical stress and triggers layer formation in myxococcus xanthus colonies,” Nat. Commun.16, 952 (2025)

    work page 2025

  26. [26]

    Topological defects in epithelia govern cell death and extrusion,

    T. Beng Saw, A. Doostmohammadi, V. Nier, L. Kocgozlu, S. Thampi, Y. Toyama, P. Marcq, C. T. Lim, J. M Yeomans, and B. Ladoux, “Topological defects in epithelia govern cell death and extrusion,” Nature544, 212–216 (2017)

    work page 2017

  27. [27]

    Topological defects control collective dynamics in neural progenitor cell cultures,

    K. Kawaguchi, R. Kageyama, and M. Sano, “Topological defects control collective dynamics in neural progenitor cell cultures,” Nature545, 327–331 (2017)

    work page 2017

  28. [28]

    Turbulent dynamics of epithelial cell cultures,

    C. Blanch-Mercader, V. Yashunsky, S. Garcia, G. Duclos, L. Giomi, and P. Silberzan, “Turbulent dynamics of epithelial cell cultures,” Phys. Rev. Lett.120, 208101 (2018)

    work page 2018

  29. [29]

    Topological excitations govern ordering kinetics in endothelial cell layers,

    I. Ruider, K. Thijssen, D. R. Vannier, V. Paloschi, A. Sciortino, A. Doostmohammadi, and A. Bausch, “Topological excitations govern ordering kinetics in endothelial cell layers,” (2024)

    work page 2024

  30. [30]

    Epithelia Realize Nematopolar Topological Defect Structures

    T. Ma, N. de Graaf Sousa, V. Grudtsyna, F. Vafa, and A. Doostmohammadi, “Epithelia realize nematopolar topological defect structures,” (2026), arXiv:2606.19844 [cond-mat.soft]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  31. [31]

    Generic theory of active polar gels: a paradigm for cytoskeletal dynamics,

    K. Kruse, J. F. Joanny, F. J¨ ulicher, J. Prost, and K. Sekimoto, “Generic theory of active polar gels: a paradigm for cytoskeletal dynamics,” Eur. Phys. J. E Soft Matter16, 5–16 (2005)

    work page 2005

  32. [32]

    Large-scale vortex lattice emerging from collectively moving microtubules,

    Y. Sumino, K. H. Nagai, Y. Shitaka, D. Tanaka, K. Yoshikawa, H. Chat´ e, and K. Oiwa, “Large-scale vortex lattice emerging from collectively moving microtubules,” Nature483, 448–452 (2012)

    work page 2012

  33. [33]

    Emergence of coexisting ordered states in active matter systems,

    L. Huber, R. Suzuki, T. Kr¨ uger, E. Frey, and A. R. Bausch, “Emergence of coexisting ordered states in active matter systems,” Science361, 255–258 (2018)

    work page 2018

  34. [34]

    Determinants of polar versus nematic organization in networks of dynamic microtubules and mitotic motors,

    J. Roostalu, J. Rickman, C. Thomas, F. N´ ed´ elec, and T. Surrey, “Determinants of polar versus nematic organization in networks of dynamic microtubules and mitotic motors,” Cell175, 796–808.e14 (2018)

    work page 2018

  35. [35]

    The interplay of polar and nematic order in active matter: implications for non-equilibrium physics and biology,

    V. Venkatesh, N. de Graaf Sousa, and A. Doostmohammadi, “The interplay of polar and nematic order in active matter: implications for non-equilibrium physics and biology,” J. Phys. A58, 263001 (2025). 17

    work page 2025

  36. [36]

    Hydrodynamics of self-propelled hard rods,

    A. Baskaran and M. C. Marchetti, “Hydrodynamics of self-propelled hard rods,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys.77, 011920 (2008)

    work page 2008

  37. [37]

    Nonlinear field equations for aligning self-propelled rods,

    A. Peshkov, I. S. Aranson, E. Bertin, H. Chat´ e, and F. Ginelli, “Nonlinear field equations for aligning self-propelled rods,” Phys. Rev. Lett.109, 268701 (2012)

    work page 2012

  38. [38]

    Hydrodynamics of soft active matter,

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, “Hydrodynamics of soft active matter,” Rev. Mod. Phys.85, 1143–1189 (2013)

    work page 2013

  39. [39]

    Phase separation and rotor self-assembly in active particle suspensions,

    J. Schwarz-Linek, C. Valeriani, A. Cacciuto, M. E. Cates, D. Marenduzzo, A. N. Morozov, and W. C. K. Poon, “Phase separation and rotor self-assembly in active particle suspensions,” Proc. Natl. Acad. Sci.Proc. Natl. Acad. Sci.109, 4052– 4057 (2012)

    work page 2012

  40. [40]

    Phase coexistence in two-dimensional passive and active dumbbell systems,

    L. F. Cugliandolo, P. Digregorio, G. Gonnella, and A. Suma, “Phase coexistence in two-dimensional passive and active dumbbell systems,” Phys. Rev. Lett.119, 268002 (2017)

    work page 2017

  41. [41]

    Alignment and propulsion of squirmer pusher ’ ¨A` ıpuller dumbbells,

    J. Clop´ es, G. Gompper, and R. G. Winkler, “Alignment and propulsion of squirmer pusher ’ ¨A` ıpuller dumbbells,” The Journal of Chemical Physics156, 194901 (2022)

    work page 2022

  42. [42]

    Strings in two-dimensional classical xy models,

    D. H. Lee and G. Grinstein, “Strings in two-dimensional classical xy models,” Phys. Rev. Lett.55, 541–544 (1985)

    work page 1985

  43. [43]

    Ordering kinetics and steady states of xy-model with ferromagnetic and nematic interaction,

    P. S. Mondal, P. K. Mishra, and S. Mishra, “Ordering kinetics and steady states of xy-model with ferromagnetic and nematic interaction,” J. Phys.36, 285101 (2024)

    work page 2024

  44. [44]

    Surface-directed dynamics in living liquid crystals,

    A. Vats, V. Banerjee, and S. Puri, “Surface-directed dynamics in living liquid crystals,” Phys. Rev. E.110, 034701 (2024)

    work page 2024

  45. [45]

    Phase diagram, confining strings, and a new universality class in nematopolar matter,

    F. Vafa and A. Doostmohammadi, “Phase diagram, confining strings, and a new universality class in nematopolar matter,” Europhys. Lett.152, 57002 (2025)

    work page 2025

  46. [46]

    String formation and arrested ordering kinetics in nematics induced by polar particles,

    P. K. Mishra, P. S. Mondal, P. Jena, and S. Mishra, “String formation and arrested ordering kinetics in nematics induced by polar particles,” New J. Phys.27, 074602 (2025)

    work page 2025

  47. [47]

    Active topological strings in renewing nematopolar fluids,

    A. Dinelli, L. Dumoulin, and K. Kruse, “Active topological strings in renewing nematopolar fluids,” (2026), arXiv:2601.18307 [cond-mat.soft]

    work page arXiv 2026

  48. [48]

    Unifying polar and nematic active matter: emergence and co-existence of half-integer and full-integer topological defects,

    D. Amiri, R. Mueller, and A. Doostmohammadi, “Unifying polar and nematic active matter: emergence and co-existence of half-integer and full-integer topological defects,” J. Phys. A55, 094002 (2022)

    work page 2022

  49. [49]

    Long-range order in a two-dimensional dynamical XY model: How birds fly together,

    J. Toner and Y. Tu, “Long-range order in a two-dimensional dynamical XY model: How birds fly together,” Phys. Rev. Lett.75, 4326–4329 (1995)

    work page 1995

  50. [50]

    R. J. LeVeque, Finite difference methods for ordinary and partial differential equations (SIAM, 2007)

    work page 2007

  51. [51]

    A numerical technique for predicting microstructure in liquid crystalline polymers,

    J. Hobdell and A. Windle, “A numerical technique for predicting microstructure in liquid crystalline polymers,” Liq. Cryst. 23, 157–173 (1997)

    work page 1997

  52. [52]

    Growth kinetics of systems with continuous symmetry,

    F. Liu and G. F. Mazenko, “Growth kinetics of systems with continuous symmetry,” Phys. Rev. B45, 6989–7001 (1992)

    work page 1992

  53. [53]

    Numerical method of lines for the relaxational dynamics of nematic liquid crystals,

    A. K. Bhattacharjee, Gautam I. Menon, and R. Adhikari, “Numerical method of lines for the relaxational dynamics of nematic liquid crystals,” Phys. Rev. E78, 026707 (2008)

    work page 2008

  54. [54]

    Scaling and vortex dynamics after the quench of a system with a continuous symmetry,

    M. Mondello and N. Goldenfeld, “Scaling and vortex dynamics after the quench of a system with a continuous symmetry,” Phys. Rev. A42, 5865–5872 (1990)

    work page 1990

  55. [55]

    Dynamical scaling in two-dimensional quenched uniaxial nematic liquid crystals,

    S. Dutta and S. K. Roy, “Dynamical scaling in two-dimensional quenched uniaxial nematic liquid crystals,” Phys. Rev. E 71, 026119 (2005)

    work page 2005

  56. [56]

    The topological theory of defects in ordered media,

    N. D. Mermin, “The topological theory of defects in ordered media,” Rev. Mod. Phys.51, 591–648 (1979)

    work page 1979

  57. [57]

    Kosterlitz ’ ¨A` ıthouless physics: a review of key issues,

    J. M. Kosterlitz, “Kosterlitz ’ ¨A` ıthouless physics: a review of key issues,” Rep. Prog. Phys.79, 026001 (2016)

    work page 2016

  58. [58]

    Topological point defects of liquid crystals in quasi-two-dimensional geometries,

    K. Harth and R. Stannarius, “Topological point defects of liquid crystals in quasi-two-dimensional geometries,” Front. Phys.Volume 8 - 2020(2020)

    work page 2020

  59. [59]

    Ordering, metastability and phase transitions in two-dimensional systems,

    J M Kosterlitz and D J Thouless, “Ordering, metastability and phase transitions in two-dimensional systems,” J. Phys. C 6, 1181 (1973)

    work page 1973

  60. [60]

    Patterning of morphogenetic anisotropy fields,

    Z. Wang, M C. Marchetti, and F. Brauns, “Patterning of morphogenetic anisotropy fields,” Proc. Natl. Acad. Sci.120, e2220167120 (2023)

    work page 2023

  61. [61]

    Free energy of a nonuniform system. I. Interfacial free energy,

    J. W. Cahn and J. E. Hilliard, “Free energy of a nonuniform system. I. Interfacial free energy,” J. Chem. Phys.28, 258–267 (1958)

    work page 1958

  62. [62]

    Numerical studies of phase separation in models of binary alloys and polymer blends,

    J. D. Gunton, R. Toral, and A. Chakrabarti, “Numerical studies of phase separation in models of binary alloys and polymer blends,” Physica Scripta1990, 12 (1990)

    work page 1990

  63. [63]

    Scalarϕ 4 field theory for active-particle phase separation,

    R. Wittkowski, A. Tiribocchi, J. Stenhammar, R. J Allen, D. Marenduzzo, and M. E. Cates, “Scalarϕ 4 field theory for active-particle phase separation,” Nat. Comm.5, 4351 (2014)

    work page 2014