Geometry and relaxation dynamics of nematic loops

Technical Summary: Geometry and Relaxation Dynamics of Nematic Loops

Problem Statement
While the topological classification of disclination lines in three-dimensional nematic liquid crystals is well-established (specifically via homotopy theory and the $Z_4$ classification for loops), this framework provides only a coarse-grained description. It successfully categorizes global topological features, such as whether loops are pierced or unpierced, but fails to encode the specific geometry of the defect profile along the loop. Recent work suggests that this local geometric structure—specifically the transport of the defect triad (or "flagpole") along the loop—contains physically relevant information that influences defect dynamics, a nuance missed by topology alone. The paper addresses the need for a mathematical framework that unifies the local defect profile with its global transport properties to predict dynamical outcomes.

Methodology
The authors employ a dual approach combining algebraic formulation and continuum numerical simulations:

1. Algebraic Formulation (Clifford Algebra):
   The authors propose describing nematic disclination loops using the three-dimensional Clifford algebra $Cl(3, 0)$. In this framework, the local defect profile (a triad describing the director rotation axis and in-plane phase) is encoded as a spinor in the even subalgebra $Cl^+(3, 0) \simeq SU(2)$.

   • Constant Profiles: They construct spinors for constant $+1/2$ and $-1/2$ wedge profiles.
   • Geometric Transformations: They utilize rotors within the algebra to model geometric variations. Specifically, they apply rotations around the local tangent ($T$), binormal ($z$), and normal ($N$) vectors.
     • Rotation around $T$ introduces self-twist (torsional transport).
     • Rotation around $z$ induces discrete transitions between $+1/2$ and $-1/2$ wedge profiles (e.g., Plus-Minus-Plus-Minus or PMPM).
     • Rotation around $N$ induces transitions between wedge and twist profiles.
   • Holonomy: The transport of the triad along the loop is mathematically represented as an $SU(2)$ holonomy, which the Clifford product exposes as a transparent composition of rotor factors.

2. Numerical Simulations:
   The authors perform continuum simulations of nematic liquid crystals using the Landau–de Gennes free energy and Beris–Edwards equations, solved via a hybrid lattice Boltzmann method.

   • Setup: A circular disclination loop with specific geometric profiles (defined by the Clifford spinors) is inserted into an otherwise uniform nematic.
   • Conditions: The system is initialized in the isotropic phase outside the loop torus. Simulations track the relaxation dynamics, defect fraction, and topological evolution over time.

Key Results
The simulations reveal that the geometric holonomy of the initial loop dictates the nucleation and topology of transient defect structures:

1. Topological Blob Formation: In all cases, the insertion of a loop nucleates a "topological blob" of defects. This blob grows rapidly (resembling turbulent blob dynamics in fluids) before dissipating as the system relaxes toward a defect-free nematic equilibrium. The defect fraction decays as $t^{-1/2}$ at late times.

2. Self-Twisted Loops (Torsional Holonomy):

   • Loops with a constant $-1/2$ wedge profile and a self-twist $\tau(\theta) = n\theta/3$ (where $n$ is an integer twist rate index) nucleate additional disclination lines that link with the original loop.
   • A simple arithmetic relation is established between the twist rate index $n$ and the number of linked loops $\ell$:
     $$\ell = \left\lfloor \frac{2n + 1}{3} \right\rfloor$$
   • The newly nucleated segments inside the loop predominantly exhibit a comet-like ($+1/2$) character, while the original loop retains its $-1/2$ trefoil character.

3. Discrete Profile Transitions (PMPM Loops):

   • Loops with an even number of discrete profile transitions (e.g., PMPM configurations with zero net self-twist) do not generate linked loops.
   • Instead, they nucleate defect structures that are "threaded" through the initial loop (interpenetrating but unlinked), analogous to threading in polymer melts.
   • These configurations contain both comet-like and trefoil-like segments within the enclosed region, indicating that the initial profile geometry sculpts the local character of the nucleated defects.

Significance and Claims
The paper claims to establish a direct, predictive link between the geometric $SU(2)$ holonomy of a disclination loop and its subsequent relaxation dynamics.

• Beyond Topology: The results demonstrate that geometric features are not merely mathematical refinements of topology but are physically predictive quantities with detectable signatures in the medium's dynamics.
• Algebraic Utility: The Clifford algebra formulation provides a systematic, unified language to decompose loop geometry into algebraically transparent components (self-twist vs. discrete transitions), allowing for the calculation of holonomy.
• Novelty: To the authors' knowledge, this is the first demonstration of a direct link between $SU(2)$ holonomy and specific defect nucleation dynamics (linking vs. threading).
• Scope: The authors suggest this framework may be extendable to more complex order parameter manifolds (cholesterics, smectics) and active media, though the current work is limited to passive nematics and circular loops. They explicitly note that the transient structures observed are not stabilized in these simulations, leaving the potential for stabilization via confinement or activity as a question for future inquiry.