# Encoding Gaussian curvature in glassy and elastomeric liquid crystal solids

## Abstract

We describe shape transitions of thin, solid nematic sheets with smooth, preprogrammed, in-plane director fields patterned across the surface causing spatially inhomogeneous local deformations. A metric description of the local deformations is used to study the intrinsic geometry of the resulting surfaces upon exposure to stimuli such as light and heat. We highlight specific patterns that encode constant Gaussian curvature of prescribed sign and magnitude. We present the first experimental results for such programmed solids, and they qualitatively support theory for both positive and negative Gaussian curvature morphing from flat sheets on stimulation by light or heat. We review logarithmic spiral patterns that generate cone/anti-cone surfaces, and introduce spiral director fields that encode non-localized positive and negative Gaussian curvature on punctured discs, including spherical caps and spherical spindles. Conditions are derived where these cap-like, photomechanically responsive regions can be anchored in inert substrates by designing solutions that ensure compatibility with the geometric constraints imposed by the surrounding media. This integration of such materials is a precondition for their exploitation in new devices. Finally, we consider the radial extension of such director fields to larger sheets using nematic textures defined on annular domains.

### 1. Introduction

It is well known that inhomogeneous local deformations in thin solid sheets, for instance differential growth in leaves, can lead to the formation of Gaussian curvature and complex shape transitions [1]. We are concerned with modern responsive materials that can be preprogrammed to undergo prescribed spatially inhomogeneous expansions and contractions in response to external stimuli, offering exciting possibilities for the design and production of switchable surfaces for a variety of applications [2–6]. In particular, liquid crystal (LC) solids, either glassy or elastomeric, have orientational order about a director ** n**, but unlike their classical fluid analogues, they cannot flow. Nematic liquid crystalline glasses and elastomers are particularly promising candidates as responsive materials. They are richer than isotropic systems, for instance those with spatially varying capacity to swell, in that both the direction

**of anisotropy as well as the degree of mechanical response are variable. Stimuli inducing such changes are typically light or heat, both of which can reduce the degree of order and cause contraction. The contraction λ<1 along**

*n***, and elongation by λ**

*n*^{−ν}>1 in the directions perpendicular to

**, are large in elastomers, where the optothermal Poisson ratio [7] $\nu =\frac{1}{2}$ relates the perpendicular and parallel responses (with constancy of volume). Glasses have λ∼ 0.90–1.0 only, and their $\nu \in (\frac{1}{2},2)$ (volume increase).**

*n*Liquid crystal elastomers have a subtle mechanics, because director rotation can accompany some imposed deformation and there can be little or no stress—a remarkable response termed soft elasticity. Conceivably, soft elasticity could mask the Gaussian curvature we explore, except that (i) heating all the way to the isotropic phase eliminates order and hence anisotropy that can be redirected in order to relieve stress, and (ii) we find states that are unstretched with respect to the new metric, and there is little incentive for director rotation. We henceforth ignore soft elasticity and treat nematic glasses and elastomers in the same way, until we come to discuss actuation and device implications.

The spontaneous deformation gradient tensor in nematic glasses [8] and elastomers [9] is of the form

_{3}denotes the identity operator on ${\mathbb{R}}^{3}$. When

*F*is achieved, the heated or irradiated body is in a relaxed state without stresses.

The distortions associated with heating a sheet with uniform ** n** are quite spectacular, but more subtle and perhaps more useful effects are achieved with a programmed-in, non-uniform

**. Gaussian curvature arises, because the metric, specifying the intrinsic geometry of the deformed plane, spatially varies. The first example investigated [4] was that generated by concentric circles of**

*n***. One immediately sees, if circumferences contract by a factor of λ<1 and radii extend by λ**

*n*^{−ν}>1, the circles now sit on what has morphed into a cone, concentric with its tip where the localized Gaussian curvature resides. Experimentally, effects are large, even in glasses [10], and still more so in elastomers [11]. These director circles are topological charge

*q*=1 defects (disclinations) in a two-dimensional director field. Other charges instead give delocalized Gaussian curvature [7] of both signs. Programming

**into LC sheets is possible at such precision that**

*n**q*=6 defects can be created in elastomers, and in the case of glasses, even

*q*=±10 defects have been produced which generate highly complex surfaces [12,11]. In optics, such structures are used as axially symmetric wave plates [13], for instance with charge

*q*=64, or termed ‘

*q*-plates’, producing sophisticated light polarization, with charge

*q*=100 being employed [14].

More general distributions of Gaussian curvature by variation of director in plane were attacked in the seminal work of Aharoni *et al.* [5] that describes the interplay between the nematic director field of a thin elastomeric sheet and the resulting three-dimensional configuration attained upon heating. In particular, they consider the difficult reverse problem of constructing a director field that induces a specified two-dimensional metric. In cases of one-dimensional variation (leading to surfaces of revolution), they have a direct way of finding the metric. These authors also go on to consider more ambitious problems such as two-dimensional variation. In this paper, we follow [6] and consider some specific in-plane director field patterns on thin nematic sheets that give interesting curvature distributions and consequent shapes that are then compared qualitatively with our experiments. We continue by considering spiral director patterns that produce more advanced shapes—cones, hyperbolic cones, (pseudo)spherical caps and spindles. An additional but important constraint is then added in—that the generated shape is geometrically compatible with a surrounding, inert sheet, so that the morphing structure can be anchored. Practical devices will require such attachment, as they respond to stimuli, in order to exploit induced Gaussian curvature to pump, do work or actuate. We will conclude by discussing how applications might be achieved because of this essential advance. We now review the mathematical framework [5,6] used in the subsequent analysis.

### 2. Generating intrinsically curved surfaces

#### (a) Gaussian curvature through patterning

It is assumed that the director field does not vary across the thickness of the sheet, so that the same pattern is repeated at each level of thickness. For sufficiently thin sheets, stimulation of the system will result in pure bending of the sheet at no stretch energy cost and one expects an isometric immersion of the prescribed local deformations as determined by the director field pattern.

Let $({x}_{1},{x}_{2})\in \omega \subset {\mathbb{R}}^{2}$ be Cartesian coordinates parametrizing the midsurface of the initially flat sheet and $\mathit{n}({x}_{1},{x}_{2})={n}_{1}{\hat{\mathit{e}}}^{1}+{n}_{2}{\hat{\mathit{e}}}^{2}$ be the director field pattern across the surface, where ${\hat{\mathit{e}}}^{1}$, ${\hat{\mathit{e}}}^{2}$ form the standard orthonormal basis of ${\mathbb{R}}^{2}$. The associated in-plane spontaneous deformation tensor *F* has components *F*_{αβ}=(λ−λ^{−ν})*n*_{α}*n*_{β}+λ^{−ν}*δ*_{αβ}, where *α*, *β*=1,2. The resulting two-dimensional metric of the deformed sheet upon stimulation is *a*=*F*^{T}*F*, which simplifies to

*K*of a surface is an intrinsic geometric property (

*Theorema Egregium*) that is determined by the first fundamental form

*a*

_{αβ}of the surface via

*a*

^{μρ})=(

*a*

_{μρ})

^{−1}.

Ignoring spatial variation in λ that would also contribute to the Christoffel symbols ${\mathrm{\Gamma}}_{\mu \rho}^{\tau}$ and hence *K*,^{1} we concentrate on the two-dimensional director field that is characterized by an angle scalar field *ψ*=*ψ*(*x*_{1},*x*_{2}) specifying the in-plane orientation of the director at each point on the initially flat sheet, so that ${n}_{1}=\mathrm{cos}\psi $ and ${n}_{2}=\mathrm{sin}\psi $. The Gaussian curvature determined by the nematic metric can be expressed in terms of *ψ* as

*π*/2 in-plane at every point, so that $\psi \to \psi +\pi /2$, then the resulting Gaussian curvature flips sign everywhere, because $\mathrm{sin}2\psi \to -\mathrm{sin}2\psi $ and $\mathrm{cos}2\psi \to -\mathrm{cos}2\psi $. That is,

*π*/2 radian rotation of the directors as orthogonal duals, pairs that relate to heating and cooling.

^{2}In this paper we consider only λ≤1, that is heating; we achieve changes in sign of

*K*by changing the

*n*(

*r*).

#### (b) Spherical and pseudo-spherical surfaces

We now restrict our attention to director fields of the form $\mathit{n}=\mathrm{cos}\psi ({x}_{2})\hspace{0.17em}{\hat{\mathit{e}}}^{1}+\mathrm{sin}\psi ({x}_{2})\hspace{0.17em}{\hat{\mathit{e}}}^{2}$, whose alignment angle field varies only with respect to one of the coordinates. The Gaussian curvature upon stimulation is $K=-\frac{1}{2}({\mathrm{\lambda}}^{-2}-{\mathrm{\lambda}}^{2\nu})({\psi}^{\u2033}\mathrm{sin}2\psi +2{\psi}^{\prime 2}\mathrm{cos}2\psi ).$ We can rewrite this as

*C*(

*K*)=

*K*/(λ

^{−2}−λ

^{2ν}), and solve for constant

*K*>0 to find

*c*

_{1},

*c*

_{2}are constants of integration. This pattern generates constant Gaussian curvature

*K*wherever it is well defined. Now, consider the particular solution

*c*

_{1}=1,

*c*

_{2}=−4/

*L*, and

*x*

_{2}≤2

*L*. The pattern on the square domain

*ω*=[0,2

*L*]×[0,2

*L*] is shown in figure 1

*a*. By integrating along the director field lines, we note that the integral curves of this pattern consist of semicircles of radius

*L*that are shifted along the

*x*

_{1}-axis. Thus, if we seek to encode a particular constant positive Gaussian curvature

*K*=

*K*

_{0}>0 across an initially flat sheet, then we can do so by encoding the pattern obtained by shifting a semicircle of radius $L=(1/\sqrt{{K}_{0}})\sqrt{{\mathrm{\lambda}}^{-2}-{\mathrm{\lambda}}^{2\nu}}$ as shown in figure 1

*b*.

By the observation that the orthogonal dual of a given two-dimensional director field pattern generates the exact opposite Gaussian curvature at every point, we can encode constant negative Gaussian curvature *K*=−*K*_{0} on a thin nematic sheet by simply using the orthogonal dual of the pattern that encodes positive curvature *K*_{0}>0. Returning to the example of figure 1*a*, where a pattern encoding constant positive curvature *K*=(λ^{−2}−λ^{2ν})/*L*^{2} was defined on the square domain *ω*=[0,2*L*]×[0,2*L*] by $\psi ({x}_{2})=\mathrm{arccos}(1-{x}_{2}/L)$, we immediately obtain a pattern on the same domain which encodes constant negative Gaussian curvature *K*=−(λ^{−2}−λ^{2ν})/*L*^{2}, by simply taking $\psi ({x}_{2})=\mathrm{arccos}(1-{x}_{2}/L)+\pi /2$. The resulting pattern is shown in figure 1*c*. This pattern is generated by shifting a tractrix curve along its axis.

For a surface in ${\mathbb{R}}^{3}$, the components *a*_{αβ} and *b*_{αβ} of the first and second fundamental forms satisfy a system of algebraic differential equations known as the Gauss–Codazzi–Mainardi equations. Conversely, any pair (*a*,*b*) consisting of a symmetric and positive definite matrix field (*a*_{αβ}) and a symmetric matrix field (*b*_{αβ}) that satisfy the Gauss–Codazzi–Mainardi equations determines a unique surface up to a rigid transformation in ${\mathbb{R}}^{3}$ [15]. Thus, to determine the equilibrium configuration of the mid-surface of an initially flat nematic sheet upon stimulation, we also need to know the components *b*_{αβ} of the second fundamental form that minimize the bending energy subject to the Gauss–Codazzi–Mainardi constraints.

For a fixed two-dimensional metric, the problem of identifying equilibrium configurations that minimize the bending energy reduces to the problem of minimizing the *Willmore functional*

*H*is the mean curvature of the deformed surface, among isometric immersions of the given metric [16–18].

For a metric of constant positive Gaussian curvature, it is easy to show that the Willmore functional is minimized precisely for spherical solutions. That is, a flat nematic sheet whose director field encodes constant positive curvature *K* is expected to form part of a sphere of radius $R=1/\sqrt{K}$ upon stimulation, assuming that the sheet is small enough to exclude the possibility of self-intersection [6].

In the case of a metric of constant negative Gaussian curvature, identifying minimizers of the Willmore functional is considerably less straightforward. In [19], it is shown that for a hyperbolic elastic disc that has already undergone local deformations, surfaces that are geodesic discs lying on hyperboloids of revolution of constant Gaussian curvature are minimizers of the Willmore functional among *smooth* immersions of the metric. These solutions will appear as saddle shapes in experiments and are expected to be energetically favourable for sufficiently small discs. However, it has been shown numerically that certain non-smooth wavy surfaces formed as odd periodic extensions of subsets of so-called Amsler surfaces are energetically more favourable than the smooth saddle shapes that correspond to discs lying on hyperboloids of revolution when the radius of the hyperbolic disc is sufficiently large [20]. We synthesize LC solid films with spatially programmed directors in order to realize shape-changing surfaces of constant Gaussian curvature.

#### (c) Experimental investigations

The director profile in the plane of nematic LC solid sheets can be programmed through a variety of methods, including mechanical and magnetic fields [21,22]. Using these methods, however, it is difficult to spatially control the director orientation. Here we use chemistry amenable to surface alignment: precursor molecules, to what will become the nematic solid sheet, have been specifically designed to align to treated surfaces. Using this approach, low molar mass nematic LC monomers are filled between two plates separated by a well-defined gap. The treated surfaces, on the interior of the plates, direct the ordering of the LCs along a specific orientation through the thickness of the material. By using reactive nematic mesogens, this director orientation can be trapped in an elastic solid. Spatially complex director patterns in LC cells are prepared using point-by-point photoalignment of an azobenzene dye by irradiation with polarized light [12,23]. By altering the polarization of the incident light, the in-plane orientation of the director of the LC can be spatially controlled. The resulting director field is a pixelated approximation of the desired smooth pattern with each pixel measuring 100×100 μm.

A number of glassy liquid crystalline solids have been demonstrated to be compatible with surface alignment techniques. Here we use one such composition with λ=0.94 and *ν*=0.92 [24]. Specifically, we use the composition with the lowest cross-link density from this work. This composition is representative of the larger class of nematic LC glasses that can be aligned using surface alignment techniques [25]. The director patterns depicted in figure 1*a*,*c* were chosen to assess the viability of generating Gaussian curvature on exposure to stimulus. After fabrication, the LC solid film is flat at 25°*C* and retains the expected birefringence of an aligned nematic, as seen in figure 2*a*.

The thermally induced shape change of the nematic LC glass is shown in figure 2*b*. As predicted, the pattern depicted in figure 1*a* leads to the formation of positive Gaussian curvature, whereas the pattern from figure 1*c* leads to negative Gaussian curvature. On removal of the heat, the film returns to a largely flat state. It should be noted that the positive Gaussian curvature sample exhibits a periodic buckling around a pair of oppositely faced edges of the film. This is likely owing to the relatively sharp change in director angle with respect to the resolution of the patterning technique near these edges of the film. This buckling highlights the limitation on the curvature that can be achieved in nematic LC glasses with comparatively small strains.

In order to improve the quality of the surfaces that are formed in stimulated nematic glasses, a discoid subsection of the patterned films was removed and exposed to stimulus, as shown in figure 3. In the case of films encoded with either of the identified patterns, the predicted smooth curvature is realized upon stimulation. Indeed, shape selection of the equilibrium surface in the positive curvature case seems to be in remarkable qualitative agreement with the predicted solution of a spherical cap. In the negative curvature case too, the equilibrium surface appears largely consistent with the hyperboloid saddle solution that is predicted to be energetically favourable for a glassy film at this scale. Owing to the tendency for buckling in areas of the film where the director changes rapidly with respect to the resolution of the patterning process, it should be noted that the curvature cannot be increased by simply scaling the pattern to smaller dimensions. Instead, higher strain materials are needed.

To facilitate larger curvature realization, we prepare a comparatively high strain surface-alignable LC elastomer with λ=0.65 [23]. Using the pattern depicted in figure 1*a*, positive Gaussian curvature is encoded in the elastomeric film. As can be seen in figure 4*a*, the film encoded for positive Gaussian curvature forms part of a sphere with a slightly elliptical distortion. Figure 4*a* also shows a complexly buckled hyperbolic surface that is formed when an elastomeric disc encoded to give negative Gaussian curvature is exposed to stimulus. For smaller diameter films encoded with negative Gaussian curvature, a classic saddle shape can be observed, as shown in figure 4*b*. The surface that is formed by the larger radius hyperbolic disc of figure 4*b* might be interpreted as a distorted periodic Amsler surface. Comparatively, these deformations are significantly larger than those observed for the glassy films despite being more than three times as thick. Understanding the full spectrum of shape selection for films encoded with negative Gaussian curvature is an area of ongoing consideration.

These qualitative explorations can be made quantitative if (i) the λ values can be accurately known (results are very sensitive to λ), and (ii) accurate imaging from several vantage points is made, so that shell shapes and dimensions can be accurately fitted. This work, also requiring further synthesis, is underway.

### 3. Spiral director fields

#### (a) Metric analysis in polar coordinates

Now, consider a director field $\mathit{n}={n}_{1}{\hat{\mathit{e}}}^{1}+{n}_{2}{\hat{\mathit{e}}}^{2}$ whose components are given by

*r*,

*θ*) are the polar coordinates on

*ω*, and $\mathrm{\Psi}(r,\theta )=\psi (r\mathrm{cos}\theta ,r\mathrm{sin}\theta )$ is the alignment angle field expressed in (

*r*,

*θ*) coordinates. The metric components with respect to polar coordinates are

*Ψ*, the metric components are

For nematic angle fields of the form *Ψ*(*r*,*θ*)=*θ*+*α*(*r*), the metric components in polar coordinates depend only on the radius *r*, and the Gaussian curvature is given by

*α*is the angle that the director makes with the radial direction. For logarithmic spiral patterns, this angle is constant, so that, by (3.6), the Gaussian curvature generated by a director field whose integral curves consist of logarithmic spirals vanishes everywhere except at the point defect at the origin. Figure 5

*a*shows a logarithmic spiral pattern with

*α*=

*π*/4.

By prescribing *α*=*α*(*r*), we can obtain director fields that generate radially symmetric Gaussian curvature distributions *K*=*K*(*r*). An example with *α*(*r*)=*r*/*L* is shown in figure 5*b*. The Gaussian curvature, in this case, is given by

For constant $K\in \mathbb{R}$, (3.6) can be viewed as an ODE in *α* which is solved by

*C*(

*K*)=

*K*/(λ

^{−2}−λ

^{2ν}), and

*c*

_{1},

*c*

_{2}are real constants of integration. Depending on the choice of

*c*

_{1}and

*c*

_{2}, these solutions define curvature-inducing nematic patterns on compact discs or annular domains. If

*K*>0,

*c*

_{2}=0 and

*c*

_{1}≤1, then the solution (3.8) defines a spiral pattern on the compact disc

*α*in (3.8) attains its maximum value of

*α*=

*π*/2. For

*c*

_{1}=1, we obtain a spiral pattern on a disc of maximal radius, which encodes constant positive Gaussian curvature

*K*. This gives the largest disc on which a smooth radially symmetric director field can be used to encode a prescribed constant Gaussian curvature. Attempting to increase the radius beyond this point by varying

*c*

_{1}as a parameter leads to a bifurcation at which the solution domain becomes annular instead of discoid. If

*K*<0, then it is

*c*

_{2}=0 and

*c*

_{1}=−1 which yield the spiral pattern of maximal radius. These patterns are shown in figure 6

*a*,

*b*.

In the light of the result that the orthogonal dual of a director field has precisely the opposite Gaussian curvature distribution, we find that the spirals in figure 6 make perfect sense. In particular, the spiral in figure 6*b* is equivalent to the spiral that turns in the opposite direction, defined by an $\alpha (r)=-\frac{1}{2}\mathrm{arccos}(-C(K){r}^{2}/2-1)$. This equivalent pattern is the orthogonal dual of the positive curvature spiral of figure 6*a* and is illustrated in figure 6*c*, noting *C*(*K*) has reversed sign.

Finally, note that for *K*=0 and the choice of *c*_{2}=0, solution (3.8) reduces to a logarithmic spiral pattern defined by *α*=constant. It should be noted, however, that the resulting surfaces upon stimulation will typically not remain planar. In particular, such patterns are expected to form cone and anti-cone equilibrium configurations upon stimulation [4,7,26,27].

#### (b) Logarithmic spiral patterns: cone/anti-cone formation

It is shown in [7,4] that azimuthal (*α*=*π*/2) and radial (*α*=0) director fields generate cone and anti-cone surfaces upon stimulation, respectively. It is also argued there that logarithmic spiral patterns result in cone and anti-cone surfaces depending on whether the angle *α* is greater than or less than some threshold angle depending on λ and *ν*. Here, we will rederive this result using the nematic metric in polar coordinates. Generally, the angles between curves on a patterned nematic sheet become distorted upon stimulation. In particular, the radial curves on the initial sheet will not in general be orthogonal to the images of the azimuthal curves upon stimulation. We note that two curves *Γ*_{1} and *Γ*_{2} with polar parametric representations (*r*,*θ*)=(*r*_{1}(*t*),*θ*_{1}(*t*)) and (*r*,*θ*)=(*r*_{2}(*t*),*θ*_{2}(*t*)) are orthogonal with respect to the metric d*s*^{2}=*a*_{rr} d*r*^{2}+2*a*_{rθ} d*r* d*θ*+*a*_{θθ} d*θ*^{2}, if and only if

*r*

_{1}′=0 and

*θ*

_{1}′=1. Thus, the curve (

*r*(

*t*),

*θ*(

*t*)) on the undeformed sheet whose image becomes orthogonal to the images of the concentric circles upon stimulation satisfies

*a*

_{θr}

*r*′+

*a*

_{θθ}

*θ*′=0. That is, d

*θ*/d

*r*=−

*a*

_{θr}/

*a*

_{θθ}. The length of such a curve emanating from the origin is given by

*θ*/d

*r*=−

*a*

_{θr}/

*a*

_{θθ}, (3.11) reduces to

*r*deforms into following exposure to stimulus is given by

*Ψ*=

*θ*+

*α*(

*r*), where ${a}_{\theta \theta}={r}^{2}[{\mathrm{\lambda}}^{-2\nu}-({\mathrm{\lambda}}^{-2\nu}-{\mathrm{\lambda}}^{2}){\mathrm{sin}}^{2}\alpha ]$ by (3.4).

For logarithmic spirals, a cone is expected to form if 2*πl*_{1}>*l*_{2} and an anti-cone is expected to form if 2*πl*_{1}<*l*_{2}. When a cone is formed, the cone angle is

*πl*

_{1}=

*l*

_{2}, the sheet is expected to remain flat (

*ϕ*=

*π*/2) upon stimulation and the threshold angle

*α*

_{c}at which this holds satisfies ${\mathrm{\lambda}}^{-2\nu}-({\mathrm{\lambda}}^{-2\nu}-{\mathrm{\lambda}}^{2}){\mathrm{sin}}^{2}{\alpha}_{\mathrm{c}}={\mathrm{\lambda}}^{1-\nu}$. That is,

*α*of a logarithmic director field such that the lengths of azimuthal curves remain fixed. This is achieved by requiring $\sqrt{{a}_{\theta \theta}(r)}=r$, which is satisfied from (3.4) when

*α*

_{c}<

*α*

_{0}, then a cone is formed and as a result the symmetry of revolution is unbroken. This condition is equivalent to λ

^{1−ν}>1; i.e.

*ν*>1. In particular, a logarithmic spiral director field with

*α*=

*α*

_{0}patterned across the surface of a disc will keep the circular boundary fixed, whereas the interior of the surface rises to form a cone provided

*ν*>1. The cone angle in this case is $\varphi =\mathrm{arcsin}({\mathrm{\lambda}}^{\nu -1})$.

#### (c) Spherical spindles

Here we extend the analysis presented in §3b to study the formation of spherical caps and spherical spindles using the positive-curvature-inducing spiral patterns defined by

*c*. Spherical spindles [28,29] are surfaces of revolution of constant positive Gaussian curvature

*K*. The top half of a spherical spindle can be parametrized as

**(**

*γ**s*)=(

*γ*

_{1}(

*s*),

*γ*

_{2}(

*s*)) is defined by

*a*illustrates one such surface and its profile curve. Note that for $\rho =1/\sqrt{K}$, the spindle reduces to a sphere of radius $1/\sqrt{K}$.

The length of each curve from the apex that intersects the azimuthal circles orthogonally is

*γ*

_{1}′

^{2}+

*γ*

_{2}′

^{2}=1. Using (3.12) we also have

*a*

_{θθ}=

*r*

^{2}(

*μ*

_{1}−

*μ*

_{2}

*r*

^{2}), where

*r*(

*s*):

*r*centred at the origin in the undeformed sheet becomes ${l}_{2}=2\pi r\sqrt{{\mu}_{1}-{\mu}_{2}{r}^{2}}$ upon stimulation. Thus, we find that

*l*

_{2}with 2

*πγ*

_{1}, the spindle radius being (3.19), we identify ${\mathrm{\lambda}}^{\nu -1}{\mu}_{1}/\sqrt{K}$ with

*ρ*. Given $\rho \le 1/\sqrt{K}$, the deformed sheet is consistent with a spherical spindle precisely when λ

^{ν−1}

*μ*

_{1}<1. That is, a spindle is expected to form if

*c*=1−2/(1+λ

^{1+ν}), then a spherical cap of radius $1/\sqrt{K}$ will form instead.

For any given choice of the parameter *c*, the spiral in figure 6*a* keeps the circumference of an azimuthal circle of radius *r*=*r*_{0} fixed provided that it satisfies $\sqrt{{a}_{\theta \theta}(r)}=r$; see (3.13). That is, *r*_{0} satisfies $\sqrt{{\mu}_{1}-{\mu}_{2}{r}_{0}^{2}}=1,$ so that

*c*=1−2/(1+λ

^{1+ν}), corresponding to a spherical cap, we obtain

*ν*>1 for

*r*

_{0}to be well-defined and positive. We have thus identified a director field on a thin disc that will cause it to form a spherical cap with the same fixed circular boundary upon exposure to stimulus, as shown in figure 7

*b*–

*d*. The zenith angle

*φ*

_{z}of the resulting spherical cap is

*ν*=2, the zenith angle of the resulting spherical cap can be as large as approximately 16°. So we see that this method allows for rather prominent spherical caps to form even in the presence of the small strains observed in nematic glasses. The same methods can be used to specify patterns that generate spindles with unchanged circular boundaries.

#### (d) Hyperbolic cones

We can perform a similar analysis in the case of the negative-curvature-inducing spiral patterns. The family of surfaces that form the negative curvature analogue of spherical spindles are hyperbolic cones [28], which are cone-like surfaces of revolution of constant negative Gaussian curvature away from the non-smooth tip of the cone. A hyperbolic cone of constant negative Gaussian curvature *K*<0 admits a parametrization

**(**

*γ**s*)=(

*γ*

_{1}(

*s*),

*γ*

_{2}(

*s*)) is defined by

*γ*

_{2}can also be expressed in terms of the incomplete elliptic integral of the second kind

*E*, as

*a*illustrates a typical hyperbolic cone and its profile curve arising from a negative curvature spiral.

Using the spiral pattern specified by (3.17) with prescribed constant *K*<0, the analysis proceeds as in the case of positive Gaussian curvature. The expression in (3.24) is modified to

*μ*

_{1}and

*μ*

_{2}are as in (3.23) with

*K*<0. This yields $r(s)=\sqrt{{\mu}_{1}/|{\mu}_{2}|}\mathrm{sinh}(\frac{1}{2}\sqrt{|K|}s)$, and thus

*γ*

_{1}as a spindle radius to generate a circumference, if we identify ${\mathrm{\lambda}}^{\nu -1}{\mu}_{1}/\sqrt{|K|}$ with

*ρ*, whereby we note that the deformed sheet is consistent with a hyperbolic cone precisely if λ

^{ν−1}

*μ*

_{2}<1; i.e. if

*c*<1−2/(1+λ

^{1+ν}).

For any given choice of the parameter *c*, the spiral in figure 6*b* keeps the circumference of an azimuthal circle of radius *r*=*r*_{0} fixed provided that it satisfies $\sqrt{{a}_{\theta \theta}(r)}=r$. That is, *r*_{0} must satisfy $\sqrt{{\mu}_{1}+|{\mu}_{2}|{r}_{0}^{2}}=1$, so that

*r*

_{0}must also satisfy ${r}_{0}\le {r}_{max}(c)$, where

*K*<0. Figure 8

*b*shows the negative curvature spiral pattern defined by (3.17) with

*c*=−1. The solid (red) circle superimposed on the director field indicates the circle of radius

*r*

_{0}whose length is unchanged. The pattern in figure 8

*c*is thus expected to generate a hyperbolic cone of fixed boundary upon stimulation.

#### (e) Annular extensions of curvature-inducing spiral patterns

Finally, we note that it is possible to radially extend the constant curvature encoding spiral patterns of figure 6 to cover larger domains by the use of director fields defined on an annulus whose inner radius matches that of the disc. We can do so consistently through judicious choices of *c*_{1} and *c*_{2}≠0 in the solution (3.8) for *α*, to ensure that the directors agree on the circular boundary between the disc and the first annulus. We can then repeat the process and extend the resulting pattern further by attaching another solution (3.8) defined on a second annulus whose inner radius matches the outer radius of the first annulus using a suitable new choice of *c*_{1} and *c*_{2}. This process can be continued to produce director fields that encode constant Gaussian curvature on circular domains of any desired size. The resulting Gaussian curvature is well defined at all points on the domain with the exception of the boundaries where the different textures meet. Figure 9 depicts the nematic director field of the extended pattern and shows how the angle function *α*=*α*(*r*) varies under such an extension scheme for the *K*>0 spiral. The orthogonal dual of the extended pattern provides the corresponding result in the case of negative Gaussian curvature.

### 4. Discussion

In general, surfaces morphing to Gaussian-curved forms offer rich possibilities, most especially for delivering work as microdevices. Blocking the stroke of such morphing surfaces creates stretch rather than bend. Normally, extensile spontaneous strains of thin structures, when blocked lead to Euler instabilities and potentially strong processes become weaker bend actions. A Gaussian-curved structure is less likely to find such soft alternative routes. Impeding motion of the tip of an emerging cone, and thus resisting the contraction of circumferences lead to large associated forces, and these thus consequently do work, or perhaps pump an enclosed fluid. One might term this action ‘strong actuation’ as opposed to a weak actuation associated with bend. Another advantage is that which is often a short stroke (associated e.g. with a 4% length change in a glass) is larger in a morphing shell, for instance in the rise of a cone.

To avoid the nugatory effects of soft elasticity relieving strong tensile forces during actuation by director rotation, one must only have such stresses along the director. Often this can be achieved, for instance when a cone rising on heating is blocked, the circumferential stresses arising cannot be negated by director rotation, because stresses are already aligned along the director. However, more complex director distributions may be less reliable and it would be safer to use elastomers heated to isotropy, or glasses where there is no director rotation.

The anchoring of deforming shells to unresponsive mounts surrounding them is of vital importance for exploiting shells with tuneable Gaussian curvature. We have given explicit forms for the director distribution to achieve deformations, of both signs of curvature, elastically compatible with rigid connections to a mounting at the shell boundary. Negative curvature surfaces might be of practical advantage because they have sharp tips and are thus more acutely sensed, for instance, in haptic devices. Further, the slopes of their profile curves decrease with increasing radius, in contrast to cones. Thus, at their fixed boundaries where they could be mounted, there is less of a difference in angle and thus lower elastic energy in the attachment region. It is, however, impossible to design the pattern so as to ensure that the radius at which the boundary is fixed precisely matches the radius at which the profile curve fully flattens.

In summary, our results clearly indicate that Gaussian curvature can be realized in both low-strain, high-modulus glassy and high-strain, low-modulus elastomeric LC solid surfaces using appropriate, smooth, in-plane director fields patterned across initially flat films. Our preliminary experimental investigations show that patterned nematic solids may indeed be promising candidates for applications. In particular, the observed shape transformations were in qualitative agreement with theoretical predictions and the behaviour of the films in response to stimulus was found to be robust and reproducible. We hope that our results will encourage and stimulate further experimental research in achieving desired shape transitions in nematic LC solid surfaces. In particular, our discovery of anchorable structures should lead to further work in assessing the viability of morphing shells for specific applications.

### Data accessibility

All experimental data in this work (photographs, experimental conditions, etc.) are given in the text. All theoretical/computational data are also given in the text and relate to plots of mathematical results derived analytically.

### Authors' contributions

C.M. and M.W. conceived of the metric mechanism for mechanical response, and the route to boundary anchoring of raised structures. C.M. performed the mathematical calculations, partially in consultation with M.W., and performed all the computations involved. C.M. wrote the paper, with sections provided by M.W., T.H.W. and T.J.W. All experimental work was performed by T.H.W. and T.J.W. All authors approved the final version of the paper.

### Competing interests

We have no competing interests.

### Funding

C.M. is supported by the Engineering and Physical Sciences Research Council of the United Kingdom. T.J.W. and T.H.W. acknowledge financial support from the Materials and Manufacturing Directorate of the Air Force Research Laboratory and the Air Force Office of Scientific Research.

## Footnotes

1 Indeed, λ(*r*) variation, but in isotropic systems without a director, is the basis of Gaussian curvature by spatially dependent swelling [2,3].

2 Heating the orthogonal dual to give a λ<1 along its director corresponds to cooling the original pattern to give an elongation along its director of *Λ*=λ^{−ν}>1. Thus, heating and cooling give converse Gaussian curvatures.

### References

- 1
Dervaux J, Ben Amar M . 2008Morphogenesis of growing soft tissues.**Phys. Rev. Lett.**, 068101. (doi:10.1103/PhysRevLett.101.068101) Crossref, PubMed, ISI, Google Scholar**101** - 2
Klein Y, Efrati E, Sharon E . 2007Shaping of elastic sheets by prescription of non-euclidean metrics.**Science**, 1116–1120. (doi:10.1126/science.1135994) Crossref, PubMed, ISI, Google Scholar**315** - 3
Kim J, Hanna JA, Byun M, Santangelo CD, Hayward RC . 2012Designing responsive buckled surfaces by halftone gel lithography.**Science**, 1201–1205. (doi:10.1126/science.1215309) Crossref, PubMed, ISI, Google Scholar**335** - 4
Modes C, Bhattacharya K, Warner M . 2011Gaussian curvature from flat elastica sheets.**Proc. R. Soc. A**, 1121–1140. (doi:10.1098/rspa.2010.0352) Link, Google Scholar**467** - 5
Aharoni H, Sharon E, Kupferman R . 2014Geometry of thin nematic elastomer sheets.**Phys. Rev. Lett.**, 257801. (doi:10.1103/PhysRevLett.113.257801) Crossref, PubMed, ISI, Google Scholar**113** - 6
Mostajeran C . 2015Curvature generation in nematic surfaces.**Phys. Rev. E**, 062405. (doi:10.1103/PhysRevE.91.062405) Crossref, ISI, Google Scholar**91** - 7
Modes CD, Warner M . 2011Blueprinting nematic glass: systematically constructing and combining active points of curvature for emergent morphology.**Phys. Rev. E**, 021711. (doi:10.1103/PhysRevE.84.021711) Crossref, ISI, Google Scholar**84** - 8
Van Oosten C, Harris K, Bastiaansen CW, Broer D . 2007Glassy photomechanical liquid-crystal network actuators for microscale devices.**Eur. Phys. J. E, Soft Matter Biol. Phys.**, 329–336. (doi:10.1140/epje/i2007-10196-1) Crossref, PubMed, ISI, Google Scholar**23** - 9
Warner M, Terentjev EM . 2003**Liquid crystal elastomers**,**vol. 120**. Oxford, UK: Oxford University Press. Google Scholar - 10
de Haan LT, Sánchez-Somolinos C, Bastiaansen CM, Schenning AP, Broer DJ . 2012Engineering of complex order and the macroscopic deformation of liquid crystal polymer networks.**Angew. Chem. Int. Ed.**, 12 469–12 472. (doi:10.1002/anie.201205964) Crossref, ISI, Google Scholar**51** - 11
Ware TH, Perry ZP, Middleton CM, Iacono ST, White TJ . 2015Programmable liquid crystal elastomers prepared by thiol–ene photopolymerization.**ACS Macro Lett.**, 942–946. (doi:10.1021/acsmacrolett.5b00511) Crossref, ISI, Google Scholar**4** - 12
McConney ME, Martinez A, Tondiglia VP, Lee KM, Langley D, Smalyukh II, White TJ . 2013Topography from topology: photoinduced surface features generated in liquid crystal polymer networks.**Adv. Mater.**, 5880–5885. (doi:10.1002/adma.201301891) Crossref, PubMed, ISI, Google Scholar**25** - 13
Tabiryan NV, Nersisyan SR, Steeves DM, Kimball BR . 2010The promise of diffractive waveplates.**Opt. Photon. News**, 40–45. (doi:10.1364/OPN.21.3.000040) Crossref, Google Scholar**21** - 14
D'Ambrosio V *et al.*2013Photonic polarization gears for ultra-sensitive angular measurements.**Nat. Commun.**, 2432. (doi:10.1038/ncomms3432) Crossref, PubMed, ISI, Google Scholar**4** - 15
Ciarlet PG . 2005An introduction to differential geometry with applications to elasticity.**J. Elast.**, 1–215. (doi:10.1007/s10659-005-4738-8) Crossref, Google Scholar**78** - 16
Lewicka M, Reza Pakzad M . 2011Scaling laws for non-euclidean plates and the w2, 2 isometric immersions of riemannian metrics.**ESAIM: Control, Optim. Calc. Var.**, 1158–1173. (doi:10.1051/cocv/2010039) Crossref, ISI, Google Scholar**17** - 17
Efrati E, Sharon E, Kupferman R . 2011Hyperbolic non-euclidean elastic strips and almost minimal surfaces.**Phys. Rev. E**, 046602. (doi:10.1103/PhysRevE.83.046602) Crossref, ISI, Google Scholar**83** - 18
- 19
Gemmer JA, Venkataramani SC . 2011Shape selection in non-euclidean plates.**Physica D: Nonlinear Phenomena**, 1536–1552. (doi:10.1016/j.physd.2011.07.002) Crossref, ISI, Google Scholar**240** - 20
Gemmer J, Venkataramani SC . 2013Shape transitions in hyperbolic non-euclidean plates.**Soft Matter**, 8151–8161. (doi:10.1039/C3SM50479D) Crossref, ISI, Google Scholar**9** - 21
Küpfer J, Finkelmann H . 1991Nematic liquid single crystal elastomers.**Die Makromolekulare Chemie, Rapid Commun.**, 717–726. (doi:10.1002/marc.1991.030121211) Crossref, Google Scholar**12** - 22
Schuhladen S, Preller F, Rix R, Petsch S, Zentel R, Zappe H . 2014Iris-like tunable aperture employing liquid-crystal elastomers.**Adv. Mater.**, 7247–7251. (doi:10.1002/adma.201402878) Crossref, PubMed, ISI, Google Scholar**26** - 23
Ware TH, McConney ME, Wie JJ, Tondiglia VP, White TJ . 2015Voxelated liquid crystal elastomers.**Science**, 982–984. (doi:10.1126/science.1261019) Crossref, PubMed, ISI, Google Scholar**347** - 24
Wie JJ, Lee KM, Ware TH, White TJ . 2015Twists and turns in glassy, liquid crystalline polymer networks.**Macromolecules**, 1087–1092. (doi:10.1021/ma502563q) Crossref, ISI, Google Scholar**48** - 25
Liu D, Broer DJ . 2014Liquid crystal polymer networks: preparation, properties, and applications of films with patterned molecular alignment.**Langmuir**, 13 499–13 509. (doi:10.1021/la500454d) Crossref, ISI, Google Scholar**30** - 26
Müller MM, Ben Amar M . 2008Conical defects in growing sheets.**Phys. Rev. Lett.**, 156104. (doi:10.1103/PhysRevLett.101.156104) Crossref, PubMed, ISI, Google Scholar**101** - 27
Modes C, Warner M . 2015Negative gaussian curvature from induced metric changes.**Phys. Rev. E**, 010401. (doi:10.1103/PhysRevE.92.010401) Crossref, ISI, Google Scholar**92** - 28
Abbena E, Salamon S, Gray A . 2006**Modern differential geometry of curves and surfaces with Mathematica**. New York, NY: CRC Press. Google Scholar - 29
Modes C, Warner M . 2012Responsive nematic solid shells: topology, compatibility, and shape.**Europhys. Lett.**, 36007. (doi:10.1209/0295-5075/97/36007) Crossref, Google Scholar**97**