Defect Kinematics in 2D Nematics: Contributions from Surface Topology, Intrinsic and Extrinsic Geometry, Solitons, Defect Orientations, and Elastic Anisotropy

Problem Statement
Topological defects in two-dimensional nematic materials are often described using a "particlelike" kinematic framework analogous to Coulomb gases in electromagnetism, where defects interact via logarithmic, pairwise forces. However, this analogy is incomplete. Unlike electromagnetism, which is a gauge-invariant theory, nematic materials are governed by elastic free energies that are not invariant under generic gauge transformations (local rotations of the director field). Consequently, quantities that are physically irrelevant in electromagnetism—such as global topology, intrinsic and extrinsic geometry, relative defect orientations, and elastic anisotropy—play crucial roles in nematic defect dynamics. The paper addresses the need for a unified theoretical framework that accounts for these geometric and topological factors, specifically investigating how they modify defect interactions, introduce many-body effects, and deviate from the standard pairwise Coulomb picture.

Methodology
The authors construct a geometric field theory for 2D nematics using the modern language of differential forms and exterior calculus. The core of the approach relies on the Hodge decomposition of the nematic connection 1-form $A$ (derived from the director angle $\theta$). The theory decomposes the field into three distinct contributions:

1. Defects: Singular fields encoding topological winding (charges).
2. Solitons: Smooth, topologically non-trivial harmonic fields (global excitations).
3. Orientation: Smooth, topologically trivial harmonic fields encoding relative defect orientations.

The authors utilize Hodge theory to treat solitons and orientation as "harmonic excitations" of the ground state. They analyze the free energy minimization under constraints imposed by these excitations. For systems where the energy functional remains linear (e.g., flat surfaces with specific topologies), they solve the Euler–Lagrange equations using Green's functions (propagators) of the Laplacian. For systems introducing non-linearities—specifically through extrinsic curvature and elastic anisotropy—the authors employ perturbation theory (expanding in a small parameter $\epsilon$) to derive corrections to the energy and forces, moving beyond the linear superposition principle.

Key Contributions and Results

• Geometric Field Theory Framework: The paper establishes that nematic defect dynamics can be rigorously described as a geometric field theory where the lack of gauge invariance allows for physical contributions from harmonic forms. This unifies the treatment of defects, solitons, and orientation under a single mathematical structure.
• Modification of the Propagator (Green's Function):
  • Global Topology: The interaction between defects is mediated by the Green's function of the Laplacian on the specific surface. The authors demonstrate how this function changes for surfaces with non-trivial topology (e.g., flat tori, cylinders, spheres), altering the pairwise interaction even in the absence of curvature.
  • Intrinsic Curvature: Intrinsic geometry modifies the Green's function via geodesic distances. The authors derive an expansion showing that curvature introduces corrections decaying as $\rho^2$ (where $\rho$ is the distance) with a prefactor proportional to the Gaussian curvature.
• Geometric Potential: Intrinsic curvature generates an effective "geometric potential" $V$ acting on individual defects, sourced by the Gaussian curvature $K$. This potential drives defects toward regions of matching curvature sign (positive charges to positive curvature, negative to negative), independent of defect-defect interactions.
• Harmonic Excitations (Solitons and Orientation):
  • Solitons: Smooth harmonic fields (associated with non-trivial cohomology classes, e.g., on a torus) do not modify the pairwise interaction Green's function but contribute a "self-energy" to the configuration. They can induce defect nucleation if their energy cost exceeds that of creating defects.
  • Defect Orientation: Relative orientations between defects are encoded as "spiral charges" (harmonic forms of the type $d \log \rho$). These introduce interaction terms proportional to the angle between defects ($\phi_{ij}$), analogous to an Aharonov-Bohm phase. The authors clarify how to consistently define these angles on curved surfaces using parallel transport or geodesic tangents.
• Non-Linearities and Many-Body Effects:
  • Extrinsic Curvature: When the director aligns with principal curvature directions, the energy includes a term proportional to curvature anisotropy ($\Delta \cos 2\theta$). This introduces non-linearity into the Euler–Lagrange equations. Perturbative analysis reveals that this induces many-body interactions (specifically torques depending on the sum of angles between all defects), breaking the purely pairwise nature of the interaction.
  • Elastic Anisotropy: Similarly, distinct splay and bend elastic constants introduce non-linearities. The authors show this leads to distance-independent torques and many-body interactions that depend on the relative orientations and positions of all defects in the system.

Significance and Claims
The paper claims that this geometric field theory approach provides a "clarifying perspective" on previously reported effects in the literature, unifying them under the concepts of propagator modification and harmonic excitations.

• Conceptual Unification: It demonstrates that the deviation from electromagnetism in nematics is not merely a collection of ad-hoc corrections but a structural consequence of gauge non-invariance, captured systematically by Hodge theory.
• Many-Body Interactions: A primary finding is that despite the underlying $U(1)$ symmetry being abelian, the lack of gauge invariance allows for non-linearities (via extrinsic curvature and elastic anisotropy) that generate genuine many-body interactions. This subverts the expectation that abelian systems only exhibit pairwise forces.
• Practical Utility: The authors argue that their framework simplifies calculations by reducing complex field problems to linear Poisson-type problems (solvable via Green's functions) for the leading order, with non-linear effects treated perturbatively. They suggest this approach is well-suited for extension to 3D nematics (by changing the symmetry group to $SU(2)$) and for numerical studies of defect dynamics on deformable surfaces.

The paper concludes that while the framework offers analytical insight and a structured way to handle complex geometries, the resulting many-body alignment problems are complex and likely require numerical analysis for full resolution, particularly in the context of active nematics or deformable biological membranes.