PHYSICAL REVIEW E
96
, 032701 (2017)
Electroclinic effect in chiral smectic-
A
liquid crystal elastomers
Noy Cohen
*
and Kaushik Bhattacharya
Division of Engineering and Applied Science, California Institute of Technology, Pasadena, California 91125, USA
(Received 30 June 2017; published 18 September 2017)
Chiral smectic-
A
liquid crystal elastomers are rubbery materials composed of a lamellar arrangement of
liquid crystalline mesogens. It has been shown experimentally that these materials shear when subjected to an
electric field due to the electrically induced tilt of the director. Experiments have also shown that shearing a
chiral smectic-
A
elastomer gives rise to a polarization. Roughly, the shear force tilts the directors which, in
turn, induce electric dipoles. This paper builds on previous works and models the electromechanical response of
smectic-
A
elastomers using free energy contributions that are associated with the lamellar structure, the relative
tilt between the director and the layer normal, and the coupling between the director and the electric field. To
illustrate the merit of the proposed model, two cases are considered—a deformation induced polarization and
an electrically induced deformation. The predictions according to these two models qualitatively agree with
experimental findings. Finally, a cylinder composed of helical smectic layers is also considered. It is shown that
the electromechanical response varies as a function of the helix angle.
DOI:
10.1103/PhysRevE.96.032701
I. INTRODUCTION
Smectic-
A
liquid crystal elastomers are rubbery materials
composed of mesogens that are uniformly oriented and
organized in a layered structure. In their pioneering work,
Garoff and Meyer [
1
,
2
] found that low molecular weight
smectic-
A
∗
liquid crystals comprising chiral molecules exhibit
an electroclinic effect, i.e., there is a direct coupling between
the molecular tilt and an applied electric field. In smectic
elastomers comprising chiral mesogens, an electric field acting
along the layer plane induces a shear deformation. This
phenomenon directly emanates from the chirality of the meso-
gens. Specifically, the electric field couples to the transverse
dipole of the chiral mesogens and induces a tilt between
the mesogens and the layer normal. As a consequence, the
elastomer shears and contracts along the direction of the layer
normal. Typically, such elastomers are incompressible and
subsequently this deformation is accompanied by an expansion
along the transverse directions. Another consequence of the
chirality of the mesogens was demonstrated by Lehmann
et al
.
[
3
], where it was found that flipping the direction of the electric
field reverses the direction of the mechanical tilt.
Several experimental works were carried out to study the
response of a smectic-
A
elastomer subjected to an electric
field. Spillmann, Ratna, and Naciri [
4
] mounted a thin sheet
of elastomer in a tension clamp with the smectic layers
either perpendicular or parallel to the tension clamps and
applied an electric field normal to the sheet. An electric
field
E
=
16
.
7MV
/
m was applied to a film the layers of
which are parallel to the clamps, and a contraction of
∼
1%
was measured. Next, the sample was placed such that the
layers are perpendicular to the clamps. Upon application of
the electric field, an expansion of
∼
8
.
2% was measured.
Hiraoka
et al.
[
5
,
6
] prepared a smectic elastomer film with
an initial thickness
t
0
=
600
μ
m and applied an electric field
E
=
1MV
/
M. They measured an electrically induced tilt
angle of
∼
4
◦
,
a contraction along the layer normal of
∼
0
.
2%,
*
noyco@caltech.edu
and a shear angle of
∼
0
.
08
◦
. A recent experiment by Spillmann
et al
.[
7
] demonstrated that smectic elastomers can also twist
in response to an electric field.
Previous works have also investigated the purely me-
chanical behavior of smectic elastomers. The experimental
work of Nishikawa, Finkelmann, and Brand [
8
], Nishikawa
and Finkelmann [
9
], and Spillmann, Ratna, and Naciri [
4
]
demonstrated that the mechanical modulus of smectic-
A
elastomers parallel to the layers normal is significantly higher
than that along the layer plane. This anisotropy emanates from
the lamellar microstructure of the elastomer. Interestingly, the
experiments of Nishikawa and Finkelmann [
9
] revealed that
the elastomer becomes significantly softer above a threshold
of
∼
3% strain along the layers normal.
On the theoretical side, Adams and Warner [
10
] developed
an energy-based model for smectic-
A
elastomers which as-
sumes that the layer normal and the director are locked
and may only rotate in response to an external force. This
assumption prevents a tilt and a smectic-
C
like ordering. The
work of Stenull and Lubensky [
11
] argued that the director is
not rigidly coupled to the layer normal and, as a consequence, a
shear deformation can rotate the smectic layers and the director
with respect to the layer normal. The later works of Adams
et al
.[
12
] and Stenull
et al
.[
13
] provided another model
that accounts for the stiffness of the smectic layers and the
resistance to the tilt by adding appropriate contributions to the
free energy.
From an electric viewpoint, the work of Corbett and Warner
[
14
] theoretically examined the electromechanical behavior of
a nematic elastomer subjected to electric fields. In this paper, it
was shown that the mechanical and electrical anisotropies can
cause a thin film to expand in the direction of the electric field
rather than contract. As a result of the Poisson effect, the film
shrinks along the two directions in the plane of the electrodes.
Selinger
et al
.[
15
] found that under high electric fields the
state of uniform induced tilt can become unstable and lead to
chiral modulations. In their work, they modeled and simulated
this phenomenon.
The aim of this paper is to model the coupled behavior
of smectic-
A
elastomers under finite deformations resulting
2470-0045/2017/96(3)/032701(7)
032701-1
©2017 American Physical Society
NOY COHEN AND KAUSHIK BHATTACHARYA
PHYSICAL REVIEW E
96
, 032701 (2017)
from the electroclinic effect. To this end, we characterize the
chirality of the mesogens by a chirality vector which couples
to the tilt of the director. The proposed model, presented in
Sec.
II
, is based on the minimization of the free energy of the
smectic elastomer. To demonstrate the behavior predicted by
the model, two cases are examined in Sec.
III
—an electrically
induced deformation and a deformation induced polarization.
In the first case, an electric field is applied to the smectic-
A
elastomer and the resulting deformation is determined through
the minimization of the free energy. Subjected to this loading,
the tilt angle is larger than the shear angle. This trend is in
agreement with the experimental findings of Hiraoka
et al
.
[
6
]. We also show that the orientation of the chirality of the
mesogens can either enhance or diminish the polarization.
In the second case, we follow the experimental work of
Kramer and Finkelmann [
16
] and consider the response of a
smectic-
A
elastomer subjected to a simple shear deformation.
The proposed model predicts that the prescribed deformation
rotates the director with respect to the layer normal and
therefore induces a tilt of the mesogens. Subsequently, the
inherent chirality gives rise to spontaneous dipoles which, in
turn, induce a polarization in the elastomer. In accordance
with experimental findings, the proposed model predicts that
the shear angle is larger than the tilt angle. A cylindrical setup
of a smectic-
A
elastomer with helical layers is considered
in Sec.
IV
. The electromechanical response of an elastomer
with various pitches and helix angles is computed. The main
conclusions are summarized in Sec.
V
.
II. A COUPLED ELECTROMECHANICAL MODEL
Consider a smectic-
A
elastomer that deforms from a refer-
ence to a current configuration due to a coupled electromechan-
ical loading. The material points are denoted
x
and current po-
sitions are denoted as
y
(
x
), so that the deformation gradient is
defined as
F
=∇
y
(
x
), where
∇
denotes the gradient operation
with respect to the referential coordinate system. The smectic
order is characterized by two order parameters—the director
and the chiral vector. These are denoted
ˆ
n
0
and
q
0
in the ref-
erence configuration and
ˆ
n
and
q
in the current configuration.
Note that bold letters with hats on top denote unit vectors. The
chirality vector is perpendicular to the director such that
q
·
ˆ
n
=
0. The layer normal at the current configuration is given as
ˆ
k
=
1
−
q
2
ˆ
n
−
ˆ
n
×
q
,
(1)
where
q
is the magnitude of the chirality vector
q
. Note that
q
=
0
gives
ˆ
k
=
ˆ
n
and characterizes smectic-
A
while
q
=
0
gives rise to smectic-
C
. This smectic order gives rise to a spon-
taneous elongation along the director given by the step-length
tensor
=
r
−
1
/
3
[
I
+
(
r
−
1)
ˆ
n
⊗
ˆ
n
], where
r
is a parameter
related to the strength of the ordering [
12
]. It also gives rise to
a spontaneous polarization
P
s
=
C
q
, where
C
is a parameter.
We assume that the free energy-density function of an
elastomer subjected to a uniform electric field by surface
electrodes can be viewed as the sum of five contributions:
W
=
W
tr
+
W
l
+
W
t
+
W
e
+
W
m
.
(2)
The first contribution is the classical elastic free energy of an
anisotropic elastomer [
17
]:
W
tr
=
μ
2
Tr(
F
0
F
T
−
1
)
,
(3)
where
μ
is the shear modulus and
0
=
r
−
1
/
3
[
I
+
(
r
−
1)
ˆ
n
0
⊗
ˆ
n
0
] is the step-length tensor before the
deformation. Here,
I
is the identity tensor.
The second contribution to the free energy in Eq. (
2
)
concerns the change in the layer thickness:
W
l
=
B
2
1
|
F
−
T
ˆ
k
0
|
1
−
q
2
−
1
2
,
(4)
where
B
is a measure of the stiffness of the layers. This is
identical to the free energy contribution in the works of Adams
and Warner [
10
] and Adams
et al
.[
12
], but written in terms
of the magnitude of the physical chirality vector
q
. Above,
F
T
and
F
−
1
denote the transpose and inverse of
F
, respectively.
This formulation is advantageous because it is indifferent to
the sign of the director
ˆ
n
and therefore satisfies the objectivity
requirements associated with the energy-density function.
The third term in Eq. (
2
) concerns the change in chirality,
or equivalently the tilt (or rotation) of the director with respect
to the layer normal:
W
t
=
A
2
q
2
,
(5)
where
A
is the resistance of the mesogens to tilt. This is again
the same as Adams and Warner [
10
] and Adams
et al
.[
12
], but
written in terms of the magnitude of the chirality vector. Note
that Eq. (
5
) assumes that the ground state is smectic-
A
,butit
can be modified if the ground state is smectic-
C
.
The fourth term in Eq. (
2
) pertains to the energetic
contribution of the electrical sources to the total energy (see,
for example, Shu and Bhattacharya [
18
]). Specifically,
W
e
=−
ε
0
2
E
·
E
−
P
·
E
+
1
2
ε
0
χ
(
P
−
P
s
)
·
(
P
−
P
s
)
,
(6)
where
P
is the polarization,
χ
is the susceptibility, and
ε
0
is
the permittivity in vacuum.
To determine the polarization and the deformation under
a prescribed coupled electromechanical loading, the energy-
density function is minimized with respect to
P
,
F
, and the
chirality vector
q
. By computing
∂W
∂
P
=
0, we find that
P
=
ε
0
χ
E
+
P
s
=
ε
0
χ
E
+
C
q
.
(7)
The last term in Eq. (
2
) concerns the external mechanical
loading:
W
m
=−
S
·
F
,
(8)
where
S
is the first Piola-Kirchhoff stress tensor.
For convenience, we normalize the energy-density function
by
μ
2
and substitute Eq. (
7
) into Eq. (
2
) to obtain
w
=
2
μ
W
=
Tr(
F
0
F
T
−
1
)
+
b
1
|
F
−
T
ˆ
k
0
|
1
−
q
2
−
1
2
+
aq
2
−
ε
0
(1
+
χ
)
μ
E
·
E
−
2
C
μ
q
·
E
−
2
S
μ
·
F
,
(9)
032701-2
ELECTROCLINIC EFFECT IN CHIRAL SMECTIC-
A
. . .
PHYSICAL REVIEW E
96
, 032701 (2017)
where
b
=
B
μ
and
a
=
A
μ
are the normalized layer stiffness and
the normalized resistance to tilt, respectively.
The constants in the proposed model can be determined
through experiments. The order of magnitude of the shear
modulus and the layer stiffness are
μ
∼
0
.
1–1 MPa and
B
∼
10 MPa, respectively [
12
]. Fitting to experimental findings
resulted in a resistance to tilt parameter
A
∼
0
.
1–1 MPa
[
19
,
20
]. The anisotropy parameter
r
has been measured
between 1
.
05 and 60 in elastomers composed of prolate chains
[
10
,
17
]. Lastly, we can deduce the order of magnitude of the
chirality related parameter from experimental findings. It is
found that
C
∼
10
−
3
–10
−
4
C
/
m
2
.
III. THE ELECTROMECHANICAL RESPONSE
OF SMECTIC-
A
ELASTOMERS
To examine the predictions of the proposed model, consider
a thin film of a smectic-
A
elastomer with an initial thickness
t
0
.
We prescribe a global coordinate system
{
ˆ
X
,
ˆ
Y
,
ˆ
Z
}
such that at
the reference configuration the director and the layer normal
are along the
ˆ
Z
direction, i.e.,
ˆ
n
0
=
ˆ
Z
and
ˆ
k
0
=
ˆ
Z
(see Fig.
1
).
The faces of the film are covered with flexible electrodes such
that an applied voltage induces an electric field
E
=
V
λ
E
t
0
ˆ
X
,
where
λ
E
is the stretch along the thickness of the sheet.
Mechanically, an applied stress
s
=
S
ˆ
Z
·
ˆ
Y
shears the film.
The elastomer experiences the general deformation gradient
F
=
⎛
⎝
λ
E
00
0
λ
y
λ
yz
00
λ
z
⎞
⎠
,
(10)
where
λ
z
=
h
H
,
λ
y
=
l
L
,
λ
yz
=
λ
z
tan
φ
, and
φ
is the shear
angle (see Fig.
1
). Due to the incompressibility constraint,
λ
E
=
1
λ
y
λ
z
. As a result of an electromechanical loading, the
director rotates such that at the current configuration
ˆ
n
=
sin
θ
ˆ
Y
+
cos
θ
ˆ
Z
. The chirality vector is
q
=
sin (
θ
)
ˆ
X
, where
θ
is the tilt angle between
ˆ
n
and
ˆ
k
, in accordance to Eq. (
1
).
Note that in accordance with its definition, the chirality vector
is zero when the layer normal is parallel to the director. By
making use of Eq. (
1
), we find that the layer normal at the
current configuration
ˆ
k
=
ˆ
Z
.
With the above definitions, the normalized energy-density
function in Eq. (
9
) can be written as
1
w
(
λ
y
,λ
z
,λ
yz
,θ
)
=
1
λ
y
λ
z
+
r
+
1
2
r
λ
2
y
+
r
λ
2
yz
+
λ
2
z
+
r
−
1
2
r
λ
2
y
+
r
λ
2
yz
−
λ
2
z
cos(2
θ
)
FIG. 1. A schematic of the electromechanical response of a
smectic-
A
elastomer. The chirality vector
q
points out of the plane.
−
2
rλ
yz
λ
z
sin(2
θ
)
+
b
1
cos
θ
−
1
2
+
a
sin
2
θ
−
2
cV λ
y
λ
z
sin
θ
−
V
V
d
2
λ
y
λ
z
−
2
sλ
yz
,
(11)
where
c
=
C
μt
0
and
V
2
d
=
μt
2
0
(1
+
χ
)
ε
0
. We also normalize the
polarization [Eq. (
7
)]:
p
=
1
μt
0
P
=
1
λ
E
V
V
2
p
ˆ
X
+
c
q
,
(12)
where
V
2
p
=
μt
2
0
χε
0
. The constants
V
d
and
V
p
have units of voltage
and are related via
1
V
2
d
−
1
V
2
p
=
ε
0
μt
2
0
.
In order to demonstrate the merit of the model, we adopt
the typical values
a
=
1,
b
=
5,
c
=
10
−
4
V
−
1
, and
r
=
1
.
5.
The constants
V
d
=
13
.
7 kV and
V
p
=
15 kV correspond to a
film with a thickness
t
0
=
0
.
1 mm and a susceptibility
χ
=
5.
A. Voltage induced deformation
Due to the electroclinic effect, smectic-
A
elastomers
deform in response to an applied electric field [
1
,
3
,
4
,
6
,
7
].
The aim of this section is to examine the voltage induced
deformations according to the proposed model. To this end,
we take
s
=
0 and minimize the energy-density function
w
(
λ
y
,λ
z
,λ
yz
,θ
)[Eq.(
11
)] with respect to
λ
y
,
λ
z
,
λ
yz
, and
θ
to
obtain the equilibrium state. This leads to nonlinear equations
which can be solved numerically. We can also gain insight by
simplifying the expressions for small voltages and retaining
only the leading order terms in
V
. We obtain
θ
≈
ξV,
(13)
φ
≈
r
−
1
r
ξV,
(14)
λ
y
≈
1
+
(
1
2
(
3
+
b
)
b
+
2
2
1
V
2
d
+
ξ
c
−
1
r
ξ
−
2
r
ξ
2
+
1
2
ξ
2
)
V
2
,
(15)
λ
z
≈
1
+
1
2
(
3
+
b
)
1
V
2
d
+
ξ
c
+
3
r
ξ
−
1
2
ξ
2
V
2
,
(16)
p
·
ˆ
X
=
1
V
2
p
+
cξ
V,
(17)
where
ξ
=
c
a
is the initial slope between the tilt angle and
the applied voltage. This quantity is related to the electroclinic
coefficient
α
=
∂θ
∂E
. Specifically, in the examined parallel plate
geometry
ξ
≈
α
t
0
. Subjected to relatively low voltages, Köhler
et al
.[
21
], Spillmann, Ratna, and Naciri [
4
], and Hiraoka
et al
.
[
6
] measured an electroclinic coefficient
α
=
0
.
045 m
/
MV,
α
=
0
.
006 m
/
MV, and
α
=
0
.
07 m
/
MV, respectively. These
findings correspond to tilt angles in the range 4–10
◦
. Since
the thickness of the smectic-
A
elastomers is around
t
0
∼
032701-3
NOY COHEN AND KAUSHIK BHATTACHARYA
PHYSICAL REVIEW E
96
, 032701 (2017)
FIG. 2. (a) The shear angle
φ
and (b) the tilt angle
θ
as a
function of the applied voltage
V
. The inset plots depict the predicted
response subjected to high voltages. The continuous and the dashed
curves correspond to the exact and the approximated predictions,
respectively.
0
.
01–1 mm, we approximate the order of magnitude of the
initial slope
ξ
∼
0
.
01–1 kV
−
1
.
The parameters
b
and
c
can be determined from a curve fit
of
λ
y
and
λ
z
as a function of the voltage via Eqs. (
15
) and (
16
)
in the low voltage regime. The quadratic dependency of
λ
y
and
λ
z
on the applied voltage stems from the initial microstructural
alignment of the director and the layer normal in the smectic-
A
elastomers.
Figures
2(a)
and
2(b)
depict the shear angle
φ
and the tilt
angle
θ
as a function of the applied voltage. The insets plot
the electromechanical response of the smectic-
A
elastomers
under high voltages. The continuous and the dashed curves
correspond to the exact predictions and the approximations
[Eqs. (
13
) and (
14
)], respectively. A close inspection of the
insets shows that the approximations hold for voltages up to
V
≈
5 kV. It is noted that the application of the electric field
leads to a tilt angle that is larger than the shear angle.
To explain this, note that Eqs. (
13
) and (
14
) reveal that in the
range of small voltages the ratio between the shear and the tilt
angles is
φ
θ
=
r
−
1
r
. Interestingly, in smectic-
A
elastomers with
r>
1 the ratio
φ
θ
<
1. This implies that the electrically induced
deformation is characterized by a tilt angle
θ
that is greater
than the shear angle
φ
. This observation is supported by the
experimental findings of Hiraoka
et al
.[
6
], where the measured
FIG. 3.
p
·
ˆ
X
as a function of the applied voltage
V
. The inset
depicts the predicted response as a result of high voltages. The
continuous and the dashed curves correspond to the polarization
[Eq. (
12
)] with
c
=
10
−
4
and 0 V
−
1
, respectively.
ratio
φ
θ
≈
0
.
02 led to the anisotropy parameter
r
=
1
.
02. It
should also be noted that the obtained value of the anisotropy
parameter is very low.
The predicted polarization is plotted as a function of the
voltage in Fig.
3
. As before, the inset plots the predicted
polarization as a result of high voltages. The continuous and
the dashed curves correspond to the polarization with
c
=
10
−
4
and 0 V
−
1
, respectively. It is shown that the electrically
induced tilt leads to an enhancement of the polarization of
the elastomer by as much as 1.8 times. We also find that
reversing the electric field flips the direction of the polarization,
as seen in the experimental findings of Hiraoka
et al
.[
6
]. As
previously discussed, this effect emanates from the chirality
of the mesogens. We emphasize that since the direction of the
chirality parameter
q
is linearly coupled to the polarization, the
spontaneous dipoles in a smectic-
A
elastomer with
q
→−
q
decrease the polarization.
B. Deformation induced polarization
Next, we follow the experiment of Kramer and Finkelmann
[
16
] and the theoretical work of Stenull
et al
.[
13
] and apply
a shear stress
s
along the
ˆ
Y
direction which deforms the
elastomer such that
λ
y
=
λ
z
=
1, corresponding to simple
shear. By setting
V
=
0, the normalized energy density is
w
=
w
(
λ
y
=
1
,λ
z
=
1
,λ
yz
,θ
)
.
(18)
To determine the equilibrium state, we minimize Eq. (
18
) with
respect to the shear
λ
yz
and the tilt angle
θ
.
Figures
4(a)
and
4(b)
depict the shear angle
φ
and the tilt
angle
θ
as a function of
s
, respectively. We emphasize that the
nonlinear dependence of the shear angle
φ
on the shear force
s
corresponds to a linear dependence between the shear
λ
yz
and
s
. It is demonstrated that the mesogens tilt in the direction
of the external shear force, as expected, and negating the force
simply flips the direction of the tilt angle
θ
.
Interestingly, as opposed to the previous case the proposed
model predicts that under the examined deformation the shear
angle is always greater than the tilt angle. Furthermore,
since
λ
z
=
1 is fixed and layer rotation is not permitted, the
tilt of the mesogens is not accompanied by a reduction in
032701-4
ELECTROCLINIC EFFECT IN CHIRAL SMECTIC-
A
. . .
PHYSICAL REVIEW E
96
, 032701 (2017)
FIG. 4. (a) The shear angle
φ
and (b) the tilt angle
θ
as a function
of the normalized nominal stress
s
.
the layer thickness. These observations are in accord with the
experimental work of Kramer and Finkelmann [
16
], where the
elastomer was subjected to a pure shear deformation such that
the maximum imposed shear angle was
φ
=
21
◦
.
An examination of Fig.
4(b)
reveals that the dependence
between the tilt of the mesogens and the applied shear force
can be divided into two regimes. Initially, a small increase
in the applied force leads to a large variation in the tilt angle.
Beyond a threshold
s
≈
5, the cost of further tilting the director
becomes very high. Specifically, a large increase in the shear
force is required to attain a small increment in the tilt of the
mesogens. Practically, elastomers fail at a shear angle
φ
90
◦
and may not reach the second regime.
From Eq. (
7
) we deduce that the polarization in the absence
of an electric field is directed along the chirality parameter
q
,
which is perpendicular to the plane spanned by
ˆ
n
and
ˆ
k
.
Figure
5
plots the normalized polarization component
p
·
ˆ
X
as a function of
s.
It is demonstrated that the direction
of the applied shear force determines the direction of the
polarization. Two reasons lead to this behavior—the first is
the antisymmetric dependency of the tilt angle on the applied
force. As discussed in relation to Fig.
4(b)
, the directions of
the shear force and the tilt are linearly dependent. Second, the
observed trend is enforced by the direct coupling between the
director and the chirality parameter
q
. This coupling is defined
such that the magnitude of the polarization is proportional
to sin (
θ
). We emphasize that the chirality of the mesogens
determines the direction of the polarization. Specifically, a
FIG. 5.
p
·
ˆ
X
as a function of the normalized nominal stress
s
.
shear force along the
ˆ
Y
direction induces a polarization along
the
−
ˆ
X
direction in elastomers characterized by a flipped
chirality vector (i.e.,
q
→−
q
).
IV. THE ELECTRICALLY INDUCED DEFORMATION
OF A CYLINDER WITH HELICAL LAYERS
To illustrate the capabilities of the electroclinic effect,
we explore the electromechanical response of a cylindrical
smectic-
A
elastomer with helical layers. Consider a hollow
thin-walled cylinder with an initial radius
R
, a wall thickness
t
0
, and a height
H
. The cylinder is made of a thin incompress-
ible smectic-
A
elastomer film with helical layers. We prescribe
a global cylindrical coordinate system
{
ˆ
R
,
ˆ
,
ˆ
Z
}
such that at
the reference configuration the director and the layer normal
are
ˆ
n
0
=
ˆ
k
0
=
sin(
ψ
0
)
ˆ
+
cos(
ψ
0
)
ˆ
Z
,
(19)
where
ψ
0
is the initial helix angle (see Fig.
6
). The positions
of the referential material points are denoted by
{
R, ,Z
}
.The
inner and outer faces of the cylindrical film are covered with
flexible electrodes such that an applied voltage induces a radial
electric field
E
≈
V
t
0
ˆ
R
. In addition, an internal pressure
p
is
exerted on the inner surface. As a result, the cylinder deforms
,
FIG. 6. A schematic of the electromechanical response of a
cylinder with helical layers.
032701-5
NOY COHEN AND KAUSHIK BHATTACHARYA
PHYSICAL REVIEW E
96
, 032701 (2017)
FIG. 7. (a) The twist
δ
and (b) the change in the enclosed volume
v
t
V
t
as a function of the applied voltage for various helix angles at
p
=
0.
such that at the current deformation
r
=
λR
2
+
β, θ
=
+
δZ, z
=
1
λ
Z,
(20)
where
{
r,θ ,z
}
denote the material points at the current config-
uration. Accordingly, the deformation gradient in cylindrical
coordinates is
F
=
⎛
⎜
⎝
∂r
∂R
00
0
r
R
∂θ
∂
r
∂θ
∂Z
00
∂z
∂Z
⎞
⎟
⎠
=
⎛
⎜
⎝
λ
R
r
00
0
r
R
δr
00
1
λ
⎞
⎟
⎠
.
(21)
At the current configuration, the chirality vector is
q
=
sin
(
φ
)
ˆ
R
,
(22)
and the director is
ˆ
n
=
sin(
ψ
+
φ
)
ˆ
+
cos(
ψ
+
φ
)
ˆ
Z
,
(23)
where
φ
is the angle between the director
ˆ
n
and the layer
normal
ˆ
k
and
ψ
is the helix angle at the current configuration.
By employing Eq. (
1
), we find that
ˆ
k
=
sin(
ψ
)
ˆ
+
cos(
ψ
)
ˆ
Z
.
(24)
The layer normal at the current configuration is directly
dependent on the deformation gradient via [
10
,
12
]
ˆ
k
=
F
−
T
ˆ
k
0
|
F
−
T
ˆ
k
0
|
,
(25)
and therefore one finds that the helix angle at the current
configuration is
ψ
=
arctan
R
sin
(
ψ
0
)
λr
[
cos
(
ψ
0
)
−
Rδ
sin
(
ψ
0
)
]
.
(26)
The energy-density function can be written as
w
(
h
)
(
λ,β,φ,δ
)
=
w
(
λ,β,φ,δ
)
−
2
μ
p
r
R
2
1
λ
,
(27)
where
w
(
λ,β,φ,δ
) is computed with Eq. (
9
),
p
is the pressure
inside the cylinder, and (
r
R
)
2
1
λ
accounts for the change in
the enclosed volume of the cylinder. The electromechanical
response of the cylinder is determined by minimizing the
energy
w
(
h
)
(
λ,β,φ,δ
) with respect to
λ
,
β
,
φ
, and
δ
.
Figures
7(a)
and
7(b)
depict the twist
δ
and the change in the
enclosed volume
v
t
V
t
=
πr
2
h
πR
2
H
, where
h
is the current height of
the cylinder, as a function of the applied voltage, respectively,
for various initial helix angles at
p
=
0. It is shown that
elastomers with cylindrical layers (or helical layers with an
angle
ψ
0
=
0
◦
) experience twist with a minor increase in the
enclosed volume. As the helix angle increases from
ψ
0
=
0to
45
◦
, the twist
δ
decreases while the enclosed volume increases.
Interestingly, the twist angle in the smectic-
A
elastomers with a
helix angle
ψ
0
=
45
◦
is extremely small. The proposed model
also predicts that such an elastomer experiences the largest
volumetric expansion. The trend is reversed as the helix angle
increases from
ψ
0
=
45 to 90
◦
. Specifically, the cylinder twists
in the other direction and the enclosed volume decreases.
We highlight that the volumetric expansion of the smectic-
A
elastomers is due to the radial expansion and the longitudinal
compression resulting from the amplification of the voltage.
Additionally, it is important to note that, due to the low voltages
that are practically applied in the actuation of these materials,
the twist and the volumetric expansions are small.
Next, we study the effects of pressure on the actuation of
smectic-
A
elastomers with helical layers at an angle
ψ
0
=
45
◦
.
FIG. 8. (a) The twist
δ
and (b) the change in the enclosed volume
v
t
V
t
as a function of the applied voltage for a cylinder comprising helical
layers with a
ψ
0
=
45
◦
angle for various pressures.
032701-6
ELECTROCLINIC EFFECT IN CHIRAL SMECTIC-
A
. . .
PHYSICAL REVIEW E
96
, 032701 (2017)
Figures
8(a)
and
8(b)
plot the twist
δ
and the change in the
enclosed volume
v
t
V
t
as a function of the applied voltage,
respectively, for a cylinder subjected to various pressures.
Figure
8(a)
reveals that the pressure causes the elastomer to
twist. However, the increase in voltage does not significantly
change the twist angle.
As shown in Fig.
8(b)
, the enclosed volume increases
with pressure. This stems from the radial expansion and the
longitudinal elongation of the cylinder in response to the
internal pressure. As the voltage is amplified, a linear increase
in the enclosed volume is observed. Interestingly, we find
that the linear dependence between
v
t
V
t
and
V
is independent
of the pressure in the examined range of voltages. It is
again noted that, even under an applied pressure, an increase
in voltage leads to a radial expansion and a longitudinal
compression.
V. CONCLUDING REMARKS
This paper proposes a new model that captures the behavior
of smectic-
A
elastomers subjected to an electromechanical
loading. Specifically, the energetic formulation of Adams
et al
.
[
12
] is rewritten in terms of a new variable
q
, that describes
the chirality of the mesogens composing the elastomer, and
an additional term is added to account for the effects of an
electric stimulus. The current configuration of the elastomer is
determined by minimizing the energy.
To demonstrate the merit of the model, two simple cases
are presented. In the first case, an electric field is applied
to the smectic-
A
elastomer. In accordance with experimental
findings, the proposed model predicts a shear deformation of
the sample. In the second case, a mechanical shear deformation
is imposed. As a consequence, the mesogens tilt and their
inherent chirality gives rise to a polarization. The predictions
qualitatively capture several experimental observations.
The potential of these elastomers is exhibited in Sec.
IV
,
where a thin-wall cylindrical configuration of a smectic-
A
film with helical layers is considered. Upon application of an
electric field, we show that the pitch and the helix angle control
the macroscopic behavior. Specifically, the enclosed volume
and the twist depend on the formation of the lamellar structure.
It is also shown that the application of internal pressure in the
tube can effectively change the macroscopic response.
ACKNOWLEDGMENT
We are grateful for the financial support of the US Air
Force Office of Scientific Research through the MURI Grant
No. FA9550-16-1-0566.
[1] S. Garoff and R. B. Meyer, Electroclinic Effect at the
A
−
C
Phase Change in a Chiral Smectic Liquid Crystal,
Phys. Rev.
Lett.
38
,
848
(
1977
).
[2] S. Garoff and R. B. Meyer, Electroclinic effect at the
A
−
C
phase change in a chiral smectic liquid crystal,
Phys. Rev. A
19
,
338
(
1979
).
[3] W. Lehmann, H. Skupin, C. Tolksdorf, E. Gebhard, R. Zentel, P.
Krüger, M. Lösche, and F. Kremer, Giant lateral electrostriction
in ferroelectric liquid-crystalline elastomers,
Nature (London)
410
,
447
(
2001
).
[4] C. M. Spillmann, B. R. Ratna, and J. Naciri, Anisotropic
actuation in electroclinic liquid crystal elastomers,
Appl. Phys.
Lett.
90
,
021911
(
2007
).
[5] K. Hiraoka, M. Kobayasi, R. Kazama, and H. Finkelmann,
Electromechanics of monodomain chiral smectic
c
elastomer:
Mechanical response to electric stimulation,
Macromolecules
42
,
5600
(
2009
).
[6] K. Hiraoka, T. Kishimoto, M. Kato, and T. Tashiro, Electroclinic
and electromechanical effects of a side-chain chiral smectic
a
elastomer,
Liq. Cryst.
38
,
489
(
2011
).
[7] C. M. Spillmann, J. Naciri, B. R. Ratna, R. L. B. Selinger,
and J. V. Selinger, Electrically induced twist in smectic liquid-
crystalline elastomers,
J. Phys. Chem. B
120
,
6368
(
2016
).
[8] E. Nishikawa, H. Finkelmann, and H. R. Brand, Smectic
A
liquid
single crystal elastomers showing macroscopic in-plane fluidity,
Macromol. Rapid Commun.
18
,
65
(
1997
).
[9] E. Nishikawa and H. Finkelmann, Smectic-
A
liquid single
crystal elastomers–strain induced break-down of smectic layers,
Macromol. Chem. Phys.
200
,
312
(
1999
).
[10] J. M. Adams and M. Warner, Elasticity of smectic-
A
elastomers,
Phys.Rev.E
71
,
021708
(
2005
).
[11] O. Stenull and T. C. Lubensky, Unconventional elasticity in
smectic-
A
elastomers,
Phys. Rev. E
76
,
011706
(
2007
).
[12] J. M. Adams, M. Warner, O. Stenull, and T. C. Lubensky,
Smectic-
A
elastomers with weak director anchoring,
Phys. Rev.
E
78
,
011703
(
2008
).
[13] O. Stenull, T. C. Lubensky, J. M. Adams, and M. Warner,
Smectic-
C
tilt under shear in smectic-
A
elastomers,
Phys. Rev.
E
78
,
021705
(
2008
).
[14] D. Corbett and M. Warner, Electromechanical elongation of
nematic elastomers for actuation,
Sens. Actuat. A
149
,
120
(
2009
).
[15] J. V. Selinger, J. Xu, R. L. B. Selinger, B. R. Ratna, and R.
Shashidhar, Theory of chiral modulations and fluctuations in
smectic-
A
liquid crystals under an electric field,
Phys.Rev.E
62
,
666
(
2000
).
[16] D. Kramer and H. Finkelmann, Shear-induced tilt in smectic-
A
elastomers,
Phys.Rev.E
78
,
021704
(
2008
).
[17] M. Warner and E. M. Terentjev,
Liquid Crystal Elastomers
(Oxford University Press, Oxford, 2003), Vol. 120.
[18] Y. C. Shu and K. Bhattacharya, Domain patterns and macro-
scopic behaviour of ferroelectric materials,
Philos. Mag. Part B
81
,
2021
(
2001
).
[19] M. Brehmer, R. Zentel, F. Giesselmann, R. Germer, and P.
Zugenmaier, Coupling of liquid crystalline and polymer network
properties in lc-elastomers,
Liq. Cryst.
21
,
589
(
1996
).
[20] P. Archer and I. Dierking, Experimental determination of the full
landau potential of bent-core doped ferroelectric liquid crystals,
Phys. Rev. E
72
,
041713
(
2005
).
[21] R. Köhler, R. Stannarius, C. Tolksdorf, and R. Zentel,
Electroclinic effect in free-standing smectic elastomer films,
Appl. Phys. A
80
,
381
(
2005
).
032701-7