Scholarly article on topic 'Control of colloidal placement by modulated molecular orientation in nematic cells'

Control of colloidal placement by modulated molecular orientation in nematic cells Academic research paper on "Materials engineering"

Share paper
Academic journal
Science Advances
OECD Field of science

Academic research paper on topic "Control of colloidal placement by modulated molecular orientation in nematic cells"


Control of colloidal placement by modulated molecular orientation in nematic cells

2016 © The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC). 10.1126/sciadv.1600932

Chenhui Peng, Taras Turiv, Yubing Guo, Sergij V. Shiyanovskii, Qi-Huo Wei, Oleg D. Lavrentovich*

Colloids self-assemble into various organized superstructures determined by particle interactions. There is tremendous progress in both the scientific understanding and the applications of self-assemblies of single-type identical particles. Forming superstructures in which the colloidal particles occupy predesigned sites and remain in these sites despite thermal fluctuations represents a major challenge of the current state of the art. We propose a versatile approach to directing placement of colloids using nematic liquid crystals with spatially varying molecular orientation preimposed by substrate photoalignment. Colloidal particles in a nematic environment are subject to the long-range elastic forces originating in the orientational order of the nematic. Gradients of the orientational order create an elastic energy landscape that drives the colloids into locations with preferred type of deformations. As an example, we demonstrate that colloidal spheres with perpendicular surface anchoring are driven into the regions of maximum splay, whereas spheres with tangential surface anchoring settle into the regions of bend. Elastic forces responsible for preferential placement are measured by exploring overdamped dynamics of the colloids. Control of colloidal self-assembly through patterned molecular orientation opens new opportunities for designing materials and devices in which particles should be placed in predesigned locations.

Q. CD Q.

f —i


Self-assembled colloidal structures are formed by small particles that interact in a predefined manner to create complex patterns (1). Self-assembly in isotropic solutions such as water already demonstrated a great potential to design materials for various scientific and technological purposes (2). Colloidal interactions and thus complexity of superstructures can be markedly enriched by using liquid crystals (LCs) as a dispersion medium (3). Long-range orientational order of LCs leads to anisotropic colloidal interactions that control assembly into clusters (4-13); mediate interactions of particles with flat (14), curved (15), and templated interfaces (16, 17); and cause trapping at topological defects (18-23). The spectrum of possible design pathways of colloidal self-assembly through LC dispersions is further expanded by sensitivity of the liquid crystalline order to external electromagnetic fields, temperature, or confinement.

Colloidal self-assembly in both isotropic and anisotropic media is currently well explored for the case when the particles are of the same type and occupy locations determined by interactions with other particles and bounding surfaces. The challenge is to develop robust approaches to placing particles into predesigned locations. Current methods of heterogeneous design are based on lock-and-key colloidal interactions of nonspherical particles (24), programmable interactions of DNA molecules attached to the surface of particles (25), and use of substrates with patterned magnetic (26) or hydrophobic properties (27). Here, we propose a new approach to the placement of colloidal particles that is based on LCs with predesigned molecular orientation, achieved through the substrate photopatterning.

We use the simplest type of LC, the so-called nematic, in which elongated molecules are oriented along a single direction called the director n, with the properties n = —n and n2 = 1. The flat nematic cells bounded by two glass plates are designed in such a way that the direc-

Liquid Crystal Institute and Chemical Physics Interdisciplinary Program, Kent State University, Kent, OH 44242, USA.

Corresponding author. Email:

tor field experiences periodic one-dimensional (1D) distortions with alternating stripes of splay and bend. These deformations are preim-posed by the substrate photoalignment technique developed recently (28, 29). Spherical particles demonstrate a marked selectivity in spatial locations within the preimposed landscape of director deformations, depending on their surface-anchoring properties, that is, the preferred direction of the local n. Namely, spheres with tangential surface alignment of n are attracted to the regions of bend, whereas spheres with perpendicular alignment are attracted to the regions of splay. Analysis of the colloidal dynamics allows us to determine the interaction potential of the particles with the underlying director landscape.


We explore a stripe pattern with periodic splay divn ^ Oandbendn x curln ^ 0 modulated along the y axis. The director is everywhere, parallel to the bounding plates; therefore, there is no z component, nz = 0. The director field imposed by the surface photoalignment layer is of the form

n = (nxnynz ) = (cos a, sin a, 0)

where a = py/l and l = 80 mm is the period (Fig. 1). An experimental cell represents a thin (with a thickness of h = 20 mm) flat layer of a nematic sandwiched between two glass plates with identical photo-induced alignment patterns at the top (z = h) and bottom (z = 0) surfaces. Splay and bend contribute to the Oseen-Frank free-energy density of the distorted nematic

f = 1 Ki(divn)2+1K2(n■ curln)2+1K3(n x curln)2 (2)

O t o cr

CD —i

Fig. 1. Elasticity-directed placement and assembly of normally anchored polystyrene spheres in regions with splay deformations of the patterned director. (A) A sphere with perpendicular surface anchoring placed in a uniform nematic causes an appearance of a hyperbolic hedgehog that can be located on either the left- or right-hand side; the resulting director deformations are of a dipolar type. (B) Periodic splay-bend stripe pattern; the normalized splay energy density is labeled by colors according to the scale on the right-hand side. A sphere is placed in the bend region by optical tweezers and released there, and migrates toward the splay region. The dark blue curve is the typical experimental trajectory of the sphere's center. (C) Self-assembly of spheres into linear chains in the regions of maximum splay with p >0 [bright-field microscope with unpolarized light, and nematic liquid crystal pentylcyanobiphenyl (5CB)]. (D) Polarizing microscopy images of a sphere moving from one splay region to the next splay region, with concomitant reorientation of the structural dipole from p < 0 to p >0; the white arrows point toward the hyperbolic hedgehog. (E) Typical experimental trajectory of the sphere within the director configuration mapped in the PolScope mode of observation of the experimental cell. (F) Experimentally measured x(t) dependence for a sphere moving from y = l to y = 0 and its theoretical fit by Eq. 5. (G) Experimentally measured y(t) dependence for a sphere moving from y = l to y = 0 and its theoretical fit by Eq. 6. All spheres have a radius of R = 2.5 mm; the period of the patterned director is l =80 mm. The nematic material corresponding to parts (D) to (G) is MLC6815.

Q. CD Q.

f —i

O t o cr

CD —i

with the elastic constants of splay (Kj), twist (K2), and bend (K3). For the director patterns (Eq. J) used in this work, twist (as well as saddle-splay) deformations vanish because the director does not change along the z axis. Figure JB shows the spatial distribution of the dimensionless splay energy density (l/p)2(divn)2 = cos2(py/l) corresponding to Eq. J. As detailed by the experiments below, the alternating regions of splay and bend attract colloidal spheres with different surface anchoring: spheres with normal anchoring are attracted to the splay regions (Fig. J), and spheres with tangential anchoring are attracted to the bend regions, as discussed later.

Spheres with normal anchoring

Polystyrene spheres with a radius of R = 2.5 mm were covered with a thin layer of a surfactant to produce perpendicular (also called homeotropic or normal) alignment of the director at their surfaces. When such a sphere is placed in a nematic with a uniform director field n0 = (J, 0,0), the radial-like director around the sphere is to-pologically compensated by a point defect, the so-called hyperbolic hedgehog located on either the right- or left-hand side of the sphere (3, 30). The resulting director structure is in the form of a dipole (Fig. JA). We direct the unit vector p of the dipole from the hedgehog toward the sphere. Note that once the direction of the dipole is randomly selected to be either parallel or antiparallel to the x axis of the uniform cell with n0 = (J, 0,0), it cannot be reversed into the opposite orientation because the energy barrier for such a realignment, on the order of KR » (J0-J6 to J0-J7) J, is much larger than the thermal energy kBT = 4.2 x J0-2J J, with the typical value of the Frank elastic constant K being a few tens of piconewton. In principle, in very shallow cells (with a thickness comparable to the colloid diameter), the dipolar structure can transform into a Saturn ring configuration (31); however, in all our experiments, the dipolar structure remained stable (32, 33). This implies that the homeotropic anchoring is strong, with the anchoring coefficient being larger than K/R » J0-5 J/m2.

Individual spheres with normal anchoring placed in the cell with a periodic splay and bend migrate toward the regions with the maximum splay, y = 0, ± l, ± 2l, ... (Fig. JB). When the concentration of spheres is high enough, they form linear chains along the x axis, located at the same ordinate of the maximum splay, y = 0, ± l, ± 2l,... (Fig. JC and fig. SJ).

It is important to note a polar character of colloidal assemblies in the splay regions. In a uniform nematic, the "left" and "right" orientations of p are equally probable. In the described splay-bend pattern, which lacks the left-right mirror symmetry, all the dipoles align themselves along the same direction, forming a ferroelectric type of order. For the chosen geometry of the splay-bend pattern (Eq. J and Fig. JB), the dipoles p of the spheres trapped at the splay locations y = 0, ± l, ± 2l, ... are always parallel to the x axis. If we use optical tweezers and place a sphere in any location within the pattern specified by Eq. J, then the sphere orients its dipole pp parallel or antiparallel to the local director. If the sphere is released from the optical tweezers, it moves toward the regions y = 0, ± l, ± 2l, ... while simultaneously rotating the dipole in accordance with the rotating local director orientation (Fig. JB). In particular, the tweezers can place a sphere in the splay region in such a way that the dipole pp is antiparallel to the x axis. This configuration is unstable because pp prefers to be parallel to the x axis at equilibrium. The sphere with the "wrong" orientation of pp moves across the stripe pattern, rotating pp in accordance with the

local director, until it reaches the neighboring splay region with p being parallel to the x axis, that is, completing a 180° rotation through a shift by |Ay| = l (Fig. 1, D and E, and movie S1). Note that the simple reorientation of p by 180° with the sphere keeping its y coordinate unchanged (|Ay| = 0) is never observed in our experiments; the reason is that this process is associated with a huge energy barrier KR >> kBT, similar to the case of the uniform cell discussed above.

At moderate concentrations, colloidal particles with homeotropic anchoring gather into linear chains (Fig. 1C). Note that all the chains have the same direction of the individual dipole moments p and are separated by the same distance l along the y axis. By increasing the concentration of colloids, the length of the chains can be increased, but because we are primarily interested in the individual behavior, we limit the study by using relatively low concentrations, in which the number of particles in each chain is less than about 10.

The interaction potential between a colloidal sphere with a dipolar director and a distorted director has been derived by Perga-menshchik (34). We follow the study of Pergamenshchik (34) but use a modified expression that conforms to our definition of the dipole; for a sphere of radius R << l, the interaction potential writes

Ud = -8rcbKR2(p ■ n) divn

where b is the amplitude of the dipole director distortions that depends on the strength of surface anchoring, anisotropy of elastic constants, accuracy of the dipolar approximation, etc., and K is the effective elastic constant in the so-called one-constant approximation. In the analysis below, we treat b as an unknown numerical coefficient, which we determine from the fitting of experimental data. In Eq. 3, we introduce the scalar product (p ■ n) to account for the polar character of coupling between the splay divn and p and also to preserve the nematic symmetry n = — n. For the director field in Eq. 1, n = [ cos(py/l), sin(py/l), 0], the interaction potential for a given sphere that keeps (p ■ n) = 1 [namely, Ud = -8p2bK(R2/l)cos(py/l)] reaches its minimum Ud>min = -8p2bK(R2/l) at the locations y = 0, ± 2l, ... and maximum U<j;max = 8p2bK(R2/l) at the locationsy = ± l, ± 3l,.... That is, the normally anchored spheres are expected to gather in the regions with a maximum splay, being oriented in such a way that the dipole p fits the polarity of splay n divn.

The potential (Eq. 3) also explains the behavior of a sphere placed using the optical tweezers in a splay region with the wrong polarity. Such a sphere moves along the potential energy landscape from the maximum energy at y = l toward the minimum at y = 0 (Fig. 1, D and E). The dynamics of the sphere from y = l toward y = 0 can be used to extract the information about the energetics of elastic coupling between the colloids and the surrounding director distortions, as detailed below.

The particle's dipole preserves its orientation parallel to the local director during migration. Then, the elastic force driving the sphere to the minimum energy location writes

F = -(3Ud/3y)êy = -8p3Pk(R)2 sin(py

where êy is the unit vector along the y axis. Because the Reynolds number in our experiments is small and because inertia contributions can be neglected, the driving force is balanced by the viscous

Oí Q. (D Q.

f —i

o t o cr

(D —i

drag force Fd = —6nRrfv, where rf is the diagonal viscosity tensor with the components fy and r|±, and v is the particle velocity vector with the components vy and v±; the subscripts || and ± refer to the components parallel and perpendicular to the director, respectively. The drag force can be rewritten in terms of the components of the viscosity tensor, Fd/6pR = fyvy(ex cos a + ey sin a) + r±v±(-ex sin a + ey cos a). By balancing the driving and drag forces in the laboratory coordinate system, we obtain the equations for the velocity com-p- — v± sin y and vy = vy sin (y) + v± cos(y) in

ponents vx = vy cos y — vx sin—

the following form dx

4p2 hx — h|| 3 h||hx

„ KR 2 py py ß ^sin7 cos7

KR hx sin2 py + hy

,c2 py

hx h||

8p4ß f K is about 1.9 x 10-

Solutions of Eqs. 5 and 6 establish how the sphere's coordinates, x and y, vary with time. These theoretical dependencies are used to fit the experimental trajectory in Fig. 1 (F and G) to verify the model and to deduce the unknown parameter p.

Fitting of the experimental dynamics of the colloids requires knowledge of the two components of the viscosity tensor, fy and

These two were experimentally determined for the same LC and normally anchored spheres by measuring the mean square displacements (MSDs) as a function of time for the Brownian diffusion parallel and perpendicular to the director in uniformly aligned cells, n0 = (1,0,0) (Fig. 2). The slopes of the MSDs versus time dependencies yield the diffusion coefficient Dy = (7.6 ± 0.2) x 10-15 m2/s for the direction parallel to n0 and D± = (4.1 ± 0.2) x 10-15 m2/s for the perpendicular direction. These values yield the corresponding components of the viscosity tensor, r y = (12 ± 1) mPa ■ s and = (22 ± 1) mPa ■ s, which follow from the Stokes-Einstein formula r = kBT/(6nRD).

As clearly seen from Fig. 1 (F and G), the experimental data are closely matched by the expected theoretical behavior, despite the fact that the model uses the so-called one-constant approximation. The only fitting parameter is p = 0.20 ± 0.02. Note that with this value and with K = 10-11 N, the depth of the potential minimum [Ud,max - Ud,min = 16n2pK(R2/l)] is 2.2 x 10-17 J, that is, about four orders of magnitude higher than the thermal energy kBT = 4.2 x 10-21 J. In other words, the elastic trapping is much stronger than the randomizing forces associated with thermal fluctuations. One can also estimate the stiffness of the potential trap that keeps the normally anchored spheres in the splay regions of a "correct" polarity. For y ^ 0, the stiffness of the elastic trap k = ^g^r =

Fig. 2. MSD of a dipolar colloid in uniformly aligned LC cell. Homeo-tropically anchored sphere with a radius of 2.5 mm in a nematic MLC6815 cell with uniform alignment n0, which is achieved by unidirectional rubbing of the aligning layer of PI2555 deposited on the bounding plates.

pattern. The tangentially anchored spheres migrate toward the regions with the maximum bend, that is, locations at y =± 1/2, ± 3//2,... (Fig. 3, B and C). As the concentration of spheres increases, they assemble into chains with tilted arrangements of the neighboring spheres, similar to those described for uniform nematic cells in the work of Smalyukh et al. (5). Figure 3C shows only the short chains since we are primarily interested in the individual behavior of colloids in the patterned landscape. Longer chains form in more concentrated dispersions and can extend over a length of 50 to 70 mm across the region with bend distortions; they are disrupted when their ends become close to the regions with splay deformations.

The potential of interaction between a quadrupolar particle and the distorted director field reads, following Pergamenshchik (34), as

UQ = —4pKcR3(n -V)divn

Spheres with tangential anchoring

Polystyrene spheres without surfactant coatings produce spontaneous tangential (parallel to the interface and degenerate) anchoring. The director deformations around these spheres in a uniform nematic are of quadrupolar symmetry, with two point defects, called boojums, at the poles (Fig. 3A). Placing these spheres in the splay-bend stripe pattern (Eq. 1) produces behavior qualitatively different from the normally anchored spheres (Fig. 3). Figure 3B maps the dimension-less bend energy (l/n)2(n x curln)2 = sin2(rcy/l) stored in the stripe

where % is the dimensionless coefficient of the quadrupolar moment, the value of which is unknown and depends on the surface anchoring. For the patterned director field (Eq. 1), this potential reduces to Uq = 4rc3%KR3 sin2 py, thus predicting that the spheres with % <0 will minimize the potential by migrating toward the locations with bend (y = ± l/2, ± 3l/2, ...), as observed in the experiments (Fig. 3C). The corresponding force driving the spheres to the equilibrium locations is

Fq = — (dUq=dy) = — 4p4cK( y ) sin2^

If the sphere is placed in the splay region y = 0, ± l,..., using the optical tweezers, a fluctuative displacement from that location causes the sphere to move toward the bend stripe region at y = ± l/2, ± 3l/2 ..., as shown by a trajectory in Fig. 3B. By balancing the driving and drag

03 Q. CD Q.

f —i

CT CD —i

Fig. 3. Elasticity-directed placement and assembly of tangentially anchored polystyrene spheres in regions with bend deformations of the patterned director. (A) Director deformations around a tangentially anchored sphere in a uniform nematic are of quadrupolar symmetry, with two point defects, called boojums, at the poles. (B) Periodic splay-bend stripe pattern; the normalized bend energy density is labeled by colors according to the scale on the right-hand side. A sphere placed and released in the splay region using optical tweezers migrates toward the bend region. The dark red curve is the typical experimental trajectory of the sphere's center. (C) Self-assembly of spheres into chains in the regions of maximum bend (bright-field microscope and unpolarized light). (D) Experimentally measured x(t) dependence for a sphere moving from y = / to y = //2 and its theoretical fit by Eq. 9. (E) Experimentally measured y(t) of a sphere moving from y = l to y = l/2 and its theoretical fit by Eq. 10. All spheres have radius of R = 2.5 mm; periodicity of the director pattern is / = 80 mm. The LCs used are 5CB (C) and MLC6815 (D and E).

Q. CD Q.

f —i

O t o cr

CD —i

forces, one finds the x and y components of the velocity, similarly to the case of normally anchored spheres

P3 hi -

h y KR2

dy dt:

2p3 KR2 hi. sin2 py + hy cos2 py . 2py

hll hi

Proceeding as above, by solving Eqs. 9 and 10 and by finding the theoretical trajectory of the particle, we fit the experimental dependencies x(t) andy(t) (Fig. 3, D and E). The model describes the experimental behavior well and yields the value of quadrupolar parameter % = -0.95 ± 0.10. With this value, the depth of the trapping potential Uq, max — Uq, mm = 4p31 % | can be estimated as 2.9 x 10-18 J, which is much larger than the thermal energy. Finally, we can estimate the stiffness k = |y^;/2 = 8p5K|%| R^fthepotential, which is 9 x 10-9N/m, that is, about two times weaker than the stiffness of the splay trap acting on the normally anchored spheres.

Fig. 4. Preferential placement of colloids with different surface properties in a cell with periodically patterned splay-bend director field. (A) Polarizing microscope texture of polystyrene colloids with homeotropic and tangential anchoring located in the alternating bands of splay and bend, respectively. (B) Scheme of preferential placement with the experimentally determined director pattern (PolScope observations). The homeotropically anchored spheres reside in the region of splay (y = 0), whereas tangentially anchored spheres localize in the regions of bend (y = ±l/2). Note the different orientations of chains formed by homeotropic spheres (along the x axis) and by the tangential spheres (perpendicular to the x axis). The colloids are placed in an LC mixture with zero dielectric anisotropy. Polarizer is parallel to the x axis and analyzer along the y axis.

Preferential placement of spheres with different surface anchoring

So far, we have discussed dispersions of single-type particles. When the cell is filled with a mixture of two types of particles, the surface-patterned director landscape leads to their preferential placement in different regions of the cell (Fig. 4). Namely, the normally anchored spheres migrate toward the splay (y = 0), and the tangentially anchored spheres migrate toward the bend (y =± l/2) (Fig. 4). Depending on the overall concentration of the spheres, they can form larger or smaller aggregates. The normally anchored spheres produce linear chains parallel to the direction of stripes, whereas the tangentially anchored spheres form chains that are predominantly perpendicular to this direction, which is an expected behavior, as evidenced by previous studies of aggregation of both types of spheres in uniformly aligned nematic cells (Fig. 4) (3).


Here, we demonstrate a new approach to controlling colloidal placement, by using the long-range elastic forces in spatially modulated ne-matic LC. The approach is based on the predesigned molecular orientation pattern at the bounding plates of the cell. The colloidal spheres interact with the underlying director distortions and migrate into the regions where the elastic energy of interactions is minimized. For spheres with normal anchoring, these regions correspond to the splay deformations; in contrast, tangentially anchored spheres are localized in the regions of bend.

The periodic splay-bend director pattern imposes strong limitations on the collective behavior of the colloids. For example, the spheres with normal anchoring form a ferroelectric type of macroscopic ordering, in which the structural dipole moment p of each and every particle points to the same direction. The system is very deterministic since we do not observe any "mistakes" made by the colloids, for example, wrong orientation of p or wrong placement of the sphere with the tangential anchoring (which would be in the regions

of splay). Note here that we explored four different nematic materials, and all of them showed similar behavior, which underscores a general and robust character of the colloidal interactions mediated by surface anchoring and elasticity of the preimposed director patterns.

Because the period l of the splay-bend background is much bigger than the particle radius R, the interaction potential is adequately described by the lowest multipole deformation around the particles that grows as R2/l for the spheres with locally dipolar director configuration around them and as R3/l2 for quadrupolar spheres with tangential anchoring. Thermal fluctuations do not affect the particle dynamics in our experiments. For the given strength of director gradients, defined by the period l = 80 mm, the smallest radius of the spheres that can be trapped at the desired location in space can be estimated by comparing the depth of the trapping potentials to the thermal energy. For the dipolar particles, the smallest radius estimated from the condition Ud>max - Ud>min = kBT is 34 nm. The last estimate presumes that the surface-anchoring strength W is sufficiently strong to maintain the dipolar configuration, namely, K/W ~ 10 nm or so. For the quadrupolar particles, placed in a pattern with l = 80 mm, the minimum size is about 220 nm. The feature of a relatively large size limit follows from the additional factor R/l in the quadrupolar potential as compared to the dipolar potential. By making the pattern with a smaller period l, one can decrease the minimum size of the controllable particles below the estimates above. However, to accurately describe the case of smaller l, one needs to account for (i) the finite size of the particles and (ii) the finite surface-anchoring strength at the surface of the spheres and at the bounding plates.

The presented experimental data on dynamics and placement of the colloidal particles are well described by the elastic model based on the Frank-Oseen functional and the dipolar and quadrupolar type of director deformations around the particles. One of the reasons why the experiment and the theory match well is the relatively large size of the colloids (5 mm). For submicrometer particles and for high concentration of particles, the description might require a more elaborate approach that would, for example, account for the finite anchoring strength. Furthermore, the theoretical description might be supplemented by the consideration of thermal fluctuations. In the uniformly aligned LC environment, the displacement of particles caused by Brownian motion is both anisotropic (35) and, at time scales on the order of the relaxation time of the director, anomalous (36). It would be of interest to explore how the prepatterned director deformations might influence the regimes of colloidal Brownian motion.


To conclude, the photopatterning surface alignment offers a versatile method to control colloidal assemblies by directing particles into predesigned locations and trapping them in these locations with forces that exceed those of entropic nature. Besides the practical utility, the study clearly demonstrates and experimentally quantifies the mechanisms of selective interactions between the colloids and the surrounding landscape of a distorted director. Here, we demonstrated the principle using spheres and a relatively simple periodic director pattern with alternating splay and bend. However, it is obvious that the similar approach can be extended to particles of different shape or size and to more complex 2D and even 3D patterns of the director field. We also expect that the approach is not limited by just normal and tangential

CD —i

anchoring at the surface of the colloidal particles; tilted anchoring is also expected to show selectivity in equilibrium placement within a patterned director field. Colloidal placement in the photopatterned templates can be realized at practically any spatial scale without mechanical stimuli (for example, shearing, pressure gradients) or topographically modified surfaces (37). The described ability to determi-nistically predesign locations of colloidal particles of different types in self-assembled patterns can be of importance for the development of microfluidic, electrooptical, and sensing devices. One of the applications might be sorting of particles with different surface properties and chemical composition, which we currently explore.


We used the nematic 5CB, MLC6815, and an LC mixture with zero dielectric anisotropy (|Ae| < 10-3) formulated by MLC7026-000 and E7 (in weight proportion of 89.1:10.9) (all materials are purchased from EM Industries). The LCs were doped with 0.01 weight % of polystyrene colloids with a radius of 2.5 mm (Duke Scientific). Untreated colloids yielded tangential anchoring. The colloids treated with octadecyl-dimethyl-(3-trimethoxysilylpropyl) ammonium chloride (DMOAP) produced perpendicular director alignment and dipolar structures with a hyperbolic hedgehog on one side of the sphere. The colloidal dispersion in the LC was injected into the photopatterned cell with a thickness of h = 20 mm at room temperature (22°C).

Optical microscopy

We used the Nikon TE2000 inverted microscope with both immersion [100x; Fluor; numerical aperture (NA), 1.3] and dry (60x; M Plan ELWD; NA, 0.7) objectives. Immersion objective was used with the optical tweezers (Verdi V6 green laser; 532 nm, 10 mW) mounted with the microscope to move the particle off the equilibrium state. Polarizing microscopic images were captured using the Energent HR-20000 camera, and the dynamics images were taken by the monochromatic MotionBlitz EoSens mini2 camera with different frame rates (1 to 20 fps).

PolScope microscopy

The director fields produced by photopatterning were established with the help of a microscope (Nikon E600) equipped with the Cambridge Research Abrio LC-PolScope package. The PolScope uses a monochromatic illumination at 546 nm and maps optical retardance and orientation of the optical axis (38).


Supplementary material for this article is available at content/full/2/9/e1600932/DC1

fig. SI. Self-assembly of colloids with homeotropic anchoring in the regions of photopatterned splay distortions during the filling of the LC into the cell.

movie S1. Colloidal sphere with dipolar director is transported from one splay region to the next splay distortion due to elastic interaction.


1. G. M. Whitesides, B. Grzybowski, Self-assembly at all scales. Science 295,2418-2421 (2002).

2. V. N. Manoharan, Colloidal matter: Packing, geometry, and entropy. Science 349,1253751 (2015).

3. P. Poulin, H. Stark, T. C. Lubensky, D. A. Weitz, Novel colloidal interactions in anisotropic fluids. Science 275, 1770-1773 (1997).

4. M. Yada, J. Yamamoto, H. Yokoyama, Direct observation of anisotropic interparticle forces in nematic colloids with optical tweezers. Phys. Rev. Lett. 92, 185501 (2004).

5. I. I. Smalyukh, O. D. Lavrentovich, A. N. Kuzmin, A. V. Kachynski, P. N. Prasad, Elasticity-mediated self-organization and colloidal interactions of solid spheres with tangential anchoring in a nematic liquid crystal. Phys. Rev. Lett. 95, 157801 (2005).

6. I. Musevic, M. Skarabot, U. Tkalec, M. Ravnik, S. Zumer, Two-dimensional nematic colloidal crystals self-assembled by topological defects. Science 313, 954-958 (2006).

7. M. Ravnik, M. Skarabot, S. Zumer, U. Tkalec, I. Poberaj, D. Babic, N. Osterman, I. Musevic, Entangled nematic colloidal dimers and wires. Phys. Rev. Lett 99, 247801 (2007).

8. U. Ognysta, A. Nych, V. Nazarenko, I. Musevic, M. Skarabot, M. Ravnik, S. Zumer, I. Poberaj, D. Babic, 2D interactions and binary crystals of dipolar and quadrupolar nematic colloids. Phys. Rev. Lett. 100, 217803 (2008).

9. T. A. Wood, J. S. Lintuvuori, A. B. Schofield, D. Marenduzzo, W. C. K. Poon, A self-quenched defect glass in a colloid-nematic liquid crystal composite. Science 334, 79-83 (2011).

10. A. Martinez, H. C. Mireles, I. I. Smalyukh, Large-area optoelastic manipulation of colloidal particles in liquid crystals using photoresponsive molecular surface monolayers. Proc. Natl. Acad. Sci. U.S.A. 108, 20891-20896 (2011).

11. B. I. Senyuk, Q. Liu, S. He, R. D. Kamien, R. B. Kusner, T. C. Lubensky, I. I. Smalyukh, Topological colloids. Nature 493, 200-205 (2013).

12. H. Mundoor, B. Senyuk, I. I. Smalyukh, Triclinic nematic colloidal crystals from competing elastic and electrostatic interactions. Science 352, 69-73 (2016).

13. C. Blanc, D. Coursault, E. Lacaze, Ordering nano- and microparticles assemblies with liquid crystals. Liq. Cryst. Rev. 1, 83-109 (2013).

14. O. P. Pishnyak, S. Tang, J. R. Kelly, S. V. Shiyanovskii, O. D. Lavrentovich, Levitation, lift, and bidirectional motion of colloidal particles in an electrically driven nematic liquid crystal. Phys. Rev. Lett. 99, 127802 (2007).

15. Y. Luo, F. Serra, K. J. Stebe, Experimental realization of the "lock-and-key" mechanism in liquid crystals. Soft Matter 12, 6027-6032 (2016).

16. J. S. Lintuvuori, A. C. Pawsey, K. Stratford, M. E. Cates, P. S. Clegg, D. Marenduzzo, Colloidal templating at a cholesteric-oil interface: Assembly guided by an array of disclination lines. Phys. Rev. Lett. 110, 187801 (2013).

17. G. Foffano, J. S. Lintuvuori, A. Tiribocchi, D. Marenduzzo, The dynamics of colloidal intrusions in liquid crystals: A simulation perspective. Liq. Cryst. Rev. 2, 1-27 (2014).

18. D. Voloschenko, O. P. Pishnyak, S. V. Shiyanovskii, O. D. Lavrentovich, Effect of director distortions on morphologies of phase separation in liquid crystals. Phys. Rev. E 65, 060701 (2002).

19. D. K. Yoon, M. C. Choi, Y. H. Kim, M. W. Kim, O. D. Lavrentovich, H. T. Jung, Internal structure visualization and lithographic use of periodic toroidal holes in liquid crystals. Nat. Mater. 6, 866-870 (2007).

20. D. Pires, J.-B. Fleury, Y. Galerne, Colloid particles in the interaction field of a disclination line in a nematic phase. Phys. Rev. Lett. 98, 247801 (2007).

21. H. Yoshida, K. Asakura, J. Fukuda, M. Ozaki, Three-dimensional positioning and control of colloidal objects utilizing engineered liquid crystalline defect networks. Nat. Commun. 6, 7180 (2015).

22. X. Wang, D. S. Miller, E. Bukusoglu, J. J. de Pablo, N. L. Abbott, Topological defects in liquid crystals as templates for molecular self-assembly. Nat. Mater. 15, 106-112 (2016).

23. M. Mitov, C. Portet, C. Bourgerette, E. Snoeck, M. Verelst, Long-range structuring of nano-particles by mimicry of a cholesteric liquid crystal. Nat. Mater. 1, 229-231 (2002).

24. S. Sacanna, W. T. M. Irvine, P. M. Chaikin, D. J. Pine, Lock and key colloids. Nature 464, 575-578 (2010).

25. D. Nykypanchuk, M. M. Maye, D. van der Lelie, O. Gang, DNA-guided crystallization of colloidal nanoparticles. Nature 451, 549-552 (2008).

26. A. F. Demirors, P. P. Pillai, B. Kowalczyk, B. A. Grzybowski, Colloidal assembly directed by virtual magnetic moulds. Nature 503, 99-103 (2013).

27. F. Q. Fan, K. J. Stebe, Size-selective deposition and sorting of lyophilic colloidal particles on surfaces of patterned wettability. Langmuir 21, 1149-1152 (2005).

28. Y. Guo, M. Jiang, C. Peng, K. Sun, O. Yaroshchuk, O. Lavrentovich, Q.-H. Wei, High-resolution and high-throughput plasmonic photopatterning of complex molecular orientations in liquid crystals. Adv. Mater. 28, 2353-2358 (2016).

29. C. Peng, Y. Guo, C. Conklin, J. Vinals, S. V. Shiyanovskii, Q.-H. Wei, O. D. Lavrentovich, Liquid crystals with patterned molecular orientation as an electrolytic active medium. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 92, 052502 (2015).

30. I. Lazo, C. Peng, J. Xiang, S. V. Shiyanovskii, O. D. Lavrentovich, Liquid crystal-enabled electro-osmosis through spatial charge separation in distorted regions as a novel mechanism of electrokinetics. Nat. Commun. 5, 5033 (2014).

31. Y. Gu, N. L. Abbott, Observation of Saturn-ring defects around solid microspheres in nematic liquid crystals. Phys. Rev. Lett. 85, 4719-4722 (2000).

32. H. Stark, Physics of colloidal dispersions in nematic liquid crystals. Phys. Rep. 351, 387-474 (2001).

—i 3

33. I. Musevic, Optical manipulation and self-assembly of nematic colloids: Colloidal crystals and superstructures. Liq. Cryst. Today 19, 2-12 (2010).

34. V. M. Pergamenshchik, Elastic multipoles in the field of the nematic director distortions. Eur. Phys. J. E 37, 121 (2014).

35. J. C. Loudet, P. Hanusse, P. Poulin, Stokes drag on a sphere in a nematic liquid crystal. Science 306, 1525 (2004).

36. T. Turiv, I. Lazo, A. Brodin, B. I. Lev, V. Reiffenrath, V. G. Nazarenko, O. D. Lavrentovich, Effect of collective molecular reorientations on Brownian motion of colloids in nematic liquid crystal. Science 342, 1351-1354 (2013).

37. N. M. Silvestre, Q. Liu, B. Senyuk, 1.1. Smalyukh, M. Tasinkevych, Towards template-assisted assembly of nematic colloids. Phys. Rev. Lett. 112, 225501 (2014).

38. M. Shribak, R. Oldenbourg, Techniques for fast and sensitive measurements of two-dimensional birefringence distributions. Appl. Optics 42, 3009-3017 (2003).

Acknowledgments: We thank J. de Pablo, O. Tovkach, R. Zhang, and Y. Zhou for useful discussions. Funding: This work was supported by the National Science Foundation grants DMR-

1507637, DMR-1121288, and CMMI-1436565. Author contributions: Conceived and designed the experiments: O.D.L., C.P., and Q.-H.W. Performed the experiments: C.P., T.T., and Y.G. Analyzed the data: C.P., T.T., S.V.S., and O.D.L. All authors participated in discussing and writing the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Submitted 29 April 2016 Accepted 11 August 2016 Published 16 September 2016 10.1126/sciadv.1600932

Citation: C. Peng, T. Turiv, Y. Guo, S. V. Shiyanovskii, Q.-H. Wei, O. D. Lavrentovich, Control of colloidal placement by modulated molecular orientation in nematic cells. Sci. Adv. 2, e1600932 (2016).

Q. CD Q.

f —i

—i 3

Control of colloidal placement by modulated molecular orientation in nematic cells

Chenhui Peng, Taras Turiv, Yubing Guo, Sergij V. Shiyanovskii, Qi-Huo Wei and Oleg D. Lavrentovich (September 16, 2016) Sci Adv 2016, 2:. doi: 10.1126/sciadv.1600932

This article is publisher under a Creative Commons license. The specific license under which this article is published is noted on the first page.

For articles published under CC BY licenses, you may freely distribute, adapt, or reuse the

article, including for commercial purposes, provided you give proper attribution. D

For articles published under CC BY-NC licenses, you may distribute, adapt, or reuse the article =

for non-commerical purposes. Commercial use requires prior permission from the American g

Association for the Advancement of Science (AAAS). You may request permission by clicking e

here. cl

The following resources related to this article are available online at C (This information is current as of October 30, 2016): t

Updated information and services, including high-resolution figures, can be found in the a

online version of this article at: v n

Supporting Online Material can be found at: s i

This article cites 38 articles, 9 of which you can access for free at: 8 a

Science Advances (ISSN 2375-2548) publishes new articles weekly. The journal is published by the American Association for the Advancement of Science (AAAS), 1200 New York Avenue NW, Washington, DC 20005. Copyright is held by the Authors unless stated otherwise. AAAS is the exclusive licensee. The title Science Advances is a registered trademark of AAAS