Ordering and Defect Dynamics in Passive and Active Nematopolars

Technical Summary: Ordering and Defect Dynamics in Passive and Active Nematopolars

Problem Statement
The paper addresses the theoretical understanding of non-equilibrium collective behavior in "nematopolar" systems, where polar (head-tail asymmetric) and nematic (head-tail symmetric) interactions coexist. Such systems are observed in diverse biological and synthetic contexts, including ferroelectric liquid crystals, living liquid crystals, bacterial colonies, and microtubule-motor mixtures. While previous theoretical approaches have utilized two-field models (coupling a polar vector to a nematic tensor) or focused on active turbulence, there is a lack of systematic characterization regarding phase ordering and coarsening dynamics in minimal single-field models for dry nematopolar systems. Specifically, the interplay between competing symmetries and the role of self-propulsion (activity) in defect dynamics and domain growth laws remain unexplored.

Methodology
The authors propose and numerically investigate a minimal single-field "dry" model for nematopolar systems.

• Model: The system is described by a single polar vector field $\mathbf{p}$ evolving via purely relaxational dynamics (Model A-like) in the passive limit. The free energy functional includes a double-well potential for bulk polarization, a polar elastic term ($\propto |\nabla \mathbf{p}|^2$), and a nematic-like term derived from a tensor constructed from $\mathbf{p}$ ($\hat{P} \propto \mathbf{p}\mathbf{p}^\top - |\mathbf{p}|^2 I/2$).
• Activity: Non-equilibrium behavior is introduced via a self-advection term, $\Lambda(\mathbf{p} \cdot \nabla)\mathbf{p}$, in the equation of motion. This term allows the polarization to act simultaneously as an order parameter and a local velocity field, consistent with Toner-Tu theories for dry active matter.
• Simulation: The authors employ finite-difference schemes on 2D square lattices ($N=512, 1024$) with periodic boundary conditions.
• Observables: The study tracks the location, number, and type of topological defects (integer and half-integer charges), polar and nematic correlation functions ($C_p, C_n$), and characteristic correlation lengths ($L_p, L_n$).
• Approach: The investigation proceeds in two stages: (1) controlled simulations of isolated defect pairs and loops to analyze elementary relaxation mechanisms, and (2) large-scale simulations starting from disordered initial conditions to study global coarsening dynamics.

Key Contributions and Results

1. Defect Morphology and String-Mediated Interactions (Passive Case):

   • In the passive regime ($\Lambda=0$) with a balanced competition between polar and nematic couplings, the system spontaneously develops depolarization strings (lines of vanishing polarization) connecting half-integer defects ($\pm 1/2$) and closed depolarization loops.
   • Defect Pairs: Oppositely charged half-integer defects connected by a string exhibit non-monotonic annihilation velocities as a function of the relative polar/nematic strength ratio ($\bar{k} = k_p/k_n$). This arises from a competition between the reduction of effective elastic constants and the energy cost of the string itself.
   • Like-Charged Defects: Unlike in pure nematic systems where like charges repel indefinitely, like-charged half-integer defects in this model reach a finite equilibrium separation determined by the balance between repulsive Coulomb-like forces and the tension of the connecting string.
   • Loop Dynamics: Depolarization loops relax via two distinct mechanisms depending on system parameters: (i) shrinking (curvature-driven, dominant in nematic-dominated regimes) and (ii) rotational relaxation where the polarization field inside the loop reorients to match the outside (dominant in polar-dominated regimes). In extreme polar regimes, loops can rupture into integer-charge strings.

2. Coarsening Dynamics (Passive Case):

   • Starting from disordered states, the system undergoes coarsening. The characteristic length scale $L(t)$ grows as $L(t) \sim (t/\ln t)^{1/2}$.
   • This scaling is consistent with systems possessing continuous rotational symmetry and point-like defects (analogous to the 2D XY model and O(2)/Z2 nematic systems), confirming that the presence of strings and loops does not alter the fundamental dynamic scaling law, though it slows the absolute rate of ordering compared to pure nematic systems.
   • The total number of defects decays as $n_{tot} \sim \ln t / t$.

3. Active Dynamics and Arrested Coarsening:

   • Upon introducing self-advection ($\Lambda \neq 0$), the system exhibits Motility-Induced Charge Symmetry Breaking (MICS). The symmetry between positive and negative defects is broken, leading to the coexistence of positive integer defects (aster-like) and negative half-integer defects (trefoil-like).
   • Negative integer defects and positive half-integer defects are dynamically suppressed and decompose into the stable species.
   • Arrested Coarsening: Sufficiently strong activity leads to a stationary state where the characteristic length scale saturates (plateaus) and the defect density remains finite. The system enters an "arrested coarsening" regime where defect evolution is strongly suppressed, preventing the system from reaching a uniform ordered state.

Significance
The paper claims to provide a simple unified framework for understanding ordering and defect dynamics in nematopolar systems. By utilizing a minimal single-field model, the authors demonstrate that:

• Complex defect structures (strings, loops, and mixed-charge coexistence) can emerge from the competition between polar and nematic free energy terms without requiring multiple order parameters.
• Activity, implemented via a simple self-advection term, is a fundamental mechanism capable of reshaping defect-mediated ordering, stabilizing non-equilibrium states (arrested coarsening), and inducing symmetry breaking that has no counterpart in passive systems.
• The results offer a theoretical basis for interpreting the behavior of various biological and synthetic systems, such as ferroelectric liquid crystals, living liquid crystals, and microtubule-motor mixtures, where both types of order emerge from the same degrees of freedom.

The authors explicitly note that while their model is minimal, it captures the essential phenomenology of these systems and disentangles the roles of competing symmetries and activity. They suggest that this framework serves as a starting point for investigating broader classes of active materials with competing orientational orders.