A computational study of nematic core structure and disclination interactions in elastically anisotropic nematics.

A singular potential method in the Q tensor order parameter representation of a nematic liquid crystal is used to study the equilibrium configuration of a disclination dipole. Unlike the well studied isotropic limit (the so called one constant approximation), we focus on the case of anisotropic Frank elasticity (bend/splay elastic constant contrast). Prior research has established that the singular potential method provides an accurate description of the tensor order parameter profile in the vicinity of a disclination core of a highly anisotropic lyotropic chromonic liquid crystal. This research is extended here to two interacting disclinations forming a dipole configuration. The director angle is shown to decay in the far field inversely with distance to the dipole as is the case in the isotropic limit, but with a different angular dependence. Therefore elastic constant anisotropy modifies the elastic screening between disclinations, with implications for the study of ensembles of defects as seen, for example, in active matter in the extended system limit.

Spontaneous flows and dynamics of full-integer topological defects in polar active matter.

Polar active matter of self-propelled particles sustain spontaneous flows through the full-integer topological defects. We study theoretically the incompressible flow profiles around ±1 defects induced by polar and dipolar active forces. We show that dipolar forces induce vortical flows around the +1 defect, while the flow around the -1 defect has an 8-fold rotational symmetry. The vortical flow changes its chirality near the +1 defect core in the absence of the friction with a substrate. We show analytically that the flow induced by polar active forces is vortical near the +1 defect and is 4-fold symmetric near the -1 defect, while it becomes uniform in the far-field. For a pair of oppositely charged defects, this polar flow contributes to a mutual interaction force that depends only on the orientation of the defect pair relative to the background polarization, and that enhances defect pair annihilation. This is in contradiction with the effect of dipolar active forces which decay inversely proportional with the defect separation distance. As such, our analyses reveals a long-ranged mechanism for the pairwise interaction between topological defects in polar active matter.

Defect self-propulsion in active nematic films with spatially varying activity.

We study the dynamics of topological defects in active nematic films with spatially varying activity and consider two set-ups: (i) a constant activity gradient and (ii) a sharp jump in activity. A constant gradient of extensile (contractile) activity endows the comet-like +1/2 defect with a finite vorticity that drives the defect to align its nose in the direction of decreasing (increasing) gradient. A constant gradient does not, however, affect the known self-propulsion of the +1/2 defect and has no effect on the -1/2 that remains a non-motile particle. A sharp jump in activity acts like a wall that traps the defects, affecting the translational and rotational motion of both charges. The +1/2 defect slows down as it approaches the interface and the net vorticity tends to reorient the defect polarization so that it becomes perpendicular to the interface. The -1/2 defect acquires a self-propulsion towards the activity interface, while the vorticity-induced active torque tends to align the defect to a preferred orientation. This effective attraction of the negative defects to the wall is consistent with the observation of an accumulation of negative topological charge at both active/passive interfaces and physical boundaries.

Flow around topological defects in active nematic films.

We study the active flow around isolated defects and the self-propulsion velocity of + 1 / 2 defects in an active nematic film with both viscous dissipation (with viscosity η ) and frictional damping Γ with a substrate. The interplay between these two dissipation mechanisms is controlled by the hydrodynamic dissipation length ℓ d = η / Γ that screens the flows. For an isolated defect, in the absence of screening from other defects, the size of the shear vorticity around the defect is controlled by the system size R . In the presence of friction that leads to a finite value of ℓ d , the vorticity field decays to zero on the lengthscales larger than ℓ d . We show that the self-propulsion velocity of + 1 / 2 defects grows with R in small systems where R < ℓ d , while in the infinite system limit or when R ≫ ℓ d , it approaches a constant value determined by ℓ d .