Nematic director ($\mathbf{n}\sim -\mathbf{n}$)

A nematic director is a unit vector specifying the local unoriented axis of a nematic liquid crystal. Head-tail symmetry identifies $\mathbf{n}$ with $-\mathbf{n}$, so the order-parameter space is a Real projective space rather than a sphere.

For a unit nematic director, the bulk elastic density is $[K_1(\nabla \cdot \mathbf{n}{)}^2+K_2(\mathbf{n}\cdot \nabla \times \mathbf{n}{)}^2+K_3|\mathbf{n}\times \nabla \times \mathbf{n}{|}^2]/2$. The three positive moduli penalize nematic splay, nematic twist, and nematic bend. Boundary elastic terms, chirality, or varying nematic order require additional contributions.

For a laterally uniform nematic director $(\cos \theta \cos \varphi ,\cos \theta \sin \varphi ,\sin \theta )$, its divergence is $\cos \theta {\theta }^{\prime }$ and its twist invariant is $-{\cos }^2\theta {\varphi }^{\prime }$. Thus $f=K_1{\cos }^2\theta +K_3{\sin }^2\theta$ and $g=K_2{\cos }^4\theta +K_3{\sin }^2\theta {\cos }^2\theta$. A normal electric field contributes $-c{\sin }^2\theta$. The Euler-Lagrange equations for two fields are $(g{\varphi }^{\prime }{)}^{\prime }=0$ and $f{\theta }^{\prime \prime }+f_{\theta }{\theta }^{\prime 2}/2-g_{\theta }{\varphi }^{\prime 2}/2+2c\sin \theta \cos \theta =0$.

Nematic bend curves the integral lines of the nematic director. For a unit field, $(\mathbf{n}\cdot \nabla )\mathbf{n}=-\mathbf{n}\times (\nabla \times \mathbf{n})$, so the bend energy measures curvature along those lines.

Nematic twist rotates the nematic director around an axis perpendicular to it. The helix $\mathbf{n}=(\cos qz,\sin qz,0)$ has zero nematic splay and nematic bend, and twist density $K_2{q}^2/2$.

Nematic splay describes local spreading or convergence of the nematic director field. Its leading elastic cost is quadratic in the divergence and invariant under the head-tail identification $\mathbf{n}\sim -\mathbf{n}$.

A nematic disclination is a topological line defect of the unoriented director field. In a two-dimensional sample its point cross-section can have half-integer winding because $\mathbf{n}$ and $-\mathbf{n}$ represent the same state.

A neutral pair of opposite nematic disclinations can annihilate under dissipative relaxation, because the total topological charge of a two-dimensional nematic disclination is zero. Its far-field elastic energy increases as $2\pi K{\lambda }_0^2{q}^2\log (R/r_0)$ with separation $R$, so the pair has an attractive elastic force proportional to $-1/R$. As opposite cores meet, the nematic order parameter may vanish locally and the two windings disappear. Integer defects may first undergo energetic fission of an integer nematic disclination, leading to encounters of elementary half-charge defects. The actual paths and annihilation time require a dynamical law; the free energy alone does not specify core mobility. Pinning, imposed boundary winding or thermal defect production can obstruct a defect-free final state.

With constant amplitude outside a core, the planar nematic divergence-elasticity identity reduces the elastic density to $K{\lambda }_0^2|\nabla \chi {|}^2/2$. A topological charge of a two-dimensional nematic disclination $q$ has $\chi =q\theta +{\chi }_0$ and $|\nabla \chi {|}^2=q^2/r^2$. Integrating over an annulus gives the displayed energy, plus finite core and boundary terms. The outer cutoff is the sample size, compensating-defect separation, or another screening length. An isolated defect in an infinite ordered plane has divergent energy. A neutral pair $q,-q$ separated by $R$ instead has $E_{\mathrm{pair}}=2\pi K{\lambda }_0^2{q}^2\log (R/r_0)$ plus core terms, since its far-field gradients cancel.

In a strictly planar nematic liquid crystal, a charge-one nematic disclination is a winding-two map into the circle ${\mathbb{RP}}^1$. Its nonzero fundamental group class prevents complete unwinding while the surrounding loop stays fixed, but permits splitting into two charge-$1/2$ defects. For separation $r_0\ll d\ll L$, the two near fields contribute $\pi K{\lambda }_0^2(1/4+1/4)\log (d/r_0)$ while their common charge-one far field contributes $\pi K{\lambda }_0^2\log (L/d)$. Compared with the unsplit defect, the change isIt is negative at sufficiently large scale separation. Thus integer-defect fission is energetically favorable despite conservation of total topological charge. This leading-log argument proves failure of global energetic stability, not absence of kinetic barriers or pinning. If escape into the third dimension is allowed, an integer winding can also become topologically trivial in ${\mathbb{RP}}^2$.

For a planar tangential charge-one nematic disclination, $\nabla \cdot Q=-({\lambda }^{\prime }/2+\lambda /r){\mathbf{e}}_r$. This follows by differentiating the tensor dyad, using $\nabla \cdot {\mathbf{e}}_{\theta }=0$ and $({\mathbf{e}}_{\theta }\cdot \nabla ){\mathbf{e}}_{\theta }=-{\mathbf{e}}_r/r$. With bulk density $a{\lambda }^2/2+b{\lambda }^4/4$ and divergence elastic density $K({\lambda }^{\prime }/2+\lambda /r{)}^2/2$, radial variation of $2\pi \int rfdr$ givesFor $a<0$, a regular melted core has $\lambda =O(r^2)$ and an ordered far field has ${\lambda }_0=\sqrt{-a/b}$. Linearizing the bulk derivative gives $\lambda /{\lambda }_0=1-K/(2|a|r^2)+O(r^{-4})$. The core scale is of order $\sqrt{K/|a|}$; the amplitude disturbance has an integrable algebraic tail, rather than exactly compact support.

For a continuous director angle around a closed counterclockwise circuit, the enclosed topological charge is $q=\Delta \theta /(2\pi )$. Since $\theta$ is defined modulo $\pi$, elementary nematic disclinations can have $q=\pm 1/2$.

A two-dimensional director winding can sometimes unwind when the director may tilt out of the plane. For a three-dimensional nematic the fundamental group of the order-parameter space is ${\pi }_1({\mathbb{RP}}^2)=\mathbb{Z}/2\mathbb{Z}$, so two disclination lines are topologically equivalent to none.

In the one-elastic-constant approximation, slow director distortions have a quadratic gradient energy with one modulus. For a planar angle field it reduces to $F=(K/2)\int |\nabla \theta {|}^2$ up to the convention-dependent effective modulus.

For a planar nematic director $\mathbf{n}=(\cos \chi ,\sin \chi )$ and constant amplitude ${\lambda }_0$, the nematic order parameter isDifferentiation gives $\nabla \cdot Q={\lambda }_0(-\sin 2\chi {\chi }_x+\cos 2\chi {\chi }_y,\cos 2\chi {\chi }_x+\sin 2\chi {\chi }_y)$. The coefficient matrix is orthogonal, proving the identity. Hence $K|\nabla \cdot Q{|}^2/2$ is the angle-gradient energy with effective modulus $K{\lambda }_0^2$. The full tensor-gradient norm instead obeys ${\sum }_{ijk}({\partial }_k{Q}_{ij}{)}^2=2{\lambda }_0^2|\nabla \chi {|}^2$ in this normalization; the two local densities must not be confused.