Defect order in active nematics on a curved surface

We investigate the effects of extrinsic curvature on the turbulent behavior of a 2D active nematic confined to the surface of a cylinder. The surface of a cylinder has no intrinsic curvature and only extrinsic curvature. A nematic field reacts to the extrinsic curvature by trying to align with the lowest principle curvature, in this case parallel to the long axis of the cylinder. When nematics are sufficiently active, there is a proliferation of defects arising from a bend or splay instability depending on the nature of the active stress. The extrinsic curvature of the cylinder breaks the rotational symmetry of this process, implying that defects are created parallel or perpendicular to the cylinder depending on whether the active nematic is contractile or extensile.

Active matter describes a class of non-equilibrium materials in which stress is generated at the subunit scale; for example by converting some chemical energy source into directed motion [1,2]. Active nematics are active materials which additionally have a broken rotation symmetry, resulting in nematic liquid crystalline phases [2][3][4][5][6]. Active nematics have been observed as a naturally occurring phenomena in the motion of tissues [7,8] and bacteria [9][10][11][12], but can also be experimentally realized from reconstituted biological components; most commonly microtubule kinesin suspensions [3,4,13,14]. These active nematics show a huge range of phenomenological behaviors dependent on the level to which they are driven [15][16][17], the geometry of the enclosure [9,18], the properties of the substrate [13,19,20], the density [11] or the boundary conditions [21]. The novel features of active matter make it of particular interest in the design of new materials with non-trivial properties [22,23]. This often requires that some control is exercised over the active matter which has been done previously through geometry and topology [22][23][24][25][26][27].
The interplay between active nematics and substrate topology were first studied experimentally by Keber et al, who put active nematics on the surface of a vesicle [25]. This lead to a flurry of interesting results focusing on active nematics on curved surfaces [28,29] and in particular spheres [25,[30][31][32][33]. Using 3D printed water droplets it is possible to study the relationship between active nematics and surface Gaussian curvature on more complex surfaces [26], which has the effect of sorting topological defects into regions of like-sign Gaussian curvature [26,27].
In this paper we study the effects of extrinsic curvature on topological defects within active nematics. This has been included in previous works, showing how it affects the lowest energy state for passive nematics [34][35][36], the laminar flow regime in active nematics [37], and in generalized active nematics for curved surfaces [27]. However its precise role in turbulent active nematics is not yet known. In order to isolate the effects of extrinsic curvature we study active nematics on the surface of a periodic cylinder; a developable surface with no Gaussian curvature. This also has the topological constraint of zero net topological charge within the nematic. We start by introducing the equations for an active nematic on a cylindrical surface, which requires only the addition of a single term to the free energy. We then go on to show how this affects the general behavior of an active nematic including the emergence of global nematic defect order. We show that this is an active effect by observing no such defect order for passive nematics on a cylinder. We then go on to propose a mechanism by which defect order is generated which we confirm by observing contractile active nematics.
Vectors on the surface of a cylinder can be written in cylindrical coordinates in the form n = n θê θ + n zê z , whereê θ andê z are basis vectors in the vertical and azimuthal directions, respectively. If we set our length-scale by the radius of the cylinder, this basis naturally becomes orthonormal, hence the surface metric can be written g ij = δ ij . In addition, the basis vectors are independent of each other, hence the Christoffel symbols are all zero, Γ k ij = 0. Finally, the surface of a cylinder is developable, therefore has no Gaussian curvature. This means that when writing the dynamical equations for a fluid on the surface of a cylinder, the intrinsic geometry of the surface has no effect on their form. Therefore we start from the generic form of the equations governing an in-compressible active nematic which are given by: Since we are in a cylindrical geometry i and j now run through z and θ, the vertical and azimuthal directions, respectively. Here ρ is the density of the nematogens, σ (t) ij is the total stress tensor and the μ is the friction coefficient. Since we are considering the in-compressible limit, ∂ i v i = 0 and we set ρ = 1 everywhere. Q ij = S(n i n j − δ ij /2) is the nematic tensor, combining the nematic order parameter, S, and a unit vector pointing parallel to the average nematic director,n = cos(ζ)ê θ + sin(ζ)ê z ; here ζ is the angle between the average nematic orientation and theθ direction. λ is the flow alignment parameter coupling the evolution of the nematic tensor to the strain rate tensor, u ij = (∂ i v j + ∂ j v i )/2. The vorticity tensor is given by ω ij = (∂ i v j − ∂ j v i )/2 and the molecular tensor is given by H ij = −∂F/∂Q ij where F is a Landau de Gennes free energy modified in order to account for the extrinsic curvature of the cylinder.
Where the parameter is a characteristic length which is proportional to the core defect radius and K is the elastic constant associated with distortions in the director field. f e describes the energy cost associated with coupling the nematic tensor to the extrinsic curvature. The extrinsic curvature tensor is given by C ij = −g i · ∂ jN ; where g i are the basis vectors andN is the unit surface normal; on the surface of a cylinder this only has one non-zero component C θθ = 1/r, where r is the radius of the cylinder. The extrinsic curvature coupling energy is written f e = K e (QC 2 ) which can be simplified to f e = K e Q θθ /r 2 on a cylinder [34]. Since we have set our length scale by the radius of the cylinder, we consider only changes in K e ; this captures the same effect as a change in curvature. This energy describes the fact that the director embedded in the surface still has to follow the curvature of the surface as it is embedded in 3D space. This has the net effect of aligning the director with the central axis of the cylinder. If the director is aligned in the azimuthal direction, the director is curved within the embedding space since the surface of the cylinder is curved in that direction; this comes with an associated cost in elastic energy. However if the director is aligned in the vertical direction it experiences no such curvature, and hence has no energy penalty. It should be noted that if the fibers within the nematic had non-zero spontaneous curvature the effect would be reversed and the fibers would align with the cylinder at an angle which accommodates their curvature; a mechanism that has been proposed for the alignment of FtsZ fibers at bacterial division sites [38][39][40].
Here α is the magnitude of the active stress and ν is the viscosity coefficient; here we only include isotropic viscosity. The active stress is what separates this system from a passive liquid crystal. This describes the generation of stresses at the microscopic length scale which can either be extensile, α < 0, or contractile, α > 0. When the magnitude of the activity is sufficiently high, the additional active stress can overcome the elastic stresses in the system and there is a proliferation of topological defects and associated flow in a state known as 'active turbulence'. The lowest order defects in a 2D nematic system have a half integer charge associated with a half turn of the director as a line is traced around the core of the defect; this results in the characteristic ±1/2 nematic defects. The orientation of the director around the core is given by ζ = k(φ − ψ) + ψ where k is the charge of the defect, φ is the angle between a reference axis, in this case thê z direction, and the position around the defect core [41]. ψ is the orientation of the defect which is polar in the case of +1/2 defects and triatic in the case of −1/2 defects.
Equations (1) and (2) are simulated with a periodic boundary in two dimensions recreating an active nematic on the surface of a periodic cylinder, see figures 1(a) and (b) and supplementary movie (https://stacks.iop.org/NJP/22/063051/mmedia). The model parameters are selected such that the system is in a state of active turbulence containing of the order of 100 defects, see SI for details.
As the extrinsic curvature coupling coefficient, K e , is increased we see an increase in the average value of Q zz , which implies a decrease in the average value of Q θθ , see figure 1(c). This suggests that much of the active nematic shows net alignment with the cylinder, even in the active turbulence regime. This effect is strongest when the magnitude of the active stress is smallest. This is because the active stress acts to destabilize the ordered nematic state and introduce defects. Defects necessitate variations in the nematic director so it inherently becomes less ordered. The value of Q zθ is uniformly low, meaning that there is no detected alignment in any other direction, see figure 1(c). The defects in an active nematic on flat space show no global ordering [19], however when K e is increased we see an emergent nematic order in which the defects on average align themselves pointing parallel or anti-parallel to the azimuthal direction, see figure 1(d).
The density of defects in the active turbulent state depends on the activity and elastic coefficient; defects can be generated on the length-scale at which the active stress balances the elastic energy, [15,17]. This gives a defect density that scales ρ D ∼ α/K, hence increasing the activity increases the number of defects. The extrinsic curvature has a small effect on the number of defects (see SI) and there is only a slight effect of the activity on the ordering of defects on a cylinder, see figure 1(e).
We can observe a passive nematic by removing the active stress, α = 0, see figures 2(a) and (b). In this case the system quickly approaches the minimum energy configuration with a uniform director aligned with the cylinder, as K e is increased the speed at which the system equilibrates is faster, see figure 2(c). The system is initialized with a randomized director which contains many defects which annihilate over the course of the simulation. This coarse-graining of the defects also processes faster when K e is increased, see figure 2(d). This is because the extrinsic curvature breaks the rotational symmetry of the system, hence there is no spontaneous symmetry breaking required. This speeds up the coarse-graining process, see figure 2(d). The defects in the passive liquid crystal are analysed when the system has reached 400 defects. This process is repeated 500 times with random starting conditions and the results are pooled to construct a probability density for the defect orientation, see figure 2(e). Here is is clear to see that there is no emergent defect order in a relaxing passive nematic on the surface of a cylinder.
The passive system is entirely driven by relaxation of the elastic energy. The elastic energy associated with a defect is described by the modified Landau-De Gennes free energy given above, see equation (3). The only part of this energy which breaks rotational symmetry with respect to the cylinder is the extrinsic curvature energy f e . Therefore we can write the elastic torque on a defect core to be τ = ∂ ψ Δ 2 dAf = ∂ ψ Δ 2 dAf e , where Δ 2 is a small neighborhood around the core of a defect and f represents the elastic energy density. By substituting in the expression for the director around the core of a defect it is easy to arrive at the following result.
Hence there are no elastic torques on the defects which explains why we observe no orientational ordering of defects within the passive system. Therefore the ordering we observe in the active system must come from the active stress.
The principal difference between the passive liquid crystals and active nematics in the turbulent regime is the generation of topological defects. When the active stress is large, the ordered nematic state becomes unstable and we see defect proliferation and the so called 'active turbulence' regime [17]. The nature of this instability depends primarily on the sign of the active stress, α [15,16]. Simple periodic deformations of the nematic director can be written in the form Δζ = cos(2πk i r i ), where k is the wave vector of the deformation and r is the position within the nematic texture. For an extensile system, α < 0, the system is unstable to bend deformations; this is when k is aligned with the nematic director (k i = kn i ), see figure 3(a). When this perturbation becomes unstable it leads to the generation of a pair of defects aligned perpendicular to the original nematic field, see figure 3(a). When K e is increased, the nematic aligns with theẑ direction, hence k i = kê z and the bend deformations will be of the form Δζ = cos(2πkz), leading to the generation of a defect pair aligned in the ±θ direction. This is confirmed by looking at the orientation of defects immediately after they are created, see figure 3(b). The stability of the ordered nematic also depends on the value of the flow alignment parameter, λ [15][16][17]20]. When λ is increased the bend deformation is no longer the only unstable mode that can lead to defect proliferation in extensile active nematics, hence we observe a decrease in the global defect order, see figure 3(c). This mechanism relies on defect generation, explaining why increasing the activity beyond the defect proliferation threshold only slightly affects defect order non-monotonically, see figure 1(e). Increasing activity decreases the lifetime of a defect, giving it less time to reorientate after generation, hence defect order increases. Further increasing the activity slightly decreases the effect since an increased defect density leads to less order in the nematic director hence less regions aligned with theẑ direction, see figure 1(c). Therefore new defects are less likely to be generated aligned with the ±θ direction.
This hypothesis relies on the fact that extensile active nematics are predominantly unstable to bend deformations. If the sign of the active stress is changed, α > 0, the active nematic becomes contractile in nature and the system is now predominantly unstable to splay deformations [15,16]. This is when the wave vector for the perturbation is perpendicular to the nematic director (k i = kn j ij ). When activity is sufficiently high, the splay deformation is unstable and leads to the generation of a pair of defects oriented parallel to the initial nematic director, see figure 4(a). As before, the extrinsic curvature aligns the nematic director with theẑ direction, meaning the wave vector is now pointing in theθ direction and splay deformations are of the form Δζ = cos(2πkθ). Therefore, when the system is contractile, α > 0, we see that the extrinsic curvature coupling causes an increased number of defects aligned perpendicular to theẑ direction, see figure 4(b). Again, this anisotropy is observed in the defects at the moment of generation, see SI. In the contractile case the response to the flow alignment parameter is reversed and the effect is accentuated when λ is increased, see figure 4(c). This is because when λ is decreased, the bend deformation becomes relatively less stable and defects are generated at other orientations. We observe a π/2 rotation of the nematic defect order depending on the sign of the activity, with contractile systems having defects predominantly aligned with theẑ direction and extensile systems having defects predominantly aligned with theθ direction.
By placing an active nematic on a simple curved surface it is possible to generate nematic defect order, even though the curvature does not exert a torque directly on the core of the defect. The curved surface causes the nematic to align with the smallest principle curvature, resulting in a preferred orientation for the director field. Defects are generated in pairs aligned perpendicular or parallel to the nematic if they come from bend or splay deformations, respectively. The instability of bend or splay deformations depends primarily on the nature of the active stress, with extensile active nematics being unstable to bend deformations and contractile active nematics being unstable to splay deformations. The flow alignment parameter affects the stability of the bend and splay modes differently and can be used to accentuate this phenomenon. We speculate that if this analysis was expanded to include multiple nematic elastic coefficients, the relative bend and splay elastic coefficients would also affect this transition [42].
The mechanism outlined here can be extended to more general geometries such as a torus [26,27] or textured sheets, and can be used to predict how active nematics would behave on such manifolds. In addition it can be used to measure the extrinsic curvature coupling coefficient of experimentally realized active nematics. On the surface of a cylinder the extrinsic curvature coupling is analogous to alignment to an external field or anisotropic substrate. Indeed a global defect ordered state has been observed experimentally in active nematics on an anisotropic substrate and may be explained by the effects described here combined with those caused by anisotropic hydrodynamics [13,19,20]. Finally, extrinsic curvature coupling could play an important role in the organization of the cytoskeleton. This is the network of biopolymers that both gives the cell many of its mechanical properties and is able to generation motion of the cell and its components. It has already been suggested that extrinsic curvature coupling plays an important role in the location of the FtsZ division machinery at the center of a dividing bacterial cell [38][39][40].