2. Mathematical formulation

Consider a two-dimensional nematic liquid crystal outside a simply connected body, we denote this domain $D$ with boundary $\partial D$, as sketched in figure 1. Assuming the one-constant approximation, the director angle, $\theta (x,y)$, is described by the bulk free energy in (1.2). At the boundaries, we assume the Rapini–Papoular form of the surface anchoring energy, $\mathcal {F}_{surface} := W\sin ^2(\theta -\phi )/2$, where $W$ is the anchoring strength and $\phi$ is the preferred orientation, which is generally a function of the local boundary orientation (Rapini & Papoular Reference Rapini and Papoular1969). For tangential anchoring, $\phi$ is the surface tangent angle. Combining the bulk and surface energies, the total energy of the configuration is defined as

(2.1)\begin{equation} \mathcal{E} := \frac{K}{2}\iint_{\mathcal{D}}|\boldsymbol{\nabla} \theta|^2\,\mathrm{d} A + \frac{W}{2}\int_{\partial D}\sin^2(\theta-\phi)\,\mathrm{d} s, \end{equation}

where $\mathrm {d} A$ and $\mathrm {d} s$ are the infinitesimal surface area and arclength elements, respectively.

Figure 1.

Sketch of the domain considered in this paper: a rigid body immersed in a two-dimensional nematic liquid crystal. The liquid crystal is described by a director field $\boldsymbol {n} = (\cos \theta, \sin \theta, 0)$ with director angle $\theta (z)\in [0,{\rm \pi} )$ for $z\in D$. The boundary of the body is shown as a solid curve, $\partial D$, with unit normal and tangent vectors $\boldsymbol {\hat {\nu }}(s)$ and $\boldsymbol {\hat {s}}(s)$, respectively. The effective boundary (used in the effective boundary technique, § 3) is shown as a dashed curve, $\partial D_w$.

A variational principle applied to (2.1) yields the equilibrium equations for the director angle, $\theta (x,y)$ (see Appendix A),

(2.2)\begin{equation} \nabla^2\theta = 0\quad \text{in }D, \end{equation}

subject to the weak anchoring boundary condition,

(2.3)\begin{equation} -K\frac{\partial \theta}{\partial\hat{\nu}} + \frac{W}{2}\sin[2(\theta-\phi)] =0\quad \text{on }\partial D, \end{equation}

where $\boldsymbol {\hat {\nu }}:= -\boldsymbol {x}_s^\perp$ is the liquid crystal pointing unit normal, as depicted in figure 1.

One of the goals of this work is to introduce a methodology for computing the force exerted on a body immersed in a liquid crystal. The liquid crystal imposes stress according to the Ericksen stress tensor $\sigma := -p I - K\boldsymbol {\nabla } \boldsymbol {n}^{\rm T}\boldsymbol {\nabla } \boldsymbol {n}$ with hydrostatic pressure $p:= -K\lVert \boldsymbol {\nabla } \boldsymbol {n}\rVert ^2/2$ (de Gennes & Prost Reference de Gennes and Prost1993; Stewart Reference Stewart2004). On $\partial D$, there is a additional stress, which originates from the weak anchoring surface energy (Virga Reference Virga2018). A virtual work argument yields the total surface traction acting on the boundary $\partial D$ (see Appendix A),

(2.4)\begin{align} \boldsymbol{t} &= \mathcal{F}_{bulk}\boldsymbol{\hat{\nu}}-\left(\boldsymbol{\hat{\nu}}\boldsymbol{\cdot}\frac{\partial \mathcal{F}_{bulk}}{\partial\boldsymbol{\nabla}\theta} \right)\boldsymbol{\nabla}\theta + \left(\mathcal{F}_{surface}\boldsymbol{\hat{s}}-\frac{\partial \mathcal{F}_{surface}}{\partial \phi} \boldsymbol{\hat{\nu}} \right)_{s}\nonumber\\ &=K\left(\frac{1}{2}|\boldsymbol{\nabla}\theta|^2\boldsymbol{\hat{\nu}}-\frac{\partial \theta}{\partial\hat{\nu}} \boldsymbol{\nabla}\theta\right)+ \frac{W}{2}\{\sin(\theta-\phi)^2 \boldsymbol{\hat{s}}+\sin[2(\theta-\phi)]\boldsymbol{\hat{\nu}} \}_{s}, \end{align}

where $\boldsymbol {\hat {\nu }}:= -\boldsymbol {x}_s^\perp$ and $\boldsymbol {\hat {s}}:= \boldsymbol {x}_s$ are the unit normal and tangent vectors, respectively, and subscript $s$ denotes an arclength derivative, which is defined anticlockwise.

Overall, given a harmonic director field that satisfies (2.3), the energy and surface traction associated with the liquid crystal can be computed using (2.1) and (2.4), respectively. In this paper, we utilize a complex analysis formulation to derive analytical solutions to these equations, we shall introduce this formulation shortly.

2.1. Dimensionless variables

We non-dimensionalize all length scales with respect to a characteristic length scale associated with the immersed body, $a$, and introduce the dimensionless free energy and traction as $\hat {\mathcal {E}}:= \mathcal {E}/K$ and $\hat {\boldsymbol {t}}:= a^2 \boldsymbol {t}/K$, respectively. The resulting equations are governed by a dimensionless anchoring strength $w := aW/K$. The dimensionless free energy of the liquid crystal, $\hat {\mathcal {E}}$, may now be written as a boundary integral using the divergence theorem,

(2.5)\begin{equation} \hat{\mathcal{E}} = \frac{1}{2}\int_{\partial D}-\theta \frac{\partial \theta}{\partial\hat{\nu}}+w\sin^2(\theta-\phi)\, \mathrm{d} s. \end{equation}

Henceforth, we shall only work in these dimensionless variables.

2.2. Complex variables formulation

We introduce the complex coordinate $z:= x+\mathrm {i} y$ and define the complex director angle as

(2.6)\begin{equation} \varOmega(z) := \tau(x,y) - \mathrm{i} \theta(x,y), \end{equation}

where $\tau (x,y)=\operatorname {Re}\varOmega (z)$ is a harmonic conjugate of $\theta (x,y)=-\operatorname {Im}\varOmega (z)$ (i.e. $\theta _x=\tau _y$ and $\theta _y=-\tau _x$). The gradient components of the director angle, $\boldsymbol {\nabla }\theta =(\theta _x,\theta _y)$, are related to the complex angle above via

(2.7)\begin{equation} \theta_x - \mathrm{i} \theta_y = \mathrm{i} \varOmega'(z). \end{equation}

Since $\theta (x,y)$ is harmonic in $D$, $\theta _x-\mathrm {i}\theta _y$ must be holomorphic in $D$ and $\varOmega (z)$ is at least locally holomorphic. In general, $\varOmega$ may not be single-valued around the body, thus the period around $\partial D$ must be defined, that is

(2.8)\begin{equation} \oint_{\partial D} \mathrm{d} \varOmega \equiv \frac{1}{\mathrm{i}}\int_{\partial D}\theta_x -\mathrm{i} \theta_y \mathrm{d} z = \varUpsilon-2{\rm \pi}\mathrm{i} M, \end{equation}

for some given real constant $\varUpsilon$ and half-integer $M$. Note that $M$ corresponds to the topological charge of the body, $\partial D$.

In these complex variables, the bulk pointing normal vector, $\boldsymbol {\hat {\nu }}$, is represented as $y_s-\mathrm {i} x_s$ and the normal derivative of the director angle is written as $\theta _{\nu }=\operatorname {Re}[\varOmega '(z)z_s]$. The boundary condition (2.3) is equivalent to a constraint on $\varOmega (z)$,

(2.9)\begin{equation} \left(|\mathrm{e}^{\varOmega(z)}|^2\right)_s+w\operatorname{Im}\left[ \mathrm{e}^{2\mathrm{i}\phi}\, \mathrm{e}^{2\varOmega(z)}\right]=0\quad \text{on }\partial D. \end{equation}

For example, $\varOmega (z)=-\mathrm {i}\alpha$ corresponds to a director angle $\theta (x,y)=\alpha$ uniform in space; the principal-value logarithm $\varOmega (z) = -M\log (\mathrm {e}^{\mathrm {i}\alpha }z)$ corresponds to an isolated defect of topological charge $M$ located at $z=0$ and oriented with $\arg z= \alpha$; and $\varOmega (z) = [{\rm \pi} /(2\mathrm {i}\log R)-1]\log z$ corresponds to a nematic contained in the annulus $1<|z|< R$ with strong normal and tangential anchoring on the inner ($|z|=1$) and outer ($|z|=R$) boundaries, respectively. This last configuration is often called the ‘magical spiral’ (de Gennes & Prost Reference de Gennes and Prost1993) and has $\varUpsilon = {\rm \pi}^2/\log R$ and $M=1$ in (2.8).

The net free energy, (2.5), may now be written in terms of $\varOmega (z)$ as

(2.10)\begin{equation} \hat{\mathcal{E}} = \frac{1}{4}\oint_{\partial D}\operatorname{Im}\left[ \left(\varOmega(z)- \overline{\varOmega(z)}\right) \varOmega'(z)z_s\right]+ w \operatorname{Re}\left[1-\mathrm{e}^{\varOmega(z)-\overline{\varOmega(z)}}\, \mathrm{e}^{2\mathrm{i}\phi} \right]\, \mathrm{d} s, \end{equation}

where the bar denotes a complex conjugate. Additionally, the surface traction $\hat {\boldsymbol {t}}\equiv (\hat {t}_x,\hat {t}_y)$ given by (2.4) may be written as

(2.11)\begin{equation} \hat{t}_x-\mathrm{i} \hat{t}_y= \frac{1}{2\mathrm{i}}\varOmega'(z)^2\frac{\partial z}{\partial s} +\frac{w}{8}\left[(2+\mathrm{e}^{{-}2\mathrm{i}\phi}\, \mathrm{e}^{\overline{\varOmega(z)}-\varOmega(z)}-3\,\mathrm{e}^{2\mathrm{i}\phi}\, \mathrm{e}^{\varOmega(z)-\overline{\varOmega(z)}})\frac{\partial \bar z}{\partial s}\right]_s. \end{equation}

Overall, given a locally holomorphic function, $\varOmega (z)$, which satisfies the weak anchoring condition (2.9) with a period (2.8), the liquid crystal director angle, $\theta (x,y)=-\operatorname {Im}\varOmega (z)$; free energy, (2.10); and surface traction, (2.11), can all be determined.

2.3. Net body force and torque

The net dimensionless force, $(\hat {F}_x,\hat {F}_y)$, and torque, $\hat {T}$, acting on the body are found by integrating the surface traction, $(\hat {t}_x,\hat {t}_y)$ given by (2.11), and surface moment, $x \hat {t}_y-y\hat {t}_x$, around $\partial D$, respectively. Integrating by parts, and imposing the boundary condition (2.9) and the fact that $\exp [2\mathrm {i}(\theta -\phi )]\equiv \exp (\bar \varOmega -\varOmega -2\mathrm {i}\phi )$ is single-valued, results in the complex contour integrals

(2.12a)$$\begin{gather} \hat{F}_x-\mathrm{i} \hat{F}_y =\oint_{\partial D} \hat{t}_x-\mathrm{i} \hat{t}_y\, \mathrm{d} s = \frac{1}{2\mathrm{i}}\oint_{\partial D} \varOmega'(z)^2\, \mathrm{d} z, \end{gather}$$

and

(2.12b)$$\begin{gather}\hat{T} =\oint_{\partial D}x \hat{t}_y -y \hat{t}_x\, \mathrm{d} s =\frac{1}{2}\operatorname{Re}\left[\oint_{\partial D} z\varOmega'(z)^2\, \mathrm{d} z\right]+\varUpsilon, \end{gather}$$

where $\varUpsilon$ is the period defined in (2.8) and the centre of rotation is assumed to be $z=0$, without loss of generality.

Assuming there are no defects within the liquid crystal, i.e. $\varOmega (z)$ is holomorphic in $D$, we are free to deform the integration contours in (2.12) provided they remain outside the body. For a bounded director angle, $\varUpsilon = 0$ and $\varOmega '(z) \sim -M/z$ as $|z|\to \infty$ for some half-integer $M$ corresponding to the topological charge of the body, as introduced in (2.8). Thus, by taking the contours in (2.12) to infinity, we immediately find that there is no net force or torque acting on a immersed connected body when the director angle is bounded (i.e. $\hat {F}_x=\hat {F}_y=\hat {T}=0$). Note that a non-zero torque ($\hat {T}\neq 0$) induces a logarithmically growing director angle, that is $\theta (z)\equiv -\operatorname {Im}\varOmega (z)\sim M\arg z -\hat {T}\log |z|/[2{\rm \pi} (1-M)]$ as $|z|\to \infty$; i.e. $\varUpsilon =\hat {T}/(1-M)$ in (2.8).

A non-zero net force or torque on a body at equilibrium requires either a secondary body or the existence of a defect within the liquid crystal. For example, consider the complex director angle $\varOmega (z)= M\log [(z-\epsilon )/(z+\epsilon )]$ for a half-integer $M$. This corresponds to an unconfined director field with two defects at $z=\pm \epsilon$ of topological charge $\mp M$, respectively. The net force acting on one of these defects can be computed by integrating $\varOmega '(z)^2= 4M^2\epsilon ^2/(z^2-\epsilon ^2)^2$ around a contour only containing said defect, i.e. (2.12a). Cauchy's residue theorem then yields the forces $(F_x^{\pm },F_y^\pm ) = (\pm {\rm \pi}M^2/\epsilon,0)$ for the defects of charge $\pm M$, respectively. We hence find that the defects are attracted to each other with a force proportional to the inverse of the separation distance, and the square of their topological charge (de Gennes & Prost Reference de Gennes and Prost1993; Gartland et al. Reference Gartland, Sonnet and Virga2002; Harth & Stannarius Reference Harth and Stannarius2020), though their mobilities and migration speeds in a dynamic setting can differ substantially (Tóth, Denniston & Yeomans Reference Tóth, Denniston and Yeomans2002; Tang & Selinger Reference Tang and Selinger2019).

Analogy with d'Alembert's paradox. A classical and paradoxical result in potential fluid flow states that an inviscid fluid presents no resistance to a body which moves through it with steady translational motion (Batchelor Reference Batchelor1967). The resolution of this paradox is found with the incorporation of a viscous boundary layer on the body, which promotes the separation of vorticity from the surface in its wake. Here, we find an analogous result: there is no net force or torque acting on a immersed connected body when the director angle is bounded (i.e. $\hat {F}_x=\hat {F}_y=\hat {T}=0$). Unlike in potential flow theory, however, this result does not challenge our physical experience, and thus we do not seek further resolution in this setting.

Analogy with Stokes’ paradox. Another classical fluid mechanical paradox appears in Stokes flow. Namely that the flow due to a cylinder upon which there is a non-zero force is logarithmic in the distance from the body, and so the fluid velocity in an infinite domain is unbounded (Batchelor Reference Batchelor1967). The resolution of this apparent paradox is found by the reintroduction of inertial effects (e.g. via the Oseen equations). The analogous result here is that the director angle far from a body immersed in a LC with a non-zero torque is, similarly, logarithmic, and thus the associated elastic energy is unbounded. The question of solving for the equilibrium director field around a two-dimensional body with a non-zero torque is therefore similarly ill-posed. Reintroducing LC dynamics would regularize this singular behaviour by presenting a finite speed of propagation of information from near the body outward towards infinity.

3. Effective boundary technique for finite anchoring strength

While infinite anchoring ($w=\infty$) is mathematically appealing, real anchoring energies are finite. We now turn to the case of large, but finite, anchoring strengths. As $w\to \infty$, the strong anchoring boundary condition is recovered from (2.9), i.e.

(3.1)\begin{equation} \operatorname{Im} \,\mathrm{e}^{\varOmega(z)+\mathrm{i}\phi}=O(1/w)\quad \text{on }\partial D. \end{equation}

This boundary condition is invariant under a conformal map, thus the corresponding problem can be solved using a wide range of complex analysis tools (Ablowitz & Fokas Reference Ablowitz and Fokas2003). However, solutions to these Dirichlet problems are known to require singularities at boundaries, which correspond to ‘surface defects’ within the liquid crystal. This is problematic as defects give rise to an infinite free energy. This issue is often resolved by introducing a ‘melting region’ local to the defects (de Gennes & Prost Reference de Gennes and Prost1993), or by studying a more general alignment tensor theory in which this melting from a nematic to isotropic phase is tracked by a Maier–Saupe scalar order parameter (Hess Reference Hess1975; Sonnet, Kilian & Hess Reference Sonnet, Kilian and Hess1995; Stark & Lubensky Reference Stark and Lubensky2003). Here, we instead seek a correction to the strong anchoring constraint as $w\to \infty$.

For convenience, we denote $\varPhi (s):= \mathrm {e}^{\varOmega (z)+\mathrm {i}\phi }$, which is defined for $z(s)\in \partial D$. The boundary condition, (2.9), can then be written as the nonlinear constraint

(3.2)\begin{equation} \operatorname{Re} \varPhi(s)\operatorname{Re} \varPhi'(s)+\operatorname{Im} \varPhi(s)\operatorname{Im} \varPhi'(s)+w\operatorname{Im} \varPhi(s)\operatorname{Re} \varPhi(s)=0. \end{equation}

Using the fact that $\operatorname {Im} \varPhi = O(1/w)$ as $w\to \infty$, i.e. (3.1), we immediately find that

(3.3)\begin{equation} \operatorname{Im}\varPhi(s)+\frac{1}{w}\operatorname{Re}\varPhi'(s)=O(1/w^3), \end{equation}

as $w\to \infty$.

The linear constraint (3.3) contains the desired correction to the strong anchoring condition when the anchoring strength is large, but finite. This constraint is analogous to a Robin boundary condition for $\operatorname {Im} \varPhi$; however, unlike its Dirichlet counterpart (when $w=\infty$), it is not preserved under a conformal map and, thus, the most appealing tools of complex analysis do not immediately apply.

It is known, however, that the configuration of a nematic with weak anchoring (finite $w$) near a flat boundary is equivalent to that of a nematic with strong anchoring ($w=\infty$) on a different surface internal to the boundary, which is recessed from the surface by the extrapolation length $K/W\equiv a/w$ (Rapini & Papoular Reference Rapini and Papoular1969; de Gennes & Prost Reference de Gennes and Prost1993). Below, we use this idea to show that imposing the finite anchoring constraint (3.3) on $\partial D$ is mathematically equivalent to imposing the strong anchoring constraint (3.1) on some ‘effective domain’ boundary, $\partial D_w$, which is the physical boundary moved by the extrapolation length. The effective boundary, $\partial D_w$, is sketched in figure 1.

The above result follows from noting that (3.3) can be written as

(3.4)\begin{equation} \operatorname{Im}\varPhi(s+\mathrm{i}/w)=O(1/w^3), \end{equation}

where we define $\varPhi (s+\mathrm {i}/w):= \varPhi (s)+\varPhi '(s)\mathrm {i}/w- \varPhi ''(s)/(2w^2) +O(1/w^3)$. Substituting in $\varPhi =\mathrm {e}^{\varOmega (z)+\mathrm {i}\phi }$ then yields the boundary condition

(3.5)\begin{equation} \operatorname{Im} \, \mathrm{e}^{\varOmega(z)+\mathrm{i}\tilde{\phi}} = O(1/w^3)\quad \text{on }\partial D_{w}, \end{equation}

where $\tilde {\phi }(s):= \phi (s)-\phi ''(s)/(2w^2) + O(1/w^3)$ and $\partial D_w$ is the boundary $\partial D$ displaced by $-\boldsymbol {\hat {\nu }}(s)/w-\boldsymbol {\hat {s}}'(s)/(2w^2) +O(1/w^3)$ for the unit normal and tangent vectors $\boldsymbol {\hat {\nu }}(s)=-\boldsymbol {x}'(s)^\perp$ and $\boldsymbol {\hat {s}}(s)=\boldsymbol {x}'(s)$, respectively – i.e. $z(s)\in \partial D$ is mapped to $z(s+\mathrm {i}/w)\in \partial D_w$.

Above, we have implicitly assumed that the boundary curve, $\partial D$, is non-singular; however, domains with corners are ubiquitous in the literature (Lapointe et al. Reference Lapointe, Mason and Smalyukh2009; Davidson & Mottram Reference Davidson and Mottram2012; Muševič Reference Muševič2017). A local asymptotic analysis shows that mapping a corner in $\partial D$ to a similar corner in $\partial D_w$ is accurate up to $O(1/w^3)$ as $w\to \infty$ – i.e. the same order of accuracy as (3.5).

Overall, we have shown that finding a locally holomorphic function in $D$ that satisfies (2.9) on $\partial D$ is mathematically equivalent to finding a locally holomorphic function in $D_w$ that satisfies (3.5) on $\partial D_w$, up to $O(1/w^3)$ as $w\to \infty$. Here, we have analytically continued $\varOmega (z)$ inside the narrow region $D_w\setminus D$. It should be noted that the modified anchoring angle in (3.5), i.e. $\tilde \phi (s)\sim \phi (s)-\phi ''(s)/(2w^2)$, may differ from the original anchoring angle in (2.9), i.e. $\phi (s)$.

The equivalence of these two problems shows that any defects on the immersed boundary for $w=\infty$ are located inside the boundary for finite $w$, thus regularizing the energy. Furthermore, it provides a method for solving the finite anchoring problem on a domain $D$: first one constructs the effective domain, $D_w$, via displacing $\partial D$ by $-\boldsymbol {\hat {\nu }}(s)/w-\boldsymbol {\hat {s}}'(s)/(2w^2) +O(1/w^3)$, whilst ensuring any corner angles are preserved; then one solves the strongly anchored problem on $D_w$, for which the standard techniques of complex analysis, in particular conformal mapping, may be used. (As a first approximation, one can instead move the boundary a distance $1/w$ in the normal direction, however, the solution will then only be accurate to $O(1/w^2)$ as $w\to \infty$. In some cases, higher-order terms are required for computing the net energy, (2.10), and force, (2.12)).

The above ‘effective boundary technique’ can also be understood by considering the mapping from the physical domain, $D$, with boundary condition (2.9), to the effective domain, $D_w$, with boundary condition (3.5), as a conformal map. The transformation of the boundary condition can then be interpreted (at least for analytic boundaries) as an embedding in the Schwarz function of the boundary curve (Shapiro Reference Shapiro1992; Ablowitz & Fokas Reference Ablowitz and Fokas2003). Such encoding of information using Schwarz functions has been recently studied in the context of travelling water waves (Crowdy Reference Crowdy2023).

3.1. Tangential (homogeneous) anchoring and complex potentials

Although the results presented thus far hold for any anchoring angle, $\phi$, we are particularly interested in the case of tangential (homogeneous) anchoring with $\phi$ being tangent to $\partial D$. For tangential anchoring, $\phi (s) := \arg z'(s)\mod {\rm \pi}$ for $z(s)\in \partial D$. Furthermore, since

(3.6)\begin{equation} \arg\left[z'(s)+\frac{\mathrm{i}}{w}z''(s)-\frac{1}{2w^2}z'''(s)\right]= \phi(s)-\frac{1}{2w^2}\phi''(s)+O(1/w^3)\mod{\rm \pi}, \end{equation}

as $w\to \infty$, the modified anchoring angle, $\tilde \phi (s)\sim \phi (s)-\phi ''(s)/(2w^2)$, is tangent to the effective domain, i.e. $\tilde \phi (s) := \arg \tilde {z}'(s) \mod {\rm \pi}$ for $\tilde {z}(s)\in \partial D_w$. With this, the effective boundary technique (derived in § 3) says that imposing weak tangential anchoring on $\partial D$ is equivalent to imposing strong tangential anchoring on $\partial D_w$ up to $O(1/w^3)$ as $w\to \infty$. (This equivalency also holds for normal, homeotropic, anchoring with $\phi (s)={\rm \pi} /2+\arg z'(s) \mod {\rm \pi}$ for $z(s)\in \partial D$ or $\partial D_w$. Moreover, since a homeotropic director field is perpendicular to its homogeneous counterpart, one can transform a homeotropic anchoring problem into a homogeneous anchoring problem by taking $\varOmega _{homeo}(z) = \varOmega _{homo}(z)+\mathrm {i}{\rm \pi} /2$.)

Strong tangential anchoring is reminiscent of a vanishing boundary-flux problem in potential theory with $\varOmega (z)$ corresponding to the logarithm of the derivative of a complex potential $f(z)$, i.e. $\varOmega (z) = \log f'(z)$ (Acheson Reference Acheson1990; Ablowitz & Fokas Reference Ablowitz and Fokas2003); for example, $f(z)=z \mathrm {e}^{-\mathrm {i} \alpha }$ corresponds to a director angle $\theta =\alpha$ uniform in space. However, the problem here differs from the typical potential problem since $\varOmega (z)$ is not necessarily single-valued due to the periods in (2.8).

Inspired by the above discussion, it is natural to remove the multivaluedness of $\varOmega (z)$ then introduce an analogous complex potential. That is, we write

(3.7)\begin{equation} \varOmega(z) = \log f'(z) + g(z), \end{equation}

where $g(z)$ is any locally holomorphic function that accounts for the period in (2.8), i.e.

(3.8)\begin{equation} \oint_{\partial D_w} \mathrm{d} g = \varUpsilon - 2{\rm \pi} \mathrm{i} M, \end{equation}

and $f'(z)$ is a single-valued holomorphic function that accounts for the boundary condition on $\partial D_w$, i.e. (2.9). By writing the complex director field in this form, we shall recover a formulation for both $f(z)$ and $g(z)$ which is ubiquitous within the potential theory literature. This will allow for the use of known analytical solutions in the case of tangential anchoring (in §§ 4 and 5).

If we restrict our attention to the case when $M=0$ in (2.8) (i.e. the body, $\partial D$, has a vanishing topological charge), we then have enough freedom to choose $g(z)$ such that

(3.9)\begin{equation} \operatorname{Im} g=0\quad \text{on }\partial D_w. \end{equation}

The boundary condition (3.5) can then be integrated once to obtain

(3.10)\begin{equation} \operatorname{Im} f= O(1/w^3)\quad \text{on }\partial D_{w}, \end{equation}

where we have fixed the gauge of $f$ such that the integration constant vanishes. However, the boundary condition (3.10) does not necessarily yield a unique solution since, although $f'(z)$ is single-valued, its primitive may be multivalued. Uniqueness is restored by specifying the period of $f(z)$ around $\partial D$, i.e.

(3.11)\begin{equation} \oint_{\partial D} \mathrm{d} f =\varGamma, \end{equation}

for some real constant $\varGamma$. This period is analogous to specifying a constant circulation around a body in potential theory (Acheson Reference Acheson1990; Ablowitz & Fokas Reference Ablowitz and Fokas2003; Vitelli & Nelson Reference Vitelli and Nelson2006; Crowdy Reference Crowdy2020), whilst here it selects the locations of topological defects on the surface, as we shall see.

Overall, when $M= 0$, the homogeneous anchoring problem comes down to finding two locally holomorphic potentials, $g(z)$ and $f(z)$, which satisfy the boundary conditions (3.9) and (3.10) with given periods (2.8) and (3.11), respectively. Although this formulation may appear daunting, $g(z)$ and $f(z)$ are both equivalent to a complex potential in a typical potential theory problem (Acheson Reference Acheson1990; Ablowitz & Fokas Reference Ablowitz and Fokas2003; Crowdy Reference Crowdy2020). We can, thus, utilize the vast number of results which already exist within the literature.

4. Example 1: immersed cylinder with tangential anchoring

Consider a liquid crystal outside a cylinder with dimensionless unit radius, we denote this region $D$. The liquid crystal is assumed to be oriented with the $x$-axis in the far-field, thus $\varOmega (z)\to 0$ as $|z|\to \infty$, and subject to finite tangential (homogeneous) anchoring on the cylinder, i.e. $\varOmega (z)$ satisfies (3.3) with $\phi =\arg z'(s)\mod {\rm \pi}$ on $|z|=1$. (Although we have assumed finite anchoring with $w\gg 1$, the solution derived in this section ultimately satisfies the weak anchoring condition, (2.9), for any $w\geq 0$.) This configuration is plotted in figure 2.

Figure 2.

(Example 1.) Two-dimensional liquid crystal outside a unit cylinder ($\partial D$, black solid curve) subject to finite tangential anchoring. The effective domain boundary is a shrunken cylinder of radius $\rho (w)$ ($\partial D_w$, black dotted curve). The director field inside the effective domain, (4.3), is shown as faded blue curves for $w=10$ and period which minimizes the free energy, $\varGamma =0$.

Since $\varOmega (z)$ is holomorphic in $D$, the contour in the period-defining integral (2.8) can be deformed freely within the liquid crystal. Sending this contour off to infinity and imposing the far-field conditions, it follows that $\varOmega (z)$ must be single-valued with the periods in (2.8) vanishing (i.e. $\varUpsilon =M=0$). Following § 3.1, we introduce a locally holomorphic potential, $f(z)$, defined such that $\varOmega (z)=\log f'(z)$ and with its period around the cylinder as yet free, i.e. $\varGamma$ in (3.11). (Note that $g(z)$ is omitted here since $\varOmega (z)$ is already single-valued.)

The first step of the effective boundary technique is to construct the effective domain, $D_w$. Since the normal and tangent vectors of the cylinder are $\boldsymbol {\hat {\nu }}(s)=(\cos s,\sin s)$ and $\boldsymbol {\hat {s}}(s)=(-\sin s,\cos s)$, respectively, where $s$ is the complex argument, the unit cylinder is displaced by $-1/w+1/(2w^2)$ in the normal direction. Thus, the effective domain, $D_w$, is the region outside a cylinder of dimensionless radius $\rho (w):= 1-1/w+1/(2w^2)+O(1/w^3)$ – as shown in figure 2. Analytically continuing $f(z)$ inside the annulus $\rho (w)<|z|<1$ then yields the equivalent problem: find an $f(z)$ such that

(4.1a)\begin{align} f(z) \text{ locally holomorphic} &\quad \text{in }D_w, \end{align}

(4.1b)\begin{align} \operatorname{Im} f(z)=0 &\quad \text{on }\partial D_w, \end{align}

(4.1c)\begin{align} f(z)\sim z &\quad \text{as }|z|\to\infty \end{align}

and with period $\oint _{D_w}\, \mathrm {d} f = \varGamma$.

The potential problem (4.1) is a historic problem in the field of fluid dynamics, corresponding to the potential flow around a circular cylinder. Its solution can be found in many textbooks (Acheson Reference Acheson1990; Ablowitz & Fokas Reference Ablowitz and Fokas2003), i.e.

(4.2)\begin{equation} f(z) = z+ \frac{\rho^2}{z}+\frac{\varGamma}{2{\rm \pi}\mathrm{i}}\log\frac{z}{\rho}. \end{equation}

This potential corresponds to the director angle

(4.3)\begin{equation} \theta(x,y)\equiv{-}\arg f'(z) = 2\arg z-\arg(z-\rho \,\mathrm{e}^{\mathrm{i}\gamma})-\arg(z+\rho \,\mathrm{e}^{-\mathrm{i}\gamma}), \end{equation}

from which it is found that the defect locations are set by the period, $\varGamma$, through the simple relation $\sin \gamma := \varGamma /(4{\rm \pi} \rho )$. In fact, this solution not only satisfies the finite anchoring condition (3.3) for $w\gg 1$, but also satisfies the weak anchoring condition (3.2) for any $w\geq 0$, provided $\varGamma =0$ and we set $\rho (w):= ({\sqrt {1+4/w^2}-2/w})^{1/2}$. We assume this choice of effective radius henceforth.

The solution (4.3) has previously been derived by Burylov & Raikher (Reference Burylov and Raikher1994) using a superposition of separable solutions to Laplace's equation. Furthermore, it is equivalent to an unrestricted director field with three defects points: two $-1$ defects at $z= \rho \,\mathrm {e}^{\mathrm {i}\gamma }$ and $z=- \rho \,\mathrm {e}^{-\mathrm {i}\gamma }$ and a $+2$ defect at $z=0$. Note that these defects do not reside in the bulk liquid crystal domain for finite $w$; however, the $-1$ defects do correspond to surface defects as $\rho (w\to \infty )\to 1$.

The period $\varGamma$, or equivalently the positions of the $-1$ defects, is chosen to minimize the energy of the static liquid crystal, (2.10). This energy reduces to the complex contour integral

(4.4)\begin{equation} \hat{\mathcal{E}} = \frac{1}{4}\operatorname{Im}\left[2{\rm \pi}\mathrm{i} w+\oint_{|z|=1} w\frac{z\, f'(z)}{\bar{f}'(1/z)} -\log\bar{f}'(1/z) \frac{f''(z)}{f'(z)}\, \mathrm{d} z\right], \end{equation}

after imposing the Schwarz function of the cylinder, $\bar z =1/z$, (Shapiro Reference Shapiro1992; Ablowitz & Fokas Reference Ablowitz and Fokas2003). Inserting (4.2) into (4.4) results in a contour integral, which can be evaluated using Cauchy's residue theorem via accounting for the poles at the three defect points: $z=0$, $\rho \,\mathrm {e}^{\mathrm {i}\gamma }$ and $-\rho \,\mathrm {e}^{-\mathrm {i}\gamma }$; this yields

(4.5)\begin{equation} \hat{\mathcal{E}} =\frac{\rm \pi}{2}w(1-\rho(w)^2)-{\rm \pi}\log |1-\rho(w)^2|-{\rm \pi}\log|1+\rho(w)^2 \mathrm{e}^{2\mathrm{i}\gamma}|. \end{equation}

Further details on evaluating (4.4) can be found in Appendix B.1.

It is evident from (4.5) that $\gamma =0$ minimizes the net free energy. This result is not surprising due to the up–down symmetry of the problem; however, this is not necessarily the case for other geometries (see § 5, for example). We also note that $\hat {\mathcal {E}}={\rm \pi} +{\rm \pi} \log (w/4) - {\rm \pi}/2\log |1-\varGamma /(4{\rm \pi} )|+ O(1/w)$ as $w\to \infty$, i.e. the energy is logarithmically singular in the strong anchoring limit (a consequence of the surface defects at $z\sim \pm \mathrm {e}^{\pm \mathrm {i}\gamma }$).

Above, we assumed the immersed cylinder only locally disturbs the liquid crystal, so that the director field is oriented at a constant angle in the far-field. Another physical, but energetically unfavourable, configuration is that the cylinder appears as a point defect of topological charge $M\neq 0$ when viewed from afar, as introduced in (2.8). This would correspond to imposing $\varOmega (z)\sim -M\log z$ as $|z|\to \infty$. In this case, one can show that the resulting director angle is

(4.6)\begin{equation} \theta(z)= (2-M)\arg z-\arg(z^{2(1-M)}-\rho_M^{2(1-M)}), \end{equation}

with effective cylinder radius, $\rho _M(w)$, given by

(4.7)\begin{equation} [\rho_M(w)]^{2-2M} := \sqrt{1+\frac{4(1-M)^2}{w^2}}-\frac{2(1-M)}{w}, \end{equation}

for any $w\geq 0$. This solution corresponds to $3-2M$ defects: $2(1-M)$ defects of strength $-1$ on $|z|=\rho _M(w)$ and a defect of strength $2-M$ at $z=0$.

We finish this example by noting that, since the potential problem (4.1) is invariant under a conformal map, the director angle (4.2) describes the director field outside any connected body, provided the conformal map from the wanted extended domain, $\zeta \in D_w$, to the punctured domain, $|z|\geq \rho (w)$, can be obtained. Although such a conformal map must exist by the Riemann mapping theorem, determining such a mapping is notoriously difficult even for simple domains (Ablowitz & Fokas Reference Ablowitz and Fokas2003). Furthermore, although the effective domain, $\zeta \in D_w$, is mapped onto $|z|\geq \rho (w)$, the original domain, $\zeta \in D$, may not be mapped onto $|z|\geq 1$ since size and curvature are not necessarily preserved under a conformal map. To demonstrate this, we now consider a slightly more complicated domain: the region outside an equilateral triangle.

5. Example 2: immersed triangle with tangential anchoring

We next consider an immersed equilateral triangle with vertices at the roots of $z^3=\mathrm {e}^{3\mathrm {i}\chi }$. As before, the liquid crystal is assumed to be oriented with the $x$-axis in the far-field, i.e. $\varOmega (z)\to 0$ as $|z|\to \infty$, and subject to finite tangential (homogeneous) anchoring on the triangle, i.e. $\varOmega (z)$ satisfies (3.3) with $\phi =\arg z'(s)\mod {\rm \pi}$ on $\partial D$. Assuming there are no defects within the liquid crystal, the complex director angle, $\varOmega$, must be single-valued. We may, thus, introduce a complex potential, $f(z)$, defined such that $\varOmega (z)=\log f'(z)$ with period $\varGamma$, as in (3.11). This configuration is plotted in figure 3.

Figure 3.

(Example 2.) Two-dimensional liquid crystal outside a triangle with corners at the roots of $z^3=\mathrm {e}^{3\mathrm {i}\chi }$ ($\partial D$, black solid curve) subject to weak (finite) tangential anchoring. The effective domain boundary is a similar triangle with corners at the roots of $z^3=(1-2/w)^3\mathrm {e}^{3\mathrm {i}\chi }$ ($\partial D_w$, black dotted curve). The director field inside the effective domain, (5.4), is shown as faded blue curves for $w=10$, $\chi ={\rm \pi} /4$, and period which minimizes the free energy, $\varGamma \approx 3.14$.

For large anchoring strengths ($w\gg 1$), the effective boundary technique is of use (§ 3). The first step is to construct the effective domain, $D_w$. As the tangent vectors are constant, the triangle edges are displaced by $1/w+O(1/w^3)$ in the normal direction with the corner angles preserved. Consequently, we find that the effective domain, $D_w$, is the region outside a triangle with corners at the roots of $z^3=(1-2/w)^3\mathrm {e}^{3\mathrm {i}\chi }$, as shown in figure 3. Analytically continuing $f(z)$ into $D_w\setminus D$ yields the same potential problem as (4.1), but with $D_w$ now denoting the domain outside the effective triangle.

Since the potential problem, (4.1), is invariant under a conformal map, the solution (4.2) can be used here, provided a conformal map from the triangular domain, $z\in D_w$, to the punctured domain, $|\zeta |\geq 1$, can be found. Such a conformal map can be constructed using a Schwarz–Christoffel mapping, which maps the unit disk (or upper half-plane) onto a simple polygon (Driscoll & Trefethen Reference Driscoll and Trefethen2002; Ablowitz & Fokas Reference Ablowitz and Fokas2003).

By the Riemann-mapping theorem, there are three real degrees of freedom when constructing a conformal map (Ablowitz & Fokas Reference Ablowitz and Fokas2003). We shall fix the preimages of the effective triangle corners (i.e. the roots of $z^3=(1-2/w)^3\,\mathrm {e}^{3\mathrm {i}\chi }$) to be the roots of $\zeta ^3=\mathrm {e}^{3\mathrm {i}\chi }$. The corresponding Schwarz–Christoffel mapping is

(5.1)\begin{equation} z(\zeta) = A+B\int^\zeta (s^3-\mathrm{e}^{3\mathrm{i}\chi})^{2/3}\frac{\mathrm{d} s}{s^2}, \end{equation}

for complex constants $A$ and $B$, which are fixed by ensuring that the corners of the triangle are mapped to their respective preimages. After evaluating the constants, the mapping (5.1) can be written as

(5.2)\begin{equation} z(\zeta) = \left(1-\frac{2}{w}\right)\frac{h(\mathrm{e}^{3\mathrm{i}\chi}/\zeta^{3})}{h(1)}\zeta, \end{equation}

for the hypergeometric function $h(\zeta ):= {}_2 F_1 (-2/3,-1/3;2/3;\zeta )$ (as defined by Abramowitz & Stegun (Reference Abramowitz and Stegun1964), for example).

In the punctured domain, $|\zeta |\geq 1$, the locally holomorphic function $F(\zeta ):= f(z(\zeta ))$ satisfies $\operatorname {Im} F(\zeta )=0$ on $|\zeta |=1$, $F(\zeta )\sim (1-2/w)\zeta /h(1)$ as $|\zeta |\to \infty$, and has a period $\varGamma$ around $|\zeta |=1$. This problem is equivalent to (4.1) and, thus, has the same solution form, i.e.

(5.3)\begin{equation} F(\zeta) =\frac{1-2/w}{h(1)}\left( \zeta+\frac{1}{\zeta}\right) +\frac{\varGamma}{2{\rm \pi} \mathrm{i}}\log \zeta. \end{equation}

Taking the derivative of (5.3) with respect to $z$ yields the director angle

(5.4)\begin{equation} \theta(z)\equiv{-}\arg f'(z) =\tfrac{2}{3}\arg[\zeta(z)^3-\mathrm{e}^{3\mathrm{i}\chi}]- \arg[\zeta(z)-\mathrm{e}^{\mathrm{i}\gamma}]- \arg[\zeta(z)+\mathrm{e}^{-\mathrm{i}\gamma}], \end{equation}

for $\sin \gamma := h(1)\varGamma /[4{\rm \pi} (1-2/w)]$ and $\zeta (z)$ given by the inverse of (5.2).

As in Example 1 (§ 4), we will determine the period, $\varGamma$, so as to minimize the net free energy of the liquid crystal, (2.10). However, unlike the cylindrical case, the free energy cannot be written as a closed integral of an analytic function since a triangle does not have an analytic Schwarz function, $\bar z = S(z)$, and the boundary angle, $\arg z_s$, is only defined piecewise (Shapiro Reference Shapiro1992; Ablowitz & Fokas Reference Ablowitz and Fokas2003). Instead we separate (2.10) into three finite contour integrals, corresponding to the three edges of the triangle. These contour integrals are then evaluated numerically for given $\varGamma$, $w$ and $\chi$. The resulting energy is then minimized numerically to determine $\varGamma =\varGamma _{min}(w,\chi )$ – we use fminsearch with numerical quadrature in MATLAB (further details on evaluating the energy can be found in Appendix B.2).

It is evident from (5.4) that there are five defect points on the effective triangle: three defects at the corners, i.e. at the roots of $z^3=(1-2/w)^3\,\mathrm {e}^{3\mathrm {i}\chi }$; and two $-1$ defects at $z=z(\pm \mathrm {e}^{\pm \mathrm {i}\gamma })$. Although the locations of the corner defects are fixed for any given $w$ and $\chi$, the $-1$ defects are dependent on $\varGamma =\varGamma _{min}(w,\chi )$, thus their locations are not immediately apparent. As $w\to \infty$, we observe that at least one of the $-1$ defects lies at a corner of the effective triangle (this can be seen approximately for $w=100$ in figure 4a, for example). This can be understood by noting that most of the liquid crystal deformation occurs in the vicinity of defects, thus the energy is minimized by minimizing the number of topological defects. For smaller anchoring strengths, however, the defects do not lie in the physical domain; thus, placing one of the $-1$ defects at a corner of the effective triangle is not necessarily the energy-minimizing solution (this can be seen for $w=10$ in figure 4${\bigcirc{\kern-6pt 3}}$, for example).

Figure 4.

(Example 2.) (a) Plot of the complex arguments of the two $-1$ defects as a function of the triangle orientation, $\chi$, for $w=10$ (dashed curve) and $w=100$ (solid curve). The rightmost defect is described by the blue curve with $\alpha _1 := \arg [z(\mathrm {e}^{\mathrm {i}\gamma _{min}})]$, whilst the leftmost defect is described by the red curve with $\alpha _2 := \arg [-z(-\mathrm {e}^{-\mathrm {i}\gamma _{min}})]$. Here, $\varGamma =\varGamma _{min} \equiv 4{\rm \pi} (1-2/w)\sin \gamma _{min} /h(1)$ is the period which minimizes the net free energy, (2.10). In ${\bigcirc{\kern-6pt 1}}$–${\bigcirc{\kern-6pt 3}}$, the positions of the $-1$ defects are plotted as coloured dots for $w=10$ and the labelled $\chi$. Note that at least one of the defects lies approximately at a corner of the effective triangle for $w=100$ (i.e. $\alpha _1-\chi \approx 0$ or $\alpha _2-\chi \approx \pm {\rm \pi}/3$), however, this is not the case for $w=10$ (as seen in ${\bigcirc{\kern-6pt 3}}$, for example). (b) Plot of the net free energy as a function of $\chi$, corresponding to the solutions presented in (a). For both $w=10$ (dashed curve) and $w=100$ (solid curve), it is evident that the free energy is smallest when the triangle is pointing upwards ($\chi =-{\rm \pi} /6$) or downwards ($\chi ={\rm \pi} /6$).

Analogy with the Kutta condition. Here we observe another analogy with the classical theory of potential flow past a body, in particular the placement of the rear stagnation point at a sharp trailing edge. The selection of the circulation which aligns these two is known as the Kutta condition, and the lift which results in that setting is known as the Kutta–Joukowski lift theorem (Acheson Reference Acheson1990). We observe a very similar structure in the case of strong anchoring ($w\to \infty$); however, there is no direction associated with a nematic, thus either defect could reside at a corner (there is no ‘trailing’ edge). Moreover, since the force and torque on the body at equilibrium must be zero (due to the analogy with Stokes’ paradox, as described in § 2.3), here there is no ‘lift’, regardless of the circulation/period.

For each orientation of the triangle, $\chi$, we are able to minimize the net free energy over the period $\varGamma$; thus, each $\chi$ has an energetic cost associated with it. In figure 4(b), this minimized free energy is plotted as a function of $\chi$ for $w=10$ and $w=100$. It is apparent that that this energy is smallest when the triangle is either pointed upwards ($\chi =-{\rm \pi} /6$) or downwards ($\chi ={\rm \pi} /6$): one of the sides of the triangle is aligned with the relaxed orientation of the liquid crystal, i.e. the $x$-axis. This observation appears to hold for all values of the anchoring strength, $w$. Furthermore, for these values of $\chi$, the two $-1$ defects approximately lie at the two vertices associated with the horizontal side of the effective triangle (see figure 4${\bigcirc{\kern-6pt 1}}$, for example). This is consistent with the experimental observations of Lapointe et al. (Reference Lapointe, Mason and Smalyukh2009), who observed that regular polygonal colloids reorient themselves until one of their sides aligns with the relaxed direction of the liquid crystal. The analysis of this phenomenon requires the introduction of dynamics within the liquid crystal and, thus, is left for future work.

6. Example 3: Mobility of a Janus particle with hybrid anchoring

For a final example, let us reconsider the oriented liquid crystal outside a cylinder with dimensionless unit radius (as was considered in § 4). Unlike the previous two examples, here we assume the liquid crystal is subject to a weak hybrid anchoring condition (finite $w$), which is partially tangential (homogeneous) and partially perpendicular (homeotropic). That is, we assume $\varOmega (z)$ satisfies (3.3) on $|z|=1$ with the anchoring angle

(6.1)\begin{equation} \phi(s) = \arg z'(s) +\begin{cases} 0 & \text{if $|\arg[z(s)\, \mathrm{e}^{-\mathrm{i}\beta}]|<\alpha/2$,}\\ {\rm \pi}/2 & \text{if $|\arg[z(s) \,\mathrm{e}^{-\mathrm{i}\beta}]|>\alpha/2$,} \end{cases}\mod{\rm \pi}, \end{equation}

where $-{\rm \pi} <\arg z\leq {\rm \pi}$ is the principal-value argument and $\alpha \in [0,2{\rm \pi} ]$ and $\beta \in (-{\rm \pi}, {\rm \pi}]$ are given angles. This boundary condition corresponds to imposing tangential anchoring on an arc of angle $\alpha$ and normal anchoring on the remaining circle. The angle $\beta$ corresponds to the orientation of the tangential section of the cylinder. Similar theoretical configurations have focused on Janus particles in three-dimensions (Conradi et al. Reference Conradi, Ravnik, Bele, Zorko, Žumer and Muševič2009; Sahu et al. Reference Sahu, Kole, Ramaswamy and Dhara2020) and in a two-dimensional film (Loewe & Shendruk Reference Loewe and Shendruk2022); related work on squirming particles may be found in Lintuvuori, Würger & Stratford (Reference Lintuvuori, Würger and Stratford2017), Daddi-Moussa-Ider & Menzel (Reference Daddi-Moussa-Ider and Menzel2018), Chi et al. (Reference Chi, Gavrikov, Berlyand and Aranson2022) and Berlyand et al. (Reference Berlyand, Chi, Potomkin and Yip2022). This configuration is plotted in figure 5.

Figure 5.

(Example 3.) Two-dimensional liquid crystal outside a unit cylinder ($\partial D$, solid black curve) subject to weak (finite) tangential anchoring on an arc of angle $\alpha$ and weak normal anchoring on the remaining cylinder. The effective domain boundary is a shrunken cylinder of radius $\rho (w)$ ($\partial D_w$, black dotted curve). The director field inside the effective domain, (6.6), is shown as faded blue curves for $w=10$, $\alpha ={\rm \pi} /4$, and an orientation which minimizes the free energy, $\beta ={\rm \pi} /2$.

Assuming there are no defects within the liquid crystal and it is oriented with the $x$-axis at infinity (i.e. $\varOmega (z)\to 0$ as $|z|\to \infty$), $\varOmega (z)$ must be single-valued with the period in (2.8) vanishing (i.e. $\varUpsilon =M=0$). Furthermore, since both tangential and normal anchoring conditions are conserved under the effective boundary technique (see § 3.1), imposing finite hybrid anchoring, i.e. (3.3) with (6.1), on $|z|=1$ is equivalent to imposing strong hybrid anchoring, i.e. (3.5) with (6.1), on the effective cylinder $|z|=\rho (w)$, where $\rho (w) := (\sqrt {1+4/w^2}-2/w)^{1/2}$ is the effective radius derived in § 4. Analytically continuing $\varOmega (z)$ inside the annulus $\rho (w)<|z|<1$ yields the problem: find a $\varOmega (z)$ such that

(6.2a)\begin{align} \varOmega(z) \text{ holomorphic} &\quad \text{in $|z|>\rho(w)$,} \end{align}

(6.2b)\begin{align} \operatorname{Im} \left[z_s\exp{\varOmega(z)}\right]=0 &\quad \text{on $|z|=\rho(w)$ with $\left|\arg[z \mathrm{e}^{-\mathrm{i}\beta}]\right|<\alpha/2$,} \end{align}

(6.2c)\begin{align} \operatorname{Re} \left[z_s\exp{\varOmega(z)}\right]=0 &\quad \text{on $|z|=\rho(w)$ with $\left|\arg[z \mathrm{e}^{-\mathrm{i}\beta}]\right|>\alpha/2$,} \end{align}

(6.2d)\begin{align} \varOmega(z)\sim 1 &\quad \text{as $|z|\to\infty$}. \end{align}

Since $z_s=\mathrm {i} z$ for a circular boundary, it is convenient to introduce the complex function $F(z) := \mathrm {i} z \exp \varOmega (z)$. It follows from (6.2) that $F(z)$ must be holomorphic inside the extended domain, has a pole at infinity, and has either the real or imaginary part vanishing on $|z|=\rho (w)$. This problem is conformally invariant, thus we are free to use a conformal map to derive a solution.

Consider the conformal map from the $z$-plane to the $\zeta$-plane:

(6.3)\begin{equation} \zeta(z)=\mathrm{e}^{\mathrm{i}\alpha/4} \sqrt{\frac{z-\rho(w) \,\mathrm{e}^{\mathrm{i}(\beta-\alpha/2)}}{z-\rho(w) \,\mathrm{e}^{\mathrm{i}(\beta+\alpha/2)}}}. \end{equation}

This map is a composition of two conformal maps: a Möbius transformation, which maps the punctured $z$-plane to the upper-half–plane; and the square-root map, which maps the upper-half–plane to the first quadrant, i.e. the $\zeta$-plane.

In the $\zeta$-plane, the complex function $G(\zeta ):= F(z(\zeta ))$ satisfies

(6.4a)\begin{align} G(\zeta) \text{ holomorphic} &\quad \text{in $\operatorname{Re}\zeta>0$ & $\operatorname{Im}\zeta>0$,} \end{align}

(6.4b)\begin{align} \operatorname{Im} G(\zeta)=0 &\quad \text{on $\operatorname{Im}\zeta=0$,} \end{align}

(6.4c)\begin{align} \operatorname{Re} G(\zeta)=0 &\quad \text{on $\operatorname{Re}\zeta=0$,} \end{align}

(6.4d)\begin{align} G(\zeta)\sim \frac{\rho(w)\sin(\alpha/2)\, \mathrm{e}^{\mathrm{i}\beta}}{1-\mathrm{e}^{-\mathrm{i}\alpha/4} \zeta} &\quad \text{as $\zeta\to \mathrm{e}^{\mathrm{i}\alpha/4}$}. \end{align}

Here, the pole at $z=\infty$ has been mapped to an interior point at $\zeta =\mathrm {e}^{\mathrm {i}\alpha /4}$. The boundary conditions (6.4b) and (6.4c) can be satisfied by introducing three additional poles at the images of $\zeta =\mathrm {e}^{\mathrm {i}\alpha /4}$ across the real and imaginary axes. This is commonly referred to as the method of images (Ablowitz & Fokas Reference Ablowitz and Fokas2003) and yields the solution to (6.4),

(6.5)\begin{equation} G(\zeta) = \frac{\rho\sin(\alpha/2)\,\mathrm{e}^{\mathrm{i}\beta}}{1-\mathrm{e}^{-\mathrm{i}\alpha/4} \zeta}+\frac{\rho\sin(\alpha/2)\,\mathrm{e}^{-\mathrm{i}\beta}}{1+\mathrm{e}^{\mathrm{i}\alpha/4} \zeta}-\frac{\rho\sin(\alpha/2)\,\mathrm{e}^{-\mathrm{i}\beta}}{1-\mathrm{e}^{\mathrm{i}\alpha/4} \zeta}-\frac{\rho\sin(\alpha/2)\,\mathrm{e}^{\mathrm{i}\beta}}{1+\mathrm{e}^{-\mathrm{i}\alpha/4} \zeta}. \end{equation}

Changing variables with (6.3) results in an expression for the director angle outside the cylinder,

(6.6) \begin{align} \theta(z) &\equiv{-}\arg[-\mathrm{i} F(z)/z] =2\arg z-\arg(z-\rho \,\mathrm{e}^{-\mathrm{i}\beta})\nonumber\\ &\quad-\tfrac{1}{2}\arg(z-\rho \,\mathrm{e}^{\mathrm{i}(\beta+\alpha/2)})- \tfrac{1}{2}\arg(z-\rho \,\mathrm{e}^{\mathrm{i}(\beta-\alpha/2)}). \end{align}

This solution is equivalent to an unrestricted director field with four defects: two $-1/2$ defects corresponding to the tangential–normal anchoring interfaces at $z=\rho \,\mathrm {e}^{\mathrm {i}(\beta \pm \alpha /2)}$; one $-1$ defect at $z=\rho \,\mathrm {e}^{-\mathrm {i}\beta }$; and one $+2$ defect at $z=0$.

When $\alpha =2{\rm \pi}$, the liquid crystal is subject to finite tangential anchoring on the entire cylinder (as considered in § 4). The director angle (6.6) recovers (4.3) in this case, with the body orientation, $\beta$, corresponding to the period $\gamma := \arcsin [\varGamma /(4{\rm \pi} \rho )]$. Alternatively, when $\alpha =0$, the liquid crystal is subject to finite normal anchoring on the entire cylinder. The resulting director angle is equivalent to taking $z\mapsto \mathrm {i} z$ in (4.3) – i.e. the equipotential lines of the complex potential (4.2) rotated by ${\rm \pi} /2$. The body orientation, $\beta$, is again analogous to the period $\gamma$ in this case. Finally, when $\alpha ={\rm \pi}$, the liquid crystal is subject to finite tangential/normal anchoring on exactly one half of the cylinder. This boundary condition has previously been termed ‘Janus anchoring’ due to the analogy with Janus particles (Loewe & Shendruk Reference Loewe and Shendruk2022). For strong Janus anchoring ($w=\infty$ and $\alpha ={\rm \pi}$), the director angle (6.6) is equivalent to a subcase of the solution derived by Loewe & Shendruk (Reference Loewe and Shendruk2022), who used a superposition of separable solutions to construct a series representation.

For any given value of $\alpha$, there exists a free energy associated with each orientation angle $\beta$, i.e. (2.10). It is, thus, natural to ask what orientation minimizes this net free energy. (This is equivalent to asking what period, $\gamma$, minimizes the net free energies in the previous two examples, §§ 4 and 5.) Since the anchoring angle (6.1) is only defined piecewise, the energy (2.10) cannot be written as a closed integral of an analytic function and, thus, cannot be evaluated using Cauchy's residue theorem (as was done in § 4, when $\alpha =2{\rm \pi}$, for example). We instead evaluate the energy using numerical quadrature for given values of the anchoring strength, $w$; sector angle, $\alpha$; and orientation angle, $\beta$. We are then able to minimize over the angle $\beta$ numerically. Further details on the evaluation of the energy integral, (2.10), along with some partial analytical progress can be found in Appendix B.3.

If the anchoring is predominately normal to the cylinder (i.e. $0\leq \alpha <{\rm \pi}$), we observe that the net free energy is minimized for a vertically aligned cylinder (i.e. $\beta =\pm {\rm \pi}/2$). Alternatively, if the anchoring is predominately tangential to the cylinder (i.e. ${\rm \pi} <\alpha \leq 2{\rm \pi}$), then the net free energy is minimized by a horizontally aligned cylinder (i.e. $\beta =0$ or $2{\rm \pi}$). In the critical case (i.e. Janus anchoring, $\alpha ={\rm \pi}$), both horizontal and vertical alignment minimize the free energy. This observation is supported by figure 6, which shows the net free energy as a function of the orientation angle, $\beta$, for various values of $\alpha$ and $w=10$. It also complements the analytical expression of the energy (4.5), which holds when $\alpha =2{\rm \pi}$. This orientation preference suggests that, if the cylindrical body was free to rotate, it would align itself with either the vertical or horizontal axes depending on the value of $\alpha$. Controlling the orientation of a body using surface anchoring has previously been studied in the context of microswimmers by Chi et al. (Reference Chi, Potomkin, Zhang, Berlyand and Aranson2020), and a related method of fluid transport by passing a travelling wave of preferred director angle has elsewhere been proposed by Krieger et al. (Reference Krieger, Dias and Powers2015). However, analysing such orientation control further requires the addition of dynamics within the liquid crystal and, thus, is left for future work.

Figure 6.

(Example 3.) Plot of the net free energy, (2.10), as a function of the orientation angle, $\beta$, for anchoring strength $w=10$ and various sector angles, $\alpha$, delineated by colour in the figure legend. It is evident that $\beta =-{\rm \pi} /2$ and $\beta ={\rm \pi} /2$ minimize the net free energy when $0\leq \alpha < {\rm \pi}$ (as shown in ${\bigcirc{\kern-6pt 1}}$, for example), whilst $\beta =0$ and $\beta =2{\rm \pi}$ minimize the net free energy when ${\rm \pi} <\alpha \leq 2{\rm \pi}$ (as shown in ${\bigcirc{\kern-6pt 2}}$, for example). In the critical case ($\alpha ={\rm \pi}$), all four orientation angles minimize the net free energy (as shown in ${\bigcirc{\kern-6pt 3}}$, for example).