Eur. Phys. J. E 28, 99-111 (2009) DOI 10.1140/epje/i2008-10360-1

THE EUROPEAN PHYSICAL JOURNAL E

Regular Article

Effective interactions of colloids on nematic films

M. Oettel1'a, A. Domínguez2, M. Tasinkevych3, and S. Dietrich3'4

1 Johannes-Gutenberg-Universitat Mainz, Institut für Physik, WA 331, D-55099 Mainz, Germany

2 Física Teórica, Universidad de Sevilla, Apdo. 1065, E-41080 Sevilla, Spain

3 Max-Planck-Institut für Metallforschung, Heisenbergstr. 3, D-70569 Stuttgart, Germany

4 Institut für Theoretische und Angewandte Physik, Universitat Stuttgart, Pfaffenwaldring 57, D-70569 Stuttgart, Germany

Received 16 April 2008

Published online: 9 September 2008 - © EDP Sciences / Societa Italiana di Fisica / Springer-Verlag 2008

Abstract. The elastic and capillary interactions between a pair of colloidal particles trapped on top of a nematic film are studied theoretically for large separations d. The elastic interaction is repulsive and of quadrupolar type, varying as d-5. For macroscopically thick films, the capillary interaction is likewise repulsive and proportional to d-5 as a consequence of mechanical isolation of the system comprised of the colloids and the interface. A finite film thickness introduces a nonvanishing force on the system (exerted by the substrate supporting the film) leading to logarithmically varying capillary attractions. However, their strength turns out to be too small to be of importance for the recently observed pattern formation of colloidal droplets on nematic films.

PACS. 82.70.Dd Colloids - 68.03.Cd Surface tension and related phenomena - 61.30.-v Liquid crystals

1 Introduction

The interactions of colloidal particles trapped at fluid interfaces have been found to differ significantly from the corresponding interactions in bulk solvents. This has been studied mostly for electrically charged particles trapped at interfaces with water. On one hand, the presence of the interface gives rise to direct dipolar electrostatic repulsions between the colloids (see Refs. [1-3] for some experimental evidence), on the other hand deformations of the interface may induce longer-ranged capillary attractions (briefly reviewed in Refs. [4-6]) which is possibly the source of pattern formation observed in various experiments [7-11]. (See, however, Ref. [12] for an alternative explanation due to interface impurities.)

Recently [13,14], the experimental observation of ordered structures of glycerol droplets bound to a nematic-air interface has been reported and attributed to an effective pair potential between the colloids which contains a repulsive, elastic part due to director deformations in the supporting nematic film and an attractive, capillary part which is long-ranged and mediated by logarithmically varying deformations of the nematic-isotropic interface caused by the droplets. The schematic setup of this experiment is depicted in Figure 1. According to reference [13], the colloidal particles experience an upward force caused by elastic forces due to director deformations in the supporting nematic film. This upward force on the colloids

a e-mail: oettelm@uni-mainz.de

is supposed to give rise to the aforementioned logarithmically varying interface deformation. Applying a superposition approximation for the deformation field, one can show that the ensuing effective capillary interaction potential between two colloids is likewise varying logarithmically. This is similar to the flotation interaction of mm-sized particles at fluid interfaces for which the force on the colloids is caused by gravity (see, e.g., Ref. [15]) and also parallels the tentative explanation given for the experimentally observed attractions between sub-^m charged colloids at a water-oil interface [16] (for the controversy around this explanation see, Refs. [17-24]).

However, it is now well established [17-19,21-23,25] that interface deformations and effective colloidal interactions varying logarithmically only arise in experimental systems which are not isolated mechanically. For mechanically isolated systems it can be shown [26-29] that both the interface deformation around a single colloid and the effective interaction between two of them are shorter-ranged and the latter cannot be calculated reliably within the superposition approximation.

In the following we will extend the arguments presented in references [21,26-29] to systems with colloids at nematic interfaces. We will show that mechanical isolation of the system "nematic film-colloid-air" can be violated through a subtle interplay between the finite thickness of the film and the anchoring conditions at the colloids and at the nematic interfaces with the substrate and with the air, respectively. However, for experimental conditions as the ones described in reference [13], a quantitative estimate

of the strength of the ensuing logarithmic attraction between the colloids yields that these attractions are unob-servably small. Therefore it seems likely that this kind of asymptotic capillary forces cannot be invoked as a relevant mechanism to account for the observations reported in reference [13].

The manuscript is organized as follows: In Section 2 the coarse-grained model for the nematic phase is introduced which will serve as the basis for all subsequent calculations. In Section 2.1 we study the case of an infinitely thick nematic film. First, we compute the asymptotic form of the director field and the ensuing elastic force between two particles. Then we calculate the effective force arising from the deformation of the fluid-nematic interface caused by the elastic stresses. In Section 2.2 we consider a nematic film of finite thickness, which models more closely the experimental setup described in reference [13]; for such a system we extend the above calculations to the two opposite cases of perpendicular and parallel anchoring of the director field at the substrate surface. In Section 3 we discuss our results.

2 Coarse-grained model

In view of the mesoscopic length scales involved, we describe the bulk part of the nematic free-energy associated with the director deformations in terms of the Frank free-energy expression within the one-coupling approximation [30]

= / d3rf b(r)

= f / d3r [(V • n)2 + (Vx n)2]

2 JVne

K f 3 = — d3r Vni • Vni

2 JVne

+K f d3r V- [n(V • n) - (n • V)n].

2 JVne

Vne denotes the volume occupied by the nematic film, n is the director field (n2 = 1), and the constant K is of the order of 10-11 N [30]. The total divergence term in the last line of equation (1) exhibits the so-called "K24-structure" and is unimportant for the bulk equations describing the equilibrium configuration. The surface free energies associated with the interfaces with air and substrate, respectively, (see Fig. 1) are described in terms of the Poulini expression [31]

Fne 2 J

dA(n • eA)2

Lair-ne

i dA(n • eA)2.

' Asub-ne U Acoll-ne

Here, eA denotes the local surface normal unit vector pointing outwards from the film or the colloid. Normal alignment is favored for Wi < 0 and parallel alignment for

glycerol substrate

Fig. 1. Schematic setup of the experiment reported in reference [13]. Colloidal glycerol drops (R = 1... 7 ^m) are trapped at the surface of a thick nematic film (h ~ 60 ^m). The director field in the nematic film is sketched by the black lines and their dotted interpolations. Nematic anchoring at the glycerol substrate at the bottom of the film as well as on the surface of the glycerol drops is parallel, whereas it is perpendicular (homeotropic) at the deformed nematic-air interface. Each drop is necessarily accompanied by a topological defect (®).

Wi > 0. Typically one has W | - 10-5 N/m [31] so that the length scale K/|Wi | — 1 ^m is smaller than the range of droplet radii investigated in reference [13]. Thus, in the "strong anchoring" limit which we shall consider, the effect of the boundary terms is so strong that as a first approximation it amounts to fixing the angle between the director and the surface normal. We shall adopt W1 < 0 (normal alignment at the nematic-air interface) and W2 > 0 (parallel alignment at the nematic-glycerol interfaces). Some consequences of the deviations from the strong-anchoring limit will be discussed in Appendix A.

These surface contributions to the free energy ("wetting energies") are small corrections to the surface tensions which are mainly due to dispersion interactions. We denote these non-nematic contributions to the surface tensions as 71 (colloid-air surface tension), y2 (substrate-nematic surface tension) and 7' (nematic-air surface tension). Typically, these surface tensions are of the order of 10-2 N/m. Therefore, they are much larger than the constants Wi which determine the nematic contributions to the surface tension. We see that due to the above anchoring conditions, the surface tension of the substrate-nematic interface carries no nematic contributions due to n• eA =0 (Eq. (2)), whereas in the strong-anchoring limit the full nematic-air surface tension is 7 = 7' + W1/2 « 7' due to n • eA = 1 (Eq. (2))1.

1 A genuine contribution to the surface tension of a nematic interface arises if one takes into account the variation of the nematic tensorial order parameter through the interface, described by, e.g., the Landau-de Gennes free-energy functional (generalizing Eq. (1)). The magnitude of these contributions can be estimated by the surface tension of the interface between

The canonical stress tensor nij associated with the free-energy expression in equation (1) is given by

nk,i ^ij f 7

where nk i = dnk/dxi (summation over k). The total stress tensor nij is obtained by adding the contribution of the isotropic pressure p:

nij nij p.

2.1 Macroscopically thick nematic film

First we consider the limiting case h (i.e., very thick

nematic films (see Fig. 1))2. Due to the small values of the elastic coupling constant, K ^ yR and of the anchoring energy, \Wi\ ^ y, the equilibrium configuration of a single colloid at the nematic-air interface deviates only slightly from the reference configuration depicted in Figure 2. In this latter configuration, the interface is flat and the colloid is positioned such that the contact angle fulfills Young's law cos 0 = (yi — Y2)/y-

The total force on the whole system reads (the superscript denotes evaluation on the positive (negative) side of the oriented surface, i.e., on the side the arrows in Figure 2— which indicate the surface normals— point to (do not point to)):

/ dV(V • n)

Jvi u V2

+ / dA • (n+ - n-)

J Smen U S1

- dA(nzz + Pair - p) ez

J Smen

+ / dA • [n +(pair - p)1] ■J Si

■.ez +

dA • n.

In obtaining this equation we have applied Gauss' theorem. Furthermore we have used the relation V- n = 0 in volumina Vi and V2 which is valid because the reference configuration is taken to be in force equilibrium. This also implies that the isotropic pressures above the interface

the nematic and the isotropic phase of a liquid crystal. Typically, such a surface tension is also 0(10-5 N/m) ~ |Wi| and therefore small compared to the dispersion force contribution. Furthermore, a distorted director structure in the bulk may also give rise to surface energy contributions on the boundaries. For an example, see reference [30], p. 131 and p. 174. Also in this case it can be argued that the corresponding contributions to the surface tension do not exceed |Wi|.

2 This was implicitly assumed also by the authors of reference [13] in discussing Figure 2 therein.

nematic phase

Fig. 2. In the reference configuration the whole system is divided into volumes Vi and V2. Volume V2 (enclosed by the upper dashed curve) includes the air and the glycerol drop and volume V1 (enclosed by the lower dashed curve) includes the nematic. The arrows indicate the direction in which the surfaces (including the infinitesimally displaced ones) are oriented: S encloses the whole system, Smen is the interface between the nematic phase and air (acting as a meniscus), and S1 is the interface between the colloidal drop and the nematic phase. The director field and the topological defect (®) are indicated as in Figure 1. The radius of the three-phase contact line is denoted by p0 = R sin 6, where 6 is the contact angle of the air-nematic interface with the colloid of radius R. p denotes the lateral distance from the vertical symmetry axis of the colloidal drop.

(pair) and below it (p) are equal and that the director configuration is given by the corresponding Euler-Lagrange equilibrium equations following from the functional in equation (1). Since at the interface S1 the colloidal drop is rigidly attached to the liquid crystal, we can identify the vertical force F on the colloid and the total force Fn on the air-nematic interface by

dA • n

respectively. Mechanical isolation of the system means that the total force j>S dA • n acting on it is zero which leads to

F = Fn. (8)

For a given force on the colloid and a given stress on the interface, the interface deformation relative to the reference configuration can be determined perturbatively. To that end, we summarize briefly those results of reference [21] which are pertinent also for the present system.

1 n —

With the introduction of the two small, dimensionless parameters

eF = -F/(2npo Y) and ^ = /(2npo y), (9)

one can expand (up to second order in eF and en ) the free-energy difference F associated with the interface deformation u(p > p0) around a single colloid (see Fig. 2) and with a vertical shift which is the difference of the colloid center position relative to that in the reference configuration:

F ~ 2^7 / dp p

■' P0

1 /du\ u

2 VdpJ + 2A2 y

+ -nzz U

+n7[u(p0) - 4h]2 - F^h.

M CA 1 r

u(p) ~ po (ew - £f)ln + da

P Y ./ p

Tzz(a)ln

on similar considerations, we shall show that for a colloidal drop located at the air-nematic interface the boundary conditions for that interface impose a quadrupole-like asymptotic behavior of the director field. Macroscopically far from the colloid the director is oriented parallel to the z-axis. Accordingly, at large but finite distances r the director is given by n(r) ~ (ni,n2,1 — O(n2,n2)) and the bulk free energy corresponding to equation (1) is given by

Fne " 2

(Vnj)2 + O(n4)

Here, A = (Y/(gpm))1/2 is the capillary length associated with the interface where g is the gravitational constant and pm is the mass density of the nematic phase. This expression for the free energy contains all surface free-energy changes relative to the reference configuration involving the interfaces between air, nematic, or colloid. It also contains the contributions due to volume forces acting on the nematic (associated with A) and the energy change of the colloid upon vertical shifts (for further details see Ref. [21]). Note that to leading (quadratic) order in en, eF the free-energy change of the nematic due to the shifted interface and due to a change in the director configuration with respect to the reference configuration is captured by the term « J" nzz u. (The analogous textbook argument for electrostatics [32] can be easily generalized to the nematic case described by the free-energy expression in Eq. (1).)

Minimizing F with respect to u(p) and and focussing on the regime p ^ A yields

Here we have discarded the total divergence term in the free-energy expression (1). It is unimportant for the bulk equations and it adds a mere constant to the free energy because the director is anchored normally at the boundary (the nematic-air interface). Thus for each component i = 1, 2 the equilibrium director field fulfills the Laplace equation

An = 0. (13)

Analogously to electrostatics, the asymptotic solution for ni can be expanded in terms of multipoles. To this end, we consider the reference configuration in Figure 2. We choose as the origin of the coordinate system the center of the circle formed by the planar three-phase contact line. The solution for the director field has to fulfill the following requirements: i) rotational covariance around the z-axis3 and ii) n,(x, y, z = 0) = 0. Analyzing the multipole ansatz (with r = (ri,r2,r3))

r + Pia r3 + ID Qi«ß

■-rß

+..., (14)

a,ß=1

with C ~ 1.12. We see that in the case of an isolated system (en = eF) the logarithmic part of u(p) vanishes. The second term on the rhs of equation (11) leads to u(p ^ to) « p-n+2 if nzz « p-n and thus describes a shorter-ranged power law decay of the interface deformation.

The absence of logarithmic deformations for an isolated system has been derived here under certain simplifying conditions (small interfacial deformation everywhere, rotational symmetry). In Appendix B we demonstrate that this conclusion holds in general.

2.1.1 Asymptotic director configuration and elastic force between colloids

In references [33,34] it has been shown that a colloidal drop immersed in the bulk of a liquid crystal is accompanied by a single counterdefect such that the total topologi-cal charge is zero (here, the volume occupied by the colloid contains a topological charge which may be represented by a virtual defect inside the colloid) and the asymptotic behavior of the director field is of dipolar character. Based

it follows that rotational covariance requires q = 0,

Pia P $ia + Pmag ^ia3 and Qiaft Qft ^ia + Qmag,^ ^ia3

(ejk is the Levi-Civita tensor). Pmag and Qmag ft are dipole and quadrupole moments, respectively, for a director field of "magnetic" type, i.e., for which divni = 0 holds. The boundary condition ii) at the interface with the air further imposes P = Pmag = 0 and Qft = Q ¿g3, Qmag ft = Qmag ¿ft3. It appears to be difficult geometrically to match the asymptotic solution of "magnetic" type with a solution near the colloid which obeys parallel anchoring at the colloid surface. Therefore we discard the magnetic quadrupole, i.e., the leading asymptotic term is given by the remaining, "electric" quadrupole term

Q^? +...

(z = 13 ).

This is at variance with reference [13] (see Eq. (3) therein and the considerations in the paragraph above that equation which assume a dipole field) but it is consistent with the analysis in reference [34]. Dimensional analysis yields Q = O(R3) [34]. Note that we have derived the asymptotic behavior of the director field at the interface using

3 If D specifies the transformation matrix for such a rotation then this requirement is given by n(D • r) = D • n(r).

strong anchoring at the interface (ni = n2 = 0). In Appendix A we discuss corrections to strong anchoring which are, however, of subleading character and leave the leading behavior (Eq. (15)) unchanged.

The asymptotic elastic interaction between two colloids in the bulk at distance d accompanied by a quadrupo-lar director deformation has been analyzed in reference [35] (for weak anchoring) and in references [36,34] (using a coarse-graining method, applicable also for strong anchoring) yielding identical results. For the present configuration (distance vector perpendicular to the asymptotic director) the elastic potential is repulsive and varies

V KQ2 2

Vel « d5 « YPo £F

KQ2 1 2 Q p8 :

(p2 = r2 + r2). Consequently the interface deformation around a single colloid for a mechanically isolated system obeys u(p ^ to) « p-6 (see Eq. (11)).

For the problem of two identical colloids located at p1 and p2 (vectors are defined in the interface plane z = 0) a distance d = |p1 — p2| ^ R apart, the expression for the free energy is a straightforward generalization of the one for the single-colloid free energy given in equation (10):

+ ^ {2^/ d^hi — u]2 — Fi^hA .

¿=1,2 ^ Po dSi >

Here, F denotes the force on colloid i and 4hi is the relative position of its center. The integration domain Smen

is the whole interface plane except for the two circular disks bordered by the (reference configuration) contact lines dSi. The meniscus-induced effective potential is the difference between the equilibrium free energy of the two colloids at distance d and their free energy at macroscopic distance:

Vmen(d) = Feq(d) — Feq(d ^ to).

We have used that the dimensionless force parameter eF is proportional to Q2 which actually follows from equations (6, 8, 9), and (17) below. Note that we have simply extrapolated the bulk results for two interacting colloids which cause asymptotically quadrupolar deformations of the director. This appears to be reasonable because the asymptotic director field in the nematic phase for the interface problem is precisely that of the bulk solution in the lower half-plane and because the bulk solution is antisymmetric with respect to z ^ —z, thus respecting the boundary condition nj(x, y, z = 0) = 0. However, the precise numerical value of the quadrupole moment Q might be rather different for the case of the colloid trapped at the interface as compared to the bulk case.

2.1.2 Asymptotic behavior of the stress on the interface and meniscus-induced effective potential between colloids

The asymptotic behavior of the stress tensor component nzz at the interface follows from inserting equation (15) into equation (3):

nzz | interface = 2 (n»,z "¿,^1 "j,^)

As before, minimization with respect to U(p) and renders the equilibrium free energy.

The behavior of Vmen(d) has been analyzed in detail in references [21,29]. Here we summarize these results as far as they are relevant for the present problem. The interfacial stress nzz may be decomposed generally as

nzz(p) = nzz(|p — Pl|) + nzz(|p — P2I) +2 nzz,m(p)

= nzz,1 + nzz,2 +2 nzz,m- (20)

Here, i denotes the stress around colloid i which pertains to the problem of a single colloid. To quadratic order the asymptotic director field around two colloids is given by the superposition of the components n of the single-colloid solutions and thus to this order we recover the decomposition in equation (20) with the mixed component of stress field nzz,m given by

n = KQ2 (P — Pl) • (P — P2) (21)

nzz,m 0 Q 1 is 1 15 ' (21)

' 2 |P — Pl|5 |P — P2|5

It turns out that for a system under an external force (en = eF) the mixed term nzz,m does not contribute to the leading term in Vmen. Thus in the case of a non-vanishing external force this leading contribution to Vmen is obtained by a superposition ansatz which consists in approximating the interfacial deformation and the total stress field by the sum of the respective single-colloid quantities only (nzz « nzz'i + nzz'2, U ~ «1 + «2) [21,29]:

Vmen(po < d < A) ~ — 2n Y — £F)2 ln — (22)

(en = £f )•

For an isolated system (en = eF) and for the stress given in equation (17) the superposition approximation nzz ~ nzZ'i + nzz,2 yields Vmen x e^/d8 as the dominant term. However, this is not the leading term, which rather stems from nzz,m. This term has two peaks around the colloid centers and therefore close to the colloids it can be approximated by

E(—1)

ed ■ (p—Pi) |p—Pi|5

where ed = (p2 — P1)/d. As discussed in reference [29], the qualitative behavior of Vmen(d) is captured by the integral

Vmen(d) « YP^f/ d2pnzz,m(p) « YP2eF (po) (24)

J Smen

(en = eF )•

over n

We note that this form of the mixed stress and thus the leading behavior of Vmen(d) is formally analogous to the contributions of the electric-field components parallel to the interface in the case of charged colloids [29]. The meniscus-induced potential Vmen is repulsive and falls off asymptotically4 <x d-5, as does the likewise repulsive elastic potential Vel. We note that Vei < |eF | and Vmen < eF, so that the meniscus-induced potential is small compared to the elastic one for the parameters of the experiment reported in reference [13].

In order to explain the experimentally observed attractions, in reference [13] a perturbative picture similar to the one presented above was suggested, but the effect of the interface stress nzz was neglected completely. (A similar error has been made in Ref. [16] which the authors of Ref. [13] refer to.) In a heuristic way, only an "upward" force on the colloid (perpendicular to the interface) associated with the anchoring "wetting" energy at the nematic-particle interface has been invoked in reference [13], neglecting the force on the interface described by nzz. In this way, the unbalance of the force gives rise to a logarithmic term in the meniscus deformation and in the meniscus-induced potential. However, mechanical isolation (i.e., force balance) renders the meniscus-induced potential actually repulsive and shorter-ranged (see Eq. (24)), apparently in contrast to the experimental results. A logarithmically varying potential can only arise if mechanical isolation is violated (see App. B). Below we shall investigate whether a net force on the system "colloid and interface" may appear if the thickness of the nematic phase is finite, as it is the case in the actual experiment.

2.2 Finite thickness of the nematic film

In our discussion of a finite film thickness of the nematic phase we shall consider two cases:

1) The anchoring of the nematic director at the surface of the bottom substrate is perpendicular as it is the case at the upper interface with the air. This case bears a strong formal resemblance to charged colloids on water surfaces which have been discussed in references [21,29].

2) The anchoring at the bottom substrate surface is parallel (as in the experiment reported in Ref. [13]). At large lateral distances from the colloids, this leads to a director field which gradually rotates from parallel orientation at the bottom substrate to the perpendicular orientation at the upper interface.

For both cases the total force on the system —comprising air, nematic film, colloid, and the substrate— must be zero. This leads to (see Fig. 3(a) and compare with

Fig. 3. The reference configuration for a nematic film (a) differs from the reference configuration shown in Figure 2 by the addition of the substrate volume V2'- The interface between the substrate and the nematic film is denoted by Ssub. Panel (b) shows the string of image quadrupoles Qb = Q which are needed to fulfill the boundary conditions for perpendicular anchoring of the director field (originating from the colloid quadrupole Q) at both interfaces confining the nematic film. The distance between any two nearest-neighbor quadrupoles is 2h.

Eq. (5))

0 = j> dA • n

dV (V-n)

jvi u v2u v2' + / dA • (n+ - n-)

J Smen U S1 U Ssub

/ dA(nzz + Pair - p) ez

J Smen

+ / dA • [n +(pair - p)1] j sx

+ dA(nzz + Psub - p)ez.

Superficially one would expect a leading decay Vmen <x d-

as displayed in equation (23). However, due to the geometric factor in the numerator of equation (23), this apparent leading order vanishes upon integration.

For reasons which will become clear in the discussion of the second anchoring case, we consider the isotropic pressures in the substrate and air, psub and pa;r, respectively, not necessarily to be equal to the pressure p in the nematic film. The second equation in equation (25) follows from the equilibrium condition V-n = 0 which holds in all volumina.

2.2.1 Perpendicular anchoring at both interfaces

As discussed above, the presence of the colloid asymptotically generates a quadrupolar director field which fulfills the boundary condition at the nematic-air interface. In order to fulfill the boundary condition at the substrate-nematic interface, an image quadrupole of the same strength Q is needed which, however, leads to a violation of the nematic-air interface boundary condition and requires a second image quadrupole etc. Continuation of this process leads to a string of image quadrupoles as depicted in Figure 3(b). This string of quadrupoles generates a stress field n which vanishes for large lateral distances.

nematic film

substrate

Therefore the isotropic pressures must be equal in all volumina: psub = pair = p. From equation (25) one finds that the difference between the force F on the colloid and the integrated stress Fn at the nematic-air interface,

F - Fn = 4F = -

^ Ssub

is given by the integrated stress over the substrate surface, i.e., the total force on all quadrupoles above the substrate surface exerted by the image quadrupoles in the substrate. Expressed in terms of the force Fq-q between two quadrupoles at distance 2h, we find

—¿f/fq-q = £ n +£ n6 + ...

En 5 = Z(5) « ^

In this equation, the first sum is the total force (divided by Fq-q) on the first quadrupole above the substrate exerted by all quadrupoles in the substrate, the second sum is the total force (divided by Fq-q) on the second quadrupole above the substrate, etc. For the force between two quadrupoles at distance 2h we find, using equations ( 5, 3), and (26)

F 5 KQ2 fq-q =6nKh6 •

This result can be obtained more easily by solving the field equations Ant = 0 in cylindrical coordinates rather than by using the image quadrupoles.

x-axis. With no colloid present at the nematic-air interface, the equilibrium director field is given by

sin(—qoz) no = I 0

, cos(—qoz),

qo = n/(2h), (29)

with the consequence that both at the nematic-air and at the nematic-substrate interface a constant stress is acting:

For the unperturbed interface to be in equilibrium, the air and substrate pressures differ from the pressure in the liquid crystal: p — psub [air] = no,zz.

We now introduce a single colloid at the nematic-air interface in the reference configuration. Due to this pressure difference, the force on the colloid and the integrated stress over the nematic-air interface are given by (see Eq. (25))

F ez = / dA ■ (n — no,zz 1), ■j s 1

Fn = dA(nzz — no,zz ),

(3 ) (32)

and the total excess force on the system "colloid and nematic-air interface" is determined by

F - Fn = 4F = -

dA(nzz — no ,zz ). (33)

Thus mechanical isolation for the system "colloid and nematic-air interface" is violated and the total force 4F on this system is (up to the factor 1 .04) given by the quadrupole force Fq-q. Nevertheless the magnitude of the corresponding induced logarithmic capillary potential (see Eq. (22) with en — eF = 4F/(2nYp0)) is small, because Q ~ R3, due to dimensional arguments [34]. For parameters similar to the ones appearing in reference [13] (R/h ~ 10-1, K/(yR) - 10-4) we find Vmen ~ 10-14 ln(R/d), which is unimportant for the actual intercolloidal interaction.

The elastic potential Vei(d) between two colloids is the interaction between the second quadrupole and the first quadrupole together with its string of image quadrupoles. Using the solution given in references [36,35] we have checked that Vel (d) remains repulsive. For d < h the overall magnitude of Vel(d) is somewhat weakened, whereas for d ^ h a crossover to Vel(d) <x exp(—d/h) is observed5.

2.2.2 Parallel anchoring at the bottom substrate

We assume that the substrate induces a preferred in-plane axis for the director orientation which we take to be the

In order to calculate the director field n and the stress tensor n in the presence of the colloid, we introduce the auxiliary director deformation fields v(x, y, z) and w(x,y, z) which parametrize the deviations from the unperturbed director field n0 and which are small at large distances from the colloid:

sin(—qoz + v) cos w sin w

, cos(—qoz + v) cos w ,

cos(qoz) v + 2 sin(qoz)(v2 + w2)N w

isin(qoz) v — 2 cos(qoz)(v2 + w2),

+ O((v,w)3).

The first equality in equation (34) is a general parame-trization of the director field n in terms of the auxiliary fields v, w which fulfills n2 = 1. These auxiliary fields are taken to be zero at the substrate and the nematic-air interface, i.e., the boundary conditions are v(x, y, 0) = v(x, y, —h) = w(x, y, 0) = w(x, y, —h) = 0. The nematic free energy of the film up to order O(v2, w2) is obtained by inserting equation (34) into equation (1) for the Frank free energy after dropping the total divergence of the K24-type. Using the boundary conditions for v and w, we obtain [37] (with the notation introduced in Eq. (3))

Fo + ■

i d3r (v2i + w2i — q2 w2) .

Here, F0 is the free energy of the film without colloid. At first glance it is not evident that this free energy is positive definite. However, the boundary conditions on v and w ensure positivity [37]. Upon minimization, we find Av = 0, i.e., the deformation field v of the director can again be expanded in terms of electrostatic multipoles (Eq. (14)). On the other hand, the solution for w must fulfill the Helmholtz equation (A + q2)w = 0 and can be expanded in terms of multipoles as follows:

w(r, 0, =

Jm\u :

j=0 m=-j

P z Q zx

Pv + +

■ + Vi,

The dipole contribution should vanish for h ^to. If one assumes a power law dependence on h, dimensional analysis for Pv leads to

Pv = O(R2 (R/h)K

(with k > 0).

Zy -L (QW J5/2(qor) + QW Y5/2(qor)) + ...

image •

6 In this respect the pictorial argument given in reference [13] (see Fig. 3(b) therein) is slightly misleading (at least in an asymptotic sense): there the assumed director configuration around the colloids for the case of a finite film thickness is a tilted dipole with a nonvanishing component v x Pvxx/z3. Actually, the broken symmetry in the x-direction rather enters through a tilted quadrupole.

We note that there are two quadrupolar contributions (with dimensionless moments QW and QW ) which show an oscillatory behavior for radial distances r ^ h (i.e., q0r ^ 1). On the other hand, for large film thicknesses there is an intermediate asymptotic regime R ^ r ^ h:

w(qoR « qor « 1)= ¡J (qW (q0r)

- QW 3

Qw (qor)3

+ ••• + w

image •

X (w/m Jj+1/2(qor)+w]m Yj+1/2(90r)) , (36)

where J[Y]j+1/2(r) are the spherical Bessel functions of the first [second] kind and Yjm(#, are the usual spherical harmonics for the standard set of spherical coordinates r, 0, (The origin is again taken as the center of the circular three-phase contact line.) The coefficients wjm] are dimensionless multipole moments.

The Dirichlet boundary conditions for v and w at the two interfaces can be fulfilled as before by constructing the full solution in terms of multipoles around the colloid and the corresponding image multipoles as shown in Figure 3. Since rotational covariance is broken by the parallel substrate anchoring, the solution for v contains a nonzero dipole contribution. Nevertheless, the director field still obeys a reflection symmetry with respect to the xz-plane: v(x,y, z) = v(x, —y, z) and w(x,y, z) = —w(x, — y, z). Therefore, the leading asymptotic behavior for v is given by

For large film thicknesses the solution for the director field near the nematic-air interface in the regime R ^ r ^ h should coincide with the solution for macroscopically thick films (Eq. (15)). Near the nematic-air interface, one has v « nx, w « ny. Therefore one recovers the rotationally covariant quadrupole solution of equation (15) for Qv = Q (obtained by equating Eq. (37) with Eq. (15) for i = 1, i.e., r1 = x) and

/q3 = Q (obtained by equating Eq. (37) with Eq. (15) for i = 2, i.e., r2 = y). Since QW - q3Q and Q = O(R3), QW = O([qoR]3) is a very small number.

The Dirichlet boundary conditions for w at the substrate and at the nematic-air interface enforce that the contribution to w due to the quadrupole QW and all corresponding image quadrupoles is zero. This holds also for the contribution of all higher multipoles of degree j for which the radial dependence is given by Jj+1/2(r). The only solution of the Helmholtz equation which fulfills the boundary conditions and which, as an additional requirement, is smooth everywhere, is w = 0. Since the Bessel functions of the first kind are smooth everywhere, all respective multipole moments must be zero. (This does not hold for the multipole moments pertaining to the Bessel functions of the second kind since Yj+1/2(r) is singular at r = 0.)

The excess force on the system "colloid and nematic-air interface" follows from equation (33) and can be expressed as

The precise functional form of Pv turns out to be unimportant for the subsequent calculations. We note that an asymptotic solution with a nonvanishing x-component of the dipole moment cannot appear because it would not fulfill the boundary conditions and the reflection symmetry w(x, y, z) = —w(x, —y, z) excludes any dipolar contribution in the solution for w (see footnote 6). Therefore the leading asymptotic behavior of w takes the form

AF = --/ dA(v2z - 2qo + w2z)

J Ssub

utilizing the boundary conditions for the solutions v and w. For the multipoles appearing in the expansion of v the method described in Subsection 2.2.1 may be used. In order to obtain the quadrupolar contributions to w, we perform the summation over the image multipoles and the integration over the substrate surface numerically. Due to the slow decay of Y5/2(r) « (cos r)/r1/2 the stress integral is superficially divergent. However, a detailed asymptotic analysis yields convergence with the result

-AF = n K

3Z(3) P2 + 5ZM Q2 + h4 +12 h6 + Cw

(QW )2

where cw « 0.30 has been determined numerically. Note that the linear contribution due to the term — q0 v,z in the stress tensor (Eq. (41)) turns out to be zero, as can be easily checked by applying Gauss' theorem and the field equation Av = 0. Using the above considerations

concerning the magnitude of the multipole moments Pv, Qv and QW, we obtain

K ~ Vh

) + aQ

where aP and oq are numerical coefficients of order 1 and k has been defined in equation (38). Again the excess force on the colloid falls off rapidly with increasing film thickness such that the induced logarithmic capillary interaction (see Eq. (22)) remains very weak. For the parameters characterizing the experimental system studied in reference [13] (R/h - 10-1, K/(yR) - 10-4), we find Vmen(d) — 10-11-4K kBT ln(R/d) which appears to be undetectably small. Note that for d < h the direct elastic repulsion remains essentially unchanged because in this regime the leading term of the elastic interaction is given by the repulsion between the two quadrupoles located at the colloid sites. In this case the image quadrupoles can be neglected.

3 Discussion and conclusion

We have investigated the effective potential between two colloidal microspheres of radius R floating at asymptotically large distances d on an interface between a nematic film of thickness h and air (Fig. 1). This effective potential is the sum of an elastic interaction caused by the director distortions around the colloids and of a capillary interaction mediated by surface deformations. We have analyzed the effective potential for large d employing the coarsegrained Frank free energy (within the one constant approximation) for the director distortions and the linearized Young-Laplace equation for the interface distortions.

In the case of a macroscopically thick nematic film, the director deformation around a single colloid is of quadrupolar type. Thus the induced elastic interaction between two colloids at distance d is repulsive and of quadrupolar type <x d-5. The capillary interaction is also repulsive and decays <x d-5 but it is much weaker than the elastic interaction. This rapid decay is a consequence of the mechanical isolation of the system, i.e., of the fact that the net force on the colloidal particles and the surrounding interface vanishes.

A finite film thickness h leads to a violation of the mechanical isolation of the system "colloid and interface". (The excess force is counteracted by the film substrate such that the whole experimental system is of course in mechanical equilibrium.) However, on the thermal energy scale the strength of the ensuing logarithmic capillary potential turns out to be very small. In the case of homeotropic boundary conditions on both sides of the nematic film the strength is proportional to (R/h)12 (see Eq. (28) and the subsequent discussion). For twisted boundary conditions (parallel at the bottom substrate, perpendicular at the upper nematic-air interface) there is a qualitatively different asymptotic behavior of the director field due to a dipole contribution, which is induced by the broken rotational symmetry and vanishes for h ^to.

However, even for these boundary conditions the strength of the logarithmically varying capillary potential remains extremely small, vanishing at least <x (R/h)8 (see Eq. (43) and the subsequent discussion).

Thus our analysis based on the mechanical isolation of the experimental system under consideration rules out a significant attractive contribution — ln d of capillary type to the effective potential between two colloids at a nematic interface. The amplitude of such a logarithmic contribution vanishes rapidly for large film thickness h. Therefore, this effective pair potential does not provide a mechanism for the stability of isolated colloid clusters at a nematic interface as reported in reference [13].

Naturally one strives for other explanations for the observation of stable clusters reported in reference [13]. There the mutual center-to-center distance between neighboring colloids in the cluster has been found to be between 3R and 5R, depending on the radius R of the colloids. If the cluster stability is attributed to a minimum in the effective pair potential at this range 3R... 5R of distances, the applicability of the above asymptotic considerations at such distances is doubtful for both the elastic and capillary contributions to the effective pair potential:

- Elastic part: For the single-colloid problem, the asymptotic, quadrupolar form of the director field (Eq. (15)) is based on the assumption that the deviations of the director from the preferred direction ez are small: |nx[y]| ^ 1. At a radial distance p from the center of the colloid this implies Q/p3 ^ 1, which seems to be fulfilled for the dimensional estimate for the quadrupole moment Q = O(R3) and for the distances p = 3R... 5R under discussion. However, the absolute magnitude of the director deformations is not fixed by these arguments, so that it might be worthwhile to determine the minimum of the Frank free energy by a full numerical calculation in the presence of one or two colloids at the nematic-air interface, similar to reference [38], where the explicit director solution around a single colloid in the bulk has been determined. Our analysis indicates that only if the multipolar expansion of the director field around a single trapped colloid fails for distances p = 3R... 5R, there is a chance for attractive elastic interactions between two colloids at these distances. Failure of the multipole expansion at these intermediate distances would point to a strong deviation of the magnitude of the multipole moments Mn of order n from the dimensional estimate Mn — Rn+1. Thus it is certainly worthwhile to determine multipole renormalizations for bulk and interface configurations, taking into account the full nonlinear-ity of the director field equations.

- Capillary part: Whenever the colloid-induced nematic stress deviates from the asymptotic, quadrupolar form (Eq. (17)), corrections to the asymptotic capillary potential will arise. However, due to the smallness of the dimensionless force eF (see Eq. (9), eF <x K/(yR) ~ 10-3) the perturbative treatment of the ensuing capillary deformations is justified. Therefore, one of our main findings, i.e., that the magnitude of the capillary

potential (« eF) is always smaller than the magnitude of the elastic interaction (« |eF|), is likely to hold also for elastic stresses deviating from the asymptotic limit. Another possible, subtle effect related to the anisotropy inherent to nematic phases allows the substrate to exert also a torque on the colloidal particles and the surrounding interface. This case can be studied similarly to what we have done here. Actually, the result can be obtained straightforwardly if the analogy of capillary deformation with 2D electrostatics is employed (see, e.g., Ref. [25]): The net force creates "capillary monopoles" Qcap with an interaction energy V

torque creates "capillary dipoles action energy between them is Vmen(d) ~ ^ cap This decay is sufficiently slow to dominate asymptoti cally the nematic-mediated repulsion. However, dimen

(d) = Q2ap/(2n7)In(CA/d) (see Eq. (22)). A net

Pcap and the inter-

(|Pcap|/d)2

sional analysis yields P(

so that the amplitude of this capillary energy, which must vanish for large film thicknesses h to, has also a too small numerical value to explain the experimental results reported in reference [13].

We recall that for simplicity our calculations have been carried out in the "strong anchoring" limit. "Weak anchoring" would lead to a different director field, and thus to different values of the dimensionless forces eF and en (see Eq. (9)). It would also lead to an explicit contribution to the expression (10) for the free energy accounting for the anchoring "wetting" energies (Eq. (2)), which in the "strong anchoring" limit considered so far were subsumed in a (quantitatively negligible) additive renormalization of the surface tensions. The consequences of this more general case of no strong anchoring have to be explored yet, but we do not expect that this alters our conclusion that a logarithmically varying capillary attraction is ruled out by mechanical isolation (eF = en), because this latter result is based on very general principles.

Finally we note that for a two-dimensional system, the interactions between two upright circular colloids trapped at the nematic-isotropic interface have been studied numerically [39] using a Landau-de Gennes free energy. In three dimensions, this corresponds to long, cylindrical colloids which are aligned side-by-side at the interface. For boundary conditions that yield a similar, quadrupolar behavior of the director field around a single colloid [40], the effective interaction is found to be repulsive and consistent with a power law decay (for intermediate distance d/R up to 7). In this particular two-dimensional situation, the repulsive interactions appear to be longer-ranged « d-1 and in the numerical results for the director field there is no trace of a sizeable interface deformation which would lead to capillary attractions.

In summary, we have presented several arguments which rule out an asymptotic attraction of capillary origin in the effective interaction potential between two colloids at a nematic interface. If one discards effective pair potentials as the source of the stability of colloidal mesostruc-tures as the ones found in reference [13], the question arises whether this stability is a consequence of genuine

many-body effects [41]. Nematic surfaces (without colloids) can stabilize regular patterns of surface defects [30]. If these defects persist in the presence of colloids, the experimentally observed regular order of the colloids might be attributed to them. This issue calls for further experimentaI and theoreticaI investigations.

We acknowledge useful discussions with S. Chernyshuk, O. Lavrentovich, B. Lev, P. Patricio, and J.M. Romero-Enrique. M. O. thanks the German Science Foundation (DFG) for financial support through the Collaborative Research Centre (SFB-TR6) "Colloids in External Fields", project section D6-NWG. A. D. acknowledges financial support from the Junta de Andalucía (Spain) through the program "Retorno de Investigadores".

Appendix A. Corrections to the strong-anchoring limit

Here we consider the asymptotic director field around a single colloid at the nematic-air interface in the case of finite anchoring strength. For large distances from the colloid, the overall deviations from the preferred director n = (n1,n2 ,n3) = (0, 0,1) are small and the total free energy of the nematic phase is given by the harmonic approximation for the bulk part (Eq. (12)) and the Poulini expression for the surface part:

Fne = K I d3r E (Vni)2 + Wo I • eA)

K / d3r f(Vni)2 + (Vn2)2)

2 JVne, V '

+ W1/ dA(1 - ni - n2).

The second equation holds because the normal e^ of the nematic-air interface Aair-ne is parallel to the z-direction and n2 = 1. For particles with radii R in the pm range, such as those investigated in reference [13], the dimension-less parameter a = K/(|W1|R) is smaller than 1. Minimizing the free energy renders the Laplace equation in the bulk,

= 0 (i = 1,2), (A.2)

supplemented with the Robin boundary condition: d

aR nj(x, y, z = 0) — nj(x, y, z = 0) = 0 (i = 1, 2).

The general, asymptotic solution is given by the multipole ansatz (Eq. (14))

n = qi 1 + E^+53 Q^^ +..., (A.4)

£=1 .3,7=1

with the requirement of rotational covariance around the z-axis. In view of the boundary condition (Eq. (A.3)), we

note that the derivative in equation (A.3) applied to a multipole of order m generates a multipole term of order m + 1. Thus, to leading order in 1/r, the boundary condition can only be met if the leading multipole (nf) fulfills the strong-anchoring condition nf(z = 0) =0 and is accompanied by a subleading multipole (nff) which is connected to the leading multipole through the condition aRdznf(z = 0) = nff(z = 0). Both the monopole and the dipole do not fulfill the strong-anchoring condition nf (z = 0) = 0. Therefore, as before, we find that the multipole expansion starts with the quadrupole solution and that it is necessarily accompanied by a hexapole:

ni = ni + nf,

ni = Q^ + Qma%ia3 ^,

ri(r2 — 5z2)

(A.5) (A.6)

+ aQmag Reias

ra(r2 — 5z2)

Since dz nf (z = 0) = 0, this solution (Eqs. (A.5-A.7)) exactly satisfies the boundary condition in equation (A.3). Thus the asymptotically dominant director field consists of a quadrupole-hexapole superposition. Furthermore, the magnitude of the accompanying hexapole moment is small due to the factor a, therefore the results reported before for strong anchoring (based on the leading quadrupole only) are unaffected by a finite anchoring strength at the nematic-air interface.

Appendix B. Force balance in a general configuration

For the benefit of the reader, we discuss in this Appendix some previous results [28,25] concerning the force balance of a general equilibrium configuration and demonstrate how the amplitude of an asymptotic, logarithmically varying interfacial deformation is determined solely by this mechanical condition of force balance.

Figure 4 represents a colloidal particle in equilibrium at the deformed interface; in general the deformation is not small. The condition of mechanical equilibrium implies that locally the net force on any part of the system must vanish. Thus one has:

1) Each of the fluid phases is in equilibrium. We introduce the forces exerted by each fluid phase on the particle, on the whole fluid meniscus, and on the wall of the container as

1(2) ._ part .

dA • n1(2),

Fm(e2n := — (+)/ dA • n1(2),

J Smen

F^l := — /1(2) dA • n1(2),

(B.1) (B.2) (B.3)

Fig. 4. A general configuration of the particle and the fluid interface. The arrows indicate the orientation of the correspond-

ing surfaces: Smen is the interface, Spart is the surface of the particle in contact with the lower (upper) fluid phase, Swaii is the surface of the container of the system (sketched here as a quadrangular box for simplicity). The three-phase contact line between the particle and the interface is denoted as Co.

the surfaces (see Fig. 4). The superscript 1(2) indicates the lower (upper) phase in Figure 4. The condition of mechanical equilibrium of each phase under the influence of these three forces reads

F1(2t + Fm2) + F1(2,l

part 1 men 1 wall

The total force exerted by the fluids on the particle is F

F part + Fpart, and on the meniscus it is

part 2 part

Fmen := F

men + Fmen. The expressions for these forces reduce to those given in equations (6) and (7), respectively, upon evaluation in the reference configuration depicted in Figure 2.

If the condition of mechanical equilibrium is applied locally to an infinitesimal volume in the bulk of each of the fluid phases, it turns into V • n1(2) = 0. In the nematic phase this condition yields the field equations determining the director field.

2) The particle is in mechanical equilibrium under the combined action of Fpart and the tension exerted by the interface on the particle at the three-phase contact line Co. This tension can be expressed in terms of a line integral involving the surface tension 7:

Fcontact . 7

where ec is the unit vector tangent to the interface, normal to the contact line, and oriented towards the particle side. Therefore, the condition of mechanical equilibrium reads

contact =

respectively, in terms of the stress tensor n1(2)(r) in each fluid phase with due account for the orientation of

3) Any piece Sint C Smen of the fluid interface is in mechanical equilibrium. The force on Sint exerted by the

fluid phases and the tension exerted on this piece at its boundary, dSint, are balanced:

dA • [n2 - n1

d^ ec = 0,

FWall + FWall - 1

d^ ec = 0,

7V2u(r) = n1z(r) - n2z(r) (r G Sext)

i(p, = Ao + Bo ln

cos —

f d^'dp'p'[niz(r')— nz2z(r')]ln

2nY ■/Sext

|r — r'i

(B.10)

We can apply equation (B.7) to the piece Sint = Smen\Sext enclosed by the contact line C0 and the circle p = £. Invoking equation (B.6), one has

with ec oriented towards the exterior of Sint. If this expression is applied locally to an infinitesimal piece of interface, it turns into an equation for the interfacial deformation relating the mean curvature of the interface to the pressure jump accross it. If, in addition, the interface deviates only slightly from a flat interface (identified with the plane z = 0), this equation reduces in turn to the well-known equation (B.9) below for the local height w(x, y) over the plane z = 0.

The net force balance of the whole system follows from the three separate balance conditions (B.4, B.6), and (B.7): with S^t = Smen (so that dSint = Co U CWaii; Cwaii is the three-phase contact line between phase 1, phase 2, and the container enclosing the system), one finds that, as expected, in equilibrium the net force of the outer environment on the system must vanish7:

Fpart + / dA • [n2 — n1

J Sint

f ^ du ~ —Y eJ d0 p

(B.11)

where the last, approximate equality involves the leadingorder term in Vw of the line integral. (This approximation is justified because the interface deviates only slightly from a flat interface at the circle p = £.) Evaluating this latter term for the general solution (B.10), one finally finds that the amplitude of the logarithmic term in equation (B.10) is proportional to the force exerted by the upper and the lower fluid on the particle and on the meniscus:

Bo = --ez • [Fpart + Fmen] • (B.12)

In obtaining equation (B.12) we have used that in the region Sext (where the interface is almost flat) one has to leading order

where ec points to the exterior of the system.

The condition of mechanical equilibrium can be applied advantageously to obtain a precise statement about the amplitude of an interfacial deformation w(p, varying logarithmically with the lateral distance p ^ R from the particle with radius R. Far away from the particle, interface deformations and their gradients are small, so that the linearized equation holds. Thus there is a distance £ beyond which the linear theory is applicable:

dA • [n2 — n1] • ez

where the piece of interface Sext is enclosed by the circle p = £ and the line Cwan. The general solution to this inhomogeneous Laplace equation can be written as

/ d^ dpp[ffz2*(r) - nL(r)].

" Sext J Sext

(B.13)

In conclusion, if the stress ffZz(r) — nZz(r) decays sufficiently fast with the distance p from the particle, the solution in equation (B.10) will be dominated asymptotically by the logarithm with an amplitude given by equation (B.12), provided this amplitude does not vanish8. As one can infer from equation (B.4), this latter condition means physically that the walls exert a nonvanishing force: a logarithmic term can only arise if mechanical isolation of the system "colloid + interface" is violated. (Note that Eq. (B.12) follows actually from Eq. (B.8) if the interfacial deformations are small.)

We emphasize the generality of this result: it only requires that the interface departs slightly from a flat one for distances sufficiently far from the particle. (Otherwise it is actually not useful to speak of a logarithmically varying deformation to begin with.) In particular, the interfacial deformation close to the particle may be arbitrarily large, the particle itself need not be perfectly spherical, or it may even consist of a many-body configuration lacking any kind of symmetry.

The fixed values of the constants Am, Bm, and are determined by the boundary conditions. This expression reduces to equation (11) in the particular case of rotational symmetry and a wall located at infinity.

The reasoning can be easily generalized to the case that in addition to the short-ranged influence of the wall there are also external fields (gravity, electric force) acting on any part of the system.

8 To which extent this asymptotic regime is actually observable in a particular experimental realization depends on the precise functional form of n1(2)(r) and the values of the constant parameters Am, Bm, and entering the solution given by equation (B.10).

References

1. P. Pieranski, Phys. Rev. Lett. 45, 569 (1980).

2. R. Aveyard, B.P. Binks, J.H. Clint, P.D.I. Fletcher, T.S. Horozov, B. Neumann, V.N. Paunov, J. Annesley, S.W. Botchway, D. Nees, A.W. Parker, A.D. Ward, A.N. Burgess, Phys. Rev. Lett. 88, 246102 (2002).

3. B.J. Park, J.P. Pantina, E. Furst, M. Oettel, S. Reynaert, J. Vermant, Langmuir 24, 1686 (2008).

4. C. Zeng, H. Bissig, A.D. Dinsmore, Solid State Commun. 139, 547 (2006).

5. F. Bresme, M. Oettel, J. Phys.: Condens. Matter 19, 413101 (2007).

6. M. Oettel, S. Dietrich, Langmuir 24, 1425 (2008).

7. F. Ghezzi, J.C. Earnshaw, J. Phys.: Condens. Matter 9, L517 (1997).

8. J. Ruiz-García, B.I. Ivlev, Mol. Phys. 95, 371 (1998).

9. M. Quesada-Perez, A. Moncho-Jordá, F. Martinez-Lopez, R. Hidalgo-Alvarez, J. Chem. Phys. 115, 10897 (2001).

10. O. Gomez-Guzman, J. Ruiz-García, J. Colloid Interface Sci. 291, 1 (2005).

11. W. Chen, S.S. Tan, Z.S. Huang, T.K. Ng, W.T. Ford, P. Tong, Phys. Rev. E 74, 021406 (2006).

12. J.C. Fernández-Toledano, A. Moncho-Jordá, F. Martinez-Lopez, R. Hidalgo-Alvarez, Langmuir 20, 6977 (2004).

13. I.I. Smalyukh, S. Chernyshuk, B.I. Lev, A.B. Nych, U. Ognysta, V.G. Nazarenko, O.D. Lavrentovich, Phys. Rev. Lett. 93, 117801 (2004).

14. A.B. Nych, U.M. Ognysta, V.M. Pergamenshchik, B.I. Lev, V.G. Nazarenko, I. Musevic, M. Skarabot, O.D. Lavrentovich, Phys. Rev. Lett. 98, 057801 (2007).

15. P.A. Kralchevsky, K. Nagayama, Adv. Colloid Interface Sci. 85, 145 (2000).

16. M.G. Nikolaides, A.R. Bausch, M.F. Hsu, A.D. Dins-more, M.P. Brenner, C. Gay, D.A. Weitz, Nature 420, 299 (2002).

17. M. Megens, J. Aizenberg, Nature 424, 1014 (2003).

18. M.G. Nikolaides, A.R. Bausch, M.F. Hsu, A.D. Dinsmore, M.P. Brenner, C. Gay, D.A. Weitz, Nature 424, 1014 (2003).

19. L. Foret, A. Wurger, Phys. Rev. Lett. 92, 058302 (2004).

20. K.D. Danov, P.A. Kralchevsky, M.P. Boneva, Langmuir 20, 6139 (2004).

21. M. Oettel, A. Domínguez, S. Dietrich, Phys. Rev. E 71, 051401 (2005).

22. M. Oettel, A. Dománguez, S. Dietrich, Langmuir 22, 846 (2006).

23. K.D. Danov, P.A. Kralchevsky, Langmuir 22, 848 (2006).

24. K.D. Danov, P.A. Kralchevsky, M.P. Boneva, Langmuir 22, 2653 (2006).

25. A. Domínguez, M. Oettel, S. Dietrich, J. Chem. Phys. 128, 114904 (2008).

26. M. Oettel, A. Domínguez, S. Dietrich, J. Phys.: Condens. Matter 17, L337 (2005).

27. A. Wurger, L. Foret, J. Phys. Chem. B 109, 16435 (2005).

28. A. Domínguez, M. Oettel, S. Dietrich, J. Phys.: Condens. Matter 17, S3387 (2005).

29. A. Domínguez, M. Oettel, S. Dietrich, J. Chem. Phys. 127, 204706 (2007).

30. P.G. de Gennes, The Physics of Liquid Crystals (Clarendon, Oxford, 1993).

31. M. Kleman, O.D. Lavrentovich, Soft Matter Physics: An Introduction (Springer, New York, 2003).

32. J. Schwinger, L.L. DeRaad jr., K.A. Milton, W. Tsai, Classical Electrodynamics (Perseus, Reading, 1998) Chapt. 11.2.

33. P. Poulin, H. Stark, T.C. Lubensky, D.A. Weitz, Science 275, 1770 (1997).

34. T.C. Lubensky, D. Pettey, N. Currier, H. Stark, Phys. Rev. E 57, 610 (1998).

35. R.W. Ruhwandl, E.M. Terentjev, Phys. Rev. E 55, 2958

(1997).

36. S. Ramaswamy, R. Nityananda, V.A. Raghunathan, J. Prost, Mol. Cryst. Liq. Cryst. 288, 175 (1996).

37. J. Fukuda, B.I. Lev, K.M. Aoki, H. Yokoyama, Phys. Rev. E 66, 051711 (2002).

38. H. Stark, Eur. Phys. J. B 10, 311 (1999).

39. D. Andrienko, M. Tasinkevych, S. Dietrich, Europhys. Lett. 70, 95 (2005).

40. D. Pettey, T.C. Lubensky, D.R. Link, Liq. Cryst. 25, 579

(1998).

41. O. Lavrentovich, private communication.