Dislocation effects in smectic-A liquid crystals* P. S. Pershan Division of Engineering and Applied Physics. Harvard University. Cambridge. Massachusetts 02138 (Received 5 November 1973) A method for calculating stress-strain fields around edge dislocations in smectic-A samples is discussed. In large part the method is isomorphic with the formalism for calculating magnetic fields around lines of electric current. The force law between dislocations that follows from the analogy is equivalent to the accepted force law between dislocations in crystals. In addition to rederiving the expression for the strain field surrounding an isolated edge dislocation that was first given by de Gennes, we present the solutions for the stress-strain fields surrounding dislocations near one or two boundaries and also the strain field surrounding an edge dislocation that is curved to form a circular loop. The stress-strain fields surrounding other defects with the same symmetries and boundary conditions can be expressed in terms of the above-mentioned solutions using Green's function techniques. The relative stability of dislocations in samples with different types of boundaries and also the effects of dislocations on the elastic properties of smectic samples are also discussed in some detail. We comment briefly on the relation between the analogy discussed here and an earlier one developed by de Gennes. I. INTRODUCTION Recent interest in the physical properties of liquid crystals, particularly of the smectic type, warrants the same type of detailed understanding of the effects of dislocations on these systems as is currently available for conventional crystals.l,a Although there have been studies directed toward this end, see for example the work of Friedel and Kleman 3 and Bouligand,4 there remain further questions concerning the effects of dislocations in liquid crystals that we cannot yet answer. For example, in the presence of dislocations the existing hydrodynamic theories of liquid crystals are not strictly valid. 5,6 We require criteria for when dislocation effects can or cannot be neglected and the conditions under which the hydrodynamics are applicable. The analogy developed by de Gennes 7 between superfluids and smectics, together with the observation that dislocations in smectics are analogous to vortices in superfluids,a suggest the type of effects that might be expected. Although focal conic -type defects do occur in smectic samples under a variety of conditions; it does seem that single -domain planar samples can be prepared without the focal conic structure. For these samples simple topological considerations suggest the possible occur- rence of dislocations and there is some experimental evidence supporting this. 3,4 Furthermore, Bouligand4 has argued that focal conic defects are not sufficient to explain observed textures in smectics and that dislocations must be present, More recently Clark and Meyer 9 proposed dislocation effects to account for the low-frequency relaxation in recently observed light scattering phenomena.10,u Durand12 and also de Gennes 13 have discussed how defects in the smectic layers would influence light-scattering cross sections. de Gennes even presented one specific solution for the strain pattern surrounding an edge dislocation in the smectic-A structure, In this paper we will further analyze the stress-strain patterns surrounding edge dislocations in smectic-A liquid crystals subjected to specific boundary conditions and externally applied forces, Forces of one dislocation on another, and also the forces of boundaries on dislocations are calculated from the stress field using the Peach-Koehler result, 14,15 In Sec, IT we review that result as it is applicable to smectic -A liquid crystals and discuss an analogy between the smectic-dislocation problem and the interactions between magnetic fields and electric currents. In Sec. III we rederive de Gennes's result 13 for the strain field surrounding an edge dislocation and calculate the resultant stress field. In Sec, IV we consider a pair of dislocations, the mutual forces between them, and the energy of interaction between the pair. In Sec. V we discuss the interaction of edge dislocations with some simple boundaries and in Sec. VI we discuss thin smectic samples and the influence of two boundaries on the stability of dislocations, In Sec. VII we discuss wedge-shaped samples and the relative stability 04 arrays of dislocations, In Sec, vrn we present exact solutions for the stress-strain pattern surrounding a circular dislocation loop and discuss the effects of such loops on bulk elastic measurements. In Sec. IX we comment on screw dislocations, Finally in Sec. X we discuss possibilities for further theoretical work relative to the macroscopic elastic properties of smectic samples. .,.,---PI' \ ,/ \ \ / / I \ , I I II I / ~ ........ _ _ _ .... /~/> /1,1 _:3" (________ //,1 \ ~~~----~------------~\ .... - ----'-"''-...---------::::.... 152 <- - - ... 8 ....... -/. " \ i 2 &-1--", 3 ===- \ , i ....... """ _ _ _ :>t' / I II. FORMAL DEVELOPMENT FIG. 1. Examples to illustrate the sign convention of the path integrals, Eqs. (2.2) and (2.3). 1590 The smectic-A phase is such that one can define a vector field mer) with the property that Copyright © 1974 American Institute of PhysiCS 1590 Journal of Applied Physics, Vol. 45, No.4, April 1974 Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1591 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1591 J ) dr'm(r)=IiNp(r U r 2) , P(rl.r2 (2.1) If the smectic is in mechanical equilibrium and the ve- where P(ru r 2) specifies the path for the line integral between points r 1 and r 2 and liNp(r 10 r 2) is equal to the number of smectic layers traversed by P(ru r 2).16 The quantity liN p (r10 r) plays a role similar to the phase variable 2 • 7 de Gennes introduced in his analogy to the superflUld. There, and here also, it is necessary to assume that one can define a direction of increasing phase. Thus, if dr has a component in the direction of increasing phase, dr' m(r) > 0, otherwise it is negative. If P(r u r 2) is some closed path, locity of the smectic is zero, Eq. (2.22) of Ref. 5 gives at constant pressure (2.8) Equations (2.3b), (2.5b), and (2. 8) are isomorphic with the equations describing the magnetic fields surrounding lines of electric current. The vector field m is analogous to the magnetic field H, act> is analogous to (41T)-lB, and I is analogous to (41T/c)J{r), where J(r) is the free -carrier current density. This analogy is most profitable for calculating interaction effects between relatively isolated dislocations. Equation (2. 3b) could also be considered analogous to the London equation [vxJ5 =-(n e2/mc)B] if lim is associated with the supercurre nt J ~d the dislocation density 1 is associated 5 with B. The latter analogy was introduced by de Gennes. 7 The two analogies do not conflict with each other, but rather both are useful for understanding different aspects of smectics. It is well known that a system of wires carrying fixed electric currents cannot be in stable mechanical equilibrium if the only forces are the classical electromagnetic forces due to fixed currents. Nevertheless, the electromagnetic forces between the wires can be calculated from the electromagnetic part of the free energy22 using the principle of virtual work. In order for such a collection of current-carrying wires to be in mechanical equilibrium, there must be nonelectromagnetic forces to cancel the electromagnetic force density 1. dr 'm{r)=Np , p (2.2) where N p= 0, ± 1, ± 2, ..• depending on the number and type of defects enclosed in P. Figure 1 illustrates some examples. The direction of "phase" increase is fixed by m(r). Note that m(r) is locally normal to the smectic layers17 and we take the § direction to define the average direction for m(r). For path P u N Pi = + 1 while for paths P 2 and P a, N p2 =Nps =0. A path similar to P l but surrounding only B would result in an N p = - 1. The quantity N p is equal to the net number of dislocations (i. e., positive minus negative) enclosed in the path. Consider only dislocations of unit strength, we assign to each a unit vector f along the dislocation line subject to the condition that a path P enclosing only that dislocation in a "right-hand sense" relative to f obtains N p= + 1. Thus for dislocation A, fA is out of the page and for B, tB is into the page. Formally we can define a density of dislocation lines per unit area l( r) such that F=(l/c)JxB. (2.9) § drom(r)=j dA'l{r) or Vxm(r)=l{r) • (2.3a) On the basis of the above-mentioned analogy, one can define the "force density" on an array of dislocations as F=aJ.xct>. (2.3b) (2.10) In the absence of either dislocations or externally ap- plied strains m{r) = 2/ a, where a is the average equilibrium thickness of one smectic layer. Deviations of m(r) from this value define the local strain. Away from the dislocation, one may assume that the deviations are small and express lim{r) in terms of the hydrodynamic variable IJ. describing the layer displacements 18 : (2.4) At constant pressure and temperature the local energy density associated with IiIJ.{r) is given by19 liE = ct>jil11i1J. , (2.5a) (2.5b) (2.6a) (2.6b) or in terms of lim{r), 1iE=-act>llim l • Although the question of what one means precisely by the force on a dislocation does not have an obvious answer, it has been carefully considered in the case of crystalS. 23 Equation (2.10) is just the Peach-Koehler4 • 25 equation under conditions of zero pressure. 26 In order for a particular array of dislocations to be in mechanical equilibrium, there will necessarily have to be other forces (as in the magnetic problem) to stabilize the system. In the dislocation problem there can be, for example, the line tension of a dislocation core, local pinning effects due to surface irregularities or impurities. In any event, as in the magnetic problem we can consider the force [Eq. (2.10)] analogously to the force (2. 9) in the magnetic problem. III. STRESS-STRAIN PATTERN SURROUNDING AN ISOLATED DISLOCATION de Gennes's20 expression for the energy density e = tb{ilslJ.)2 + tK[(il l l + il 22 )1J.]2 can be written2l e = ta 2b(lim s)2 + tKa2{illliml + il21im2)2 and ct>3=-ablim s , Consider a dislocation of the type A as in Fig. 1 at the position (xu xa) = (0,0) and of indefinite extent in the X2 direction. Assume also that the sample is of infinite extent in the Xl and Xs directions and that boundary effects can be neglected (we will discuss boundary effects later). Equation (2. 3b) becomes ilsliml - illlimS = - Ii (x 1)1i (xs) and Eq. (2.8) [using (2.7)] becomes (3.1a) (3.1b) ct>l = + Ka{illllim l + il 12 lim 2) , ct>2 = + Ka{il221im2 + il 12 Iim l ). J. Appl. Phys., Vol. 45, No.4, April 1974 (2.7) - abilslims + Ka(ol)slim l = O. Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1592 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1592 Elimination of 15m3 from Eq. (3.1) gives 033l5m1-A2(ol)4I5m1= -5(x)ogl5(x g), (3.2) Although for any plane (xa=const), one can calculate the first moment of the 3 component of force per unit length of dislocation, that is the torque per unit length, where A==(K/b)1/2. Taking the Fourier transform and inverting gives I5m 1(x U x g) I: Xl (T3a dx1 = iAab , (3.5) . = 412 IT f f 2.!:~2 4 % q1 exp[i(q1 X1 +'qaX3)]dq1 d qs, or 5m 1(Xu Xs) = - [4 (lTA)l /2]-lxsl x s l- S/ 2exp[ - xf(4A IXsl )-1] • (3.3a) 1S in terms of the de Gennes first obtained this result angle 50 between the local normal to the smectic layers and the z direction, 60 = - a M= al5ml' The layer dila1 tion 6a/ a = - a6m s, (3.3b) The components of stress field CP1 and CPs due directly to the dislocation are obtained by substituting Eqs. (3.3) into Eqs. (2.7). Although the stress tensor itself is not uniquely defined, any choice yielding the same acceleration xa), I Substitution of Eq. (3.3b) gives a2 (Tss = ab[8( lTA)l /2]-lX11 x s l- / exp[ - X~(4A IXal )-1] , 2 (Tl3 =- Aab[8(lTW/ ]-1Xs lxal- s /2[1 - x~(2'\ IXal )-1] )-1], l /2Xll x l- S / 2[ -! + x~(4A IXal t1] (Tn == Aab(lTAt a X exp[-x~(4Alxsl)-1]. X exp[-x~(4AIXsl The quantity s" specifies the sign of the dislocation. In the figure s" == - S /3 == + 1; however, what follows will be applicable to other cases. There will be a force on dislocation a due to the stress field of dislocation (:3, f y' /3. Strictly speaking, the force will depend on the boundary conditions; however, in the examples discussed here it will be possible to replace boundaries by image dislocations and calculate boundary forces in terms of the ff,8. Referring to Eq. (2.10), with I" = - s,,2, f~'/3=as"cpf(hf,h~, ff'8=-as"cp~(hf,h:), (3.4b) (4.3a) (4.3b) ----------------------- -----------~hla, I h3a ) where cp~(hf, h:) is the stress field [Eq. (2.7)] at a due to {:3. If hi = hf - h~, the direct force of {:3 on a has components --------~~~(h" ~-------------------­ h3) -------------------- --------------~-- f 't, 8 = (ht/ Ihal )Us "s/3a2 b(lTA Ihsl )-1/2 X exp[-h~(4Alhal)-1]}, (4.4a) f:' 8= - (Ah/ Ih312)[~(2A Ihal )-, - 1]{h "sp 2b (lTA Ihal )-1/2 X exp[ - ~(4A Ihal )-1]). (4.4b) --r ------------------" 3 .. 1 If s"s/3> 0, that is, if both dislocations are of the same -,,+--------,,---------2 ._. FIG. 2. A pair of edge dislocations O! and {3 at coordinate (hf,hp, ')'=0!,{3. The re1ativ~ coordinates are h1=hr-hf, hg = h~ - hg and the unit vector 2 is into the page. J. Appl. Phys .• Vol. 45. No.4, April 1974 type, the direc't force 11 is repulsive; two dislocations of the same type will tend to displace themselves so as to avoid being directly over one another. Similarly, if SO/S8 < 0, the force is attractive and opposite dislocations would attract each other. The force fs is either attractive or repulsive depending on the relative positions of the two dislocations. Like dislocations (8,,88> 0) will tend to increase their relative hs if ~ < 2A I hsl and decrease it if ~ > 2A I hal. These effects are most easily visualized with the help of the potential Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1593 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1593 -= (aabI4kT)2(41Te)-1-10\ If X is of the order of the layer spacing, the region dominated by elastic effects can be of the order of 10- 3 cm along the vertical. Although interaction between like dislocations is predominantly repulsive, two dislocations with relative positions that fall within the shaded region interact very weakly. Repulsive forces are such as to push like dislocations into the region of weak interaction. Strictly the line with ha = 0 is a line of zero potential and it is only because of thermal effects that the two dislocations do not take up just that pOSition. This suggests that if ex;" ternal strains were to induce an excess of edge dislocations (see, for example, Fig. 3 of Ref. 7) they would tend to move to relative positions like those shown in Fig. 4 to form a "domain wall". This is identical to the model for polygonization of single crystals discussed by Nabarro. 30 We will discuss this point at greater length below. FIG. 3. Curves of constant energy of Interaction as a function of relative posltion between two like edge dislocations (s",St = + 1). The energy parameter is given In units of a 2b (4~)- . In these units kT/a ~ 3. 6 x10-.'l. l/J(hu ha) = h ",8 8a2b(X/7r I hal )1/2 exp[ - ~(4X I ha!)-1] • (4.5) The gradient of l/J, (4.6) In the case of unlike dislocations (8",SB< -1) the potential changes sign and the relative position with ha = 0 is unstable. Although entropy effects dominate that region, one might expect that dislocations so placed will diffuse away from such a relative position (in the absence of pinning or other boundary forces) until they reach the region dominated by the elastic effects. Following that they will certainly attract and (unless dislocations have some type of short-range repulsive force) annihilate one another. Boundary effects, however, can have significant effects on the stability or instability of different relative positions. These will be discussed below. yields Eqs. (4.4). The potential l/J is the excess smectic elastic energy per unit length of dislocation by virtue of the superposition of the strain fields of the two disloca~ tions.27 In the Appendix we demonstrate this for a special case. Figure 3 displays curves of constant l/J as a function of reduced variables h/X and h/ X for the example 8 ",S 8 =+1. For typical values b=1010 ergcm- a,28 a=2xl0-7 cm, the energy per cm of dislocation at a separation X along ha is w(O, X) -= 0.56 X 10-4 ergcm-1 • In order to evaluate the Significance of this energy, w(O,X)IT must be compared to the entropy per cm of dislocation. Nabarro estimates the entropy contribution to the free energy per atomic length of dislocation in a crystal as - kT. If we assume that this is approximately the correct value for the smectic, the entropy contribution to the free energy per unit length of dislocation is kTI a -= 2 x 10-7 erg cm-1 • The shaded region in Fig. 4.2 indicates where kTI a? w and where thermal effects will dominate the elastic effects. 29 The line separating the two regions will eventually curve over and meet the vertical axis similarly to the lines indicating larger potentials. The region in which elastic-interaction effects dominate is an elongated region along the ha axis. The width of the elongated region can roughly be defined by observing the curves of constant l/J become vertical when 1:,8=0. This defines a parabola Ihal=(2X)-1~. V. BOUNDARY EFFECTS A. Rigid surface Consider two dislocations, as shown in Fig. 2 except directly above one another (hf = ~B). For simplicity take h:=-h~. Equation (4.1) obtains that in the plane xa=O the strain Om l = Omf + Omf =O. From Eq. (2.4) one can see that in this configuration the plane Xs = 0 coincides with a smectic layer. Conversely, if one assumes a semi-infinite smectic sample (xs < 0) with a single dislocation a distance h below the surface (h a= - h) and the boundary condition that the upper surface (xa = 0) coin- ------------~'-------~-------- ..---:::-- ""'-------- (4.7) Inspection of Eq. (4. 5) shows that for l/J -= kTI a, I~ IIx -=(2e)-l/2(a 3b/4kT)-1O+ 2 • The "top" of th.e region is bounded by the condition l/J(O, hJ -= kT/ a or (h/X) J. Appl. Phys., Vol. 45, No.4, April 1974 FIG. 4. Domain wall In smectic-A. Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1594 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1594 cides with a smectic layer, the boundary value problem is solved by locating an image dislocation of the same type (s",ss= + 1) directly above the true dislocation, a distance h above the surface (hs=h). From Eqo (3.4) this surface also has 0'13 = O. The dislocation will experience a force per unit length (due to the image dislocation) repelling it from the surface (5.1) The potential or elastic energy for the dislocation so located is (502) This is just (2)-S/2 smaller than the energy associated with the potential l/I{O, h) of Eq. (4.5) and illustrated in Fig. 30 Following the discussion pertaining to Fig. 3, the repulsive effect of the surface will dominate thermal effects up to distances h -10- s cm. At the boundary plane (xs=O), the stress tensor 0'33(X1 ,0) is O'ss=ab[4(1TA)1/2]-lx1h-3/2exp(-xf/4Ah). tance h below a surface (xs = 0) for which O'S3 = 0, is obtained by placing an "image dislocation" of opposite type a distance h above the surface and directly over the true dislocation. Taking SB= -1 in the boundary plane 5m 1 = - [2{1TAh)1/2f1 exp(- ~/4Ah), corresponding to Integration gives p. (00, 0) - p. (- 00, 0) = S sa. The smectic layer structure described here is shown in Fig. 5 for S8= -1. The region far below the dislocation is relatively stress and strain free. For this solution to be applicable, we must justify neglect of a surface tension contribution to the energy. The elastic energy per unit length of dislocation is attractive towards the surface with magnitude equal to Il/I(h) I as given by Eq. (5.2). Surface tension effects are neglibible if this energy is large compared to the extra surface energy associated with the illustrated surface deformation. This requires the surface tension to be small compared to bA[2(21T)1/2]-1 which is of the order of 400 ergcm- 2. Since surface tensions for organic liquids are rarely as large as 100 erg cm -2, surface tension effects are negligible. Although an edge dislocation will be attracted towards a free surface, we are unable to discuss exactly what happens when the dislocation gets near the surface. Clearly there will be a step on the surface approximately one smectic layer high but the precise strain fields near the core and also in the bulk cannot be obtained from the elastic theory being presented here. The stress-strain field in the bulk can probably be described if one considers the surface step as equivalent to a pair of s = + 1 and s = -1 dislocations at some relative spacing (5h ll 5hs); however, the sign and magnitude of this spacing would depend on specific properties of the surface tension, dislocation core, and the curvature of the dislocation. For example, suppose the dislocation forms a closed circular loop of radius R on the surface corresponding to a disk of thickness a sitting on top of the bulk smectic. Assuming a core energy per unit length y for the circumference of the disk (y:> 0), the energy of the smectic would be lowered if the disk could contract with the excess material in the disk diffusing into the bulk 0 If the disk does not contract, that would mean there was, in the bulk, some stress opposing molecular diffusion. Since the chemical potential for smectic liquid crystals has been introduced and expanded as a power series 5 in the strain field, this problem can be treated in principleo It does seem, however, that the isothermal and isobaric assumptions of the present work are not sufficiently general. Realistic discussion of this type of problem is beyond the scope of the present manuscript. C. Stress-free surface (x 1 = const) On inspection of Eqs. (3.3a) and (3.3b) or (3.4b) one can readily obtain image solutions to problems with boundaries at constant Xl S\ 3ject to a variety of conditions. In all cases, however, the relative pOSition of image to dislocation will always be at hs = 0 (see Figo 2). From either Eq. (4.4) or Fig. 3, one can see there will be no direct interaction between dislocation and image. (5.3) Strictly speaking, this stress implies that there must be a torque per unit length on the boundary planes equal to JXl 0'33 {x)dx 1= Aab '" 4 x 10-4 dyn cm/ cm 0 Although such a torque is necessary if we rigidly adhere to the stated boundary conditions, the actual torque on the surface is sensitive to the number and distribution of dislocations in the sample and also to the slight deviations in the boundary conditions from the ideal. B. Free surface Two dislocations of opposite type directly above one another (Fig. 2 with s",ss=-l, ~=O, h~=-h~=h) obtain CJ'33=ab5ms=0 in the plane X3=0. Thus the solution for stress-strain fields of a single dislocation, a dis- h FIG. 5. Schematic illustration of the smectic layer configuration of a + 1 dislocation a distance h below a surface of zero stress 0"33' J. Appl. Phys., Vol. 45, No.4, April 1974 Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1595 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1595 3 • + 6[(7/J(X1-~'X3-nD-(-1)"hs) n + (7/J(X 1 - hu Xs +nD - (-l)"hs)]. (6.1b) The force per unit length on the dislocation due to the boundary is obtained from Eq. (2,10) with the identification ... CPs = - [S33(~' hs) - (7ss(O, 0)] , (6.2a) (6.2b) CP1 = - [Sls(hu hs) - (713{0, 0)] . The force on the dislocation is due only to the stresses arising from the images; the stress of the dislocation itself must be subtracted. Consider first the value of Sss for (Xl - ~)2lAD» 1. In the limiting case of this inequality, the summation over n can be approximated as an integral: Sss(Xl> X3) '" ab(8h At 1 [2(x l - h1)/D] xl"" o r st2 exp[ - (Xl - ~)2( 4At)-1 ]dt • 1-1·- - L -... ----1.1 With the substitution t = y-2, the integral is trivially evaluated: SS3(XU x 3)- b(aI2D)(xl - hi) IXI - hi 1-1 . (6.3) This result can be obtained by a more intuitive argument that emphasizes the role of boundary conditions and external forces in determining the total stress strain fields. 31 Since there is one more layer to the right of the dislocation than to the left, the strain oms must charige by an amount Om 3= - oal a2 = liD in passing the dislocation; S3S(Xl»~) - Sss(xl «hi) = balD. Assume that there are no external forces opposing f1 and the dislocation is in some type of equilibrium state. This can only occur if the elastic energy does not change with small changes in hl" From Eq. (2.6b) this will be true if (om s)2 is the same for Xl »~ as for Xl «h1 • This implies S33(00,XS )=-SS3(_00,X S), which together with Sss(oo, Xs) -Sss(- 00; Xs) = balD gives Eq. (6.3). From Eq. (6.2a) one can show directly that if the only stress fields are'due to the images depicted in Fig. 6, the force per unit length fl = O. From Eqs, (3.4) the individual contributions of each image to the summation (6. 2a) vanish because the relative coordinate Xl - hi is zero; the images are directly above (or below) the dislocation. This particular solution corresponds to the dislocation being in a state of neutral equilibrium with respect to ~ motion. The elastic energy does not depend on the ~ coordinate of the dislocation so long as the dislocation is far from the lateral boundaries at Xl = const, 1. e., so long as (~L - I ~ I )2 » AD. This is obviously a rather special case since any slight change in the sample thickness D (uniform or nonuniform) would alter the neutral equilibrium. We will later discuss the problem of a wedge-shaped sample inwhich D increases linearly with poSition along the Xl direction. Extra layers, implying dislocations, would presumably have to appear regularly wherever D increases by one layer spacing. The wedge shape together with the interactions between dislocations and their images would presumably stabilize the pattern. For the moment, however, consider the case of a uniform change in sample thickness (oD independent of Xl). To be definite, consider superposing on the strain field FIG. 6. Location of image dislocations for boundary conditions discussed in the text. The 2 direction is into the page. The sample is of dimension L in the i direction and Ihll« L. We will see in subsequent sections that the horizontal forces on dislocations, towards or away from surfaces of constant Xu are strongly influenced by the boundary conditions placed on the surfaces on constant xs. For most practical purposes the dislocations can be assumed to be far from the boundaries at constant xl" Note, however, that these lateral boundaries are sources (and sinks) of dislocations. Stresses applied to the top and bottom surfaces can cause dislocations to move into (or out of) the sample through the lateral boundaries. See, for example, Secs. VIA and vm below. VI. SINGLE DISLOCATION IN A SMECTIC SAMPLE OF UNIFORM THICKNESS 0 A. Rigid upper and lower boundaries Consider a dislocation (8 = + 1) with coordinates (hlh s) in a smectic-A sample bounded at xs=± ~D by surfaces of the type conSidered in Sec, V A (om l = (71S= OJ./. =0), These boundary conditions are satisfied by an infinite set of image dislocations of the same sign (8 = + 1) located at (h u ±nD+(-l)"hs), n=1,2,3, ••• • This is illustrated in Fig. 6. If we take the I5m l and 15(7IJ given by Eqs. (3.3) and (3,4), respectively, to be the contribution of individual dislocations (real and image) to the total strain oM I and the total stress SIJ' oM/(xu Xs) = om/(x1 -~, Xs - hs) + L; [l5m/(x1 - ~,X3 -nD -(-1)"h3) " +l5m/(x/ -hu x 3 +nD-(-1)"hs)1, (6.1a) J. Appl. Phys., Vol. 45, No.4, April 1974 Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1596 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1596 of Eq. (6.1a) a uniform dilation lim3= - a- 2 lia to obtain a total strain field liM: = liM3 + lim g. (6.4) mated by integrating Eq. (6.6) to obtain the difference between the elastic energy with a dislocation at h3 and one at hs=O: E(h s) - E(O) = i\a 2b[S(7Ti\D)1/2]-1 n=1 From Eq, (3,4) the total stress (6.5) and cP g(h1, hg) = ba-llia obtaining force fl = - blia. A slight increase in sample thickness liD = -Da- 1lia produces a uniform force tending to move the dislocation through the sample, thereby relaxing the applied stress. In the absence of localized centers that might pin them, dislocations will always move in smectics in such a way as to exclude stress from the material. Even in the case of a perfect smectic with no dislocations, application of an external strain limg could cause dislocations to move into the sample from the sides in order to reduce the net external force on the surface to zero. Constant stress ISggl > balD could induce a constant flow of dislocations through the sample boundaries. A key factor in analyzing this type of behavior is the relative probabilities for nucleating dislocation pairs (or loops) in the bulk relative to the surface. A practical consequence of the analogy between the smectic dislocation problem and the magnetic field problem (Sec. II) is that mobile dislocations cause smectics to have a "diaelasticity" that is similar to the diamagnetism of a material with mobile (free) electrons. Absolutely free dislocations would result in absolute diaelasticity analogous to superconductivity in the magnetic problem, g2 Partially free dislocations would result in a macroscopic stress-strain relationship different from the intrinsic one defined by Eq. (3.4), with an effective beff < b. Thus any attempt to measure the intrinsic elastic constants of a smectic must always be done under conditions for which dislocation motion will not occur, That is, they must either be done at frequencies above some relaxation frequency for dislocation motion or under conditions where all dislocations are permanently pinned (if that is possible?). The pinning forces must be sufficient to balance the above force with negligible dislocation motions. Assuming that the dislocation can somehow be stabilized against lihl motion, the boundary conditions alone can stabilize the dislocation against lihg motion, The force per unit length fg on a dislocation of the type shown in Fig, 6 is fs=arpu where CP1 is given by Eq. (6,2b): f3 = i\a2b[S(7Ti\)1/2]-1 t [(2n -1 + 2hs)-1/2 D +(2n-l- 2~sy1/2_2(2n_l)-1/2J (6.9) for [tL -hl ]2»i\D, Figure 7 shows the curves of E(h g) - E(O) as a function of 2h/D in units of kTI a. (Values of a, b, i\, and T are the same as were used throughout the paper.) Qualitatively the boundary effects are dominated by the two image dislocations at Xg = ± D - h g • For hs =0, E(O) - i\a 2b[S( 7Ti\D)1/2]-1 while for hs =!D, E(h s) = i\a 2b[S(7Ti\D)1/2]-1(1 - '/.h/D)-1/2, The potential barrier for the dislocation becomes comparable to kTI ~ when (6.10) where A = iD - hs is the distance between the dislocation and the nearest boundary. To obtain an upper limit for the fluctuations let ~ = a, for the typical numbers we have been assuming the right-hand side of (6. S) is of the order of 7 x 10-S(DIi\)1/2. Thus for D/i\ ~ 1()6 the thickness of the barrier is greater than A/i\ - (a S bISkT)2(27T)-1 which is roughly the same number as was discussed in Sec, IV, Fig. 3 in connection with a single boundary. On the other hand, for samples in which D/i\« lOS the barrier becomes significant for h/i\ =(SkTI a Sb)7Tl / 2(DIi\)3/2 -7X10-S(D/i\)3/2, For D=l /lm, D/i\=500 and hID 0.15 - (0.04)1/2, which is approximately equal to the estimate resulting from Eq. (6.8) with ~ a. In thin smectic samples (D «1 /lm) with rigid boundary conditions, dislocations will be largely confined to the central position, the exception being in the case where the surface is irregular. = = de Gennes'slg initial discussion of the strain field of 20 E(h;J-E(Ol kT/a 0= I J.Lm 10 J.Lm x 2:) [(nD+2h g)"S/2_(nD-2hs)-g/2]. n~l, (6.6) 100J.Lm 10 3,5 If h3 «D, f =-3i\a 2bh [4(7Ti\)l/2]-1D-5/2 2:) (2l+1)-5/2 s s 1=0 (6.7) and for a sufficiently thin sample the thermal fluctuations will be o where· ~ is the characteristic length along the dislocation, Typically, with D= 1 p.m, i\=a gives «lih/D)2) SO. 03(al ~). More generally, for larger D the effects of the boundaries on localizing the dislocation can be estiJ. Appl. Phys., Vol. 45, No.4, April 1974 FIG. 7. Elastic energy per unit length as a function of k3 for a dislocation located as shown in Fig. 6 and with boundary conditions as discussed in the text, (iL - ht)2» W. The energy is expressed in units of kT/a. See Eq. (6.9). Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1597 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1597 0' 1; ... A m~~~~~~ _ - - - L - _ v~v~0S/~/~v~v~~~~~~~~~~~~A' smectic drop on a glass plate (with rigid boundary conditions) and a free upper surface (with boundary conditions O'ss = 0) dislocations will be expelled towards the free surface. For thin samples (D$l j..I.m) the potential gradient is large and the expulsion should be rapid, For thicker samples (D;Go 100 j..I.m) the potential gradient is large only near the boundaries and there is an extensive region is which dislocation motions may be dominated by other effects. 29 In this example, as in the previous example, the boundary effects are qualtitatively dominated by the nearest image as the dislocation approaches a surface. The numerical estimates that were made in Sec. VIA on the dimensions of the regions in which the boundary potentials are Significant are applicable here also, C. Upper and lower boundaries with 0'33 =0 FIG. 8. A rigid boundary (AA') with a step of height d. The dislocation is at 0 and its image is at 0' (ha» d). an edge dislocation emphasized that the same solution is. applicable to a boundary that has a step of amplitude (in our notation). In general, surface steps of any height d are possible and the stress-strain field due to them should be the same as for the dislocation [Eqs, (3.3)] except that the amplitude Should be scaled by a factor 2d/ a. Figure S illustrates the positions of a surface step, a dislocation at 0, and an image dislocation at 0'. We assume hs» d, so that the exact position of 0' is not critical. Neglecting the upper boundary, the net component of force fs on the dislocation is ta A free smectic film of the type discussed by Clark and Meyer 33 might be described in terms of the geometry being discussed here with boundary conditions on both the upper and lower surfaces of O'ss = O. With these boundary conditions and a dislocation located at (0, h), the image solution is obtained by placing image dislocations at hs(n) = ± nD + (-l)nh with Sign Sn = (_l)n. The elastic energy is given by E(h s) -E(O) = Aa 2b[S(1TAD)1/2]-lX + (2n -1- n=l t (-1)n[(2n -1 + 2h)-1/2 D ~r1/2 -2(2n _1)-1/2] , (6.13) Equation (6.13) is plotted in Fig, 11. The curves are almost identical to those of Fig, 7 except for being in- (6,11) If 2d/ a> 2- S / 2 , the dislocation will be attracted to the step and, presumably pinned there, For a dislocation of opposite Sign this step would be repulsive. B. Rigid lower boundary, upper boundary 0'33 =0 For a dislocation (8 = + 1) located at shown in Fig. 9 with boundary conditions Om 1 = 0 on the lower surface and 0'33 = abom 3 =0 on the upper surface, the stressstrain solutions are obtained with dislocations at h(n) =±2nD+h of sign Sn=(-l)n; h'(n)=(2n-l)D-hwith sign S~=(_l)n, and h"(n):::: -(2n -l)D -h with sign S~'::: (_1)"+1. Although it is straightforward to write down expressions for the stress -strain fields analogous to those of Sec. VI A, the most interesting quantity is the difference between the elastic energy per unit length for the dislocation located at hs and at hs = 0 E(hS> -E(O) ::: Aa2b[S(1TAD)1/2]-1 2h)-1/2 xE(-1)n 2n - 1 +D' -\2n-l-D' .. r,( L t 2h)-1/2J • (6.12) Equation (6.12) is plotted versus 2h/D in units of kT/ a in Fig. 10. The elastic energy becomes increasingly negative as the dislocation moves towards the upper surface (2h/D= + 1). The figure illustrates that for a J. Appl. Phys., Vol. 45, No.4, April 1974 ... FIG. 9. Dis location at (0. ha) with boundarie s Xs = ± iD of the type described in the text. Image positions and Signs are indicated. Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1598 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1598 verted. Edge dislocations in this type of sample will be unstable and move to either surface. Effects dependent on sample thickness follow the same type of arguments given for the previous examples and the same numerical estimates apply. For practical purposes once dislocations get near the surface the comments in Sec. V C apply. near the dislocations. Equation (7.1) describes this region satisfactorily but the layer displacements obtained from it are not perfectly flat near the boundaries. This latter fault is the easier of the two to improve upon. In Sec. V A we saw that images placed symmetrically about a rigid boundary, properly oriented, can exactly cancel the nonuniform layer displacements at the boundary. We also saw in Sec. VIA that although samples of finite thickness require an infinite series of image dislocations, the size of the energy barrier near the surface was dominated by the images nearest to the surface. On the other hand, the stabilizing forces against motion in the Xl direction were determined by the local stress O"ss= - ¢s which depends on only the sample thickness and the number of smectic layers in that thickness. Thus for a number of purposes one might consider as an approximate solution to the wedge problem the superposition of Eq. (7.1) and only the first set of image dislocations on the far side of the two boundaries. Figure 12 illustrates this approximation. Taking the equilibrium poSitions for the dislocations x1(n)=(~+n)d, xs=O with d-l = a-1 0 ensures local values of 0"33 such that the force f1 on individual dislocations is zero. Image dislocations at x1(n)=(~+n)d and x S =±D(x1), with D(x1) = Oxl(n) smooth out most of the nonuniformity in the 15Mf described by Eq. (7.1). Although there will be a small residual f1 on each dislocation because the images on either side have Xs coordinates D[x1(n + 1)] > D[x1(n -1)], the effects are not important for thick samples. This image distribution can be interpreted either as an approximate solution to the wedge problem or an exact solution to a wedge sample in which the boundaries are not perfectly flat but ripple slightly to accommodate the small residual 15M1 at the boundaries. The average value of 15M1 is (2a)-10 and the residual ripple in 15M1 is of the order of (Da)-1/2. If Da» (a/0)2 =d 2, the ripple is small. Dislocation densities satisfying this inequality will be referred to as low density . In this case the solutions of Sec. VIA are approximately applicable and the dominant contribution to 0"13 = - ¢1 [Eq. (3.4)] at the site of a particular dislocation arises from the dislocation's own first image. Further, that force is only large when the dislocation approaches the surface (Fig. 7). The images stabilize the dislocation against vertical motion exactly as in Sec. VIA. VII. WEDGE-SHAPED SAMPLE A. Thick samples The wedge -shaped sample mentioned in Sec. VIA appears to be practical for the study of dislocations. Assuming rigid boundary conditions, if one starts from a flat sample that is relatively free from defects28 and squeezes one end so as to wedge the sample, edge dislocations would probably form. Assuming a wedge angle o« 1T, the number of dislocations per unit distance that are required to minimize the total strain energy is tr l =a-10, where d is the average distance between dislocations.7 The minimum elastic energy results from locating the dislocations at Xl = nd + ~d, where the origin (Xl = 0) is the vertex of the wedge angle and n = nO, nO + 1, nO + 2•.. ; the Xl coordinate of the first dislocation is x~ = (nO + ~)d. Ignoring the boundaries for the moment, the direct contribution of the dislocations to the strain 15MI' is 15MI' = - [4(1TA)1/2]-lxsl xs l- S/ 2 (7.1) . If d 2 «A I X sl, Eq. (7. 1) can be approximated by an integral 15MI''' - [4(1TA)1/2]-lxsl X S I-S/2d- 1 x foo ° exp[-(s2/4AIXs l)]ds -(Xl-Xl) .. - (4d)-1(X/ Ixsl)[ 1 + erf[(x l - X~)(4A Ixsi )-1/2]], (7.2) where erf is the error integral. For thick enough samples the inequality is satisfied near the boundaries. Far from the sample edge, where Xl-X~» (2Alxsl)1/2, 15Mf .. - (2d)- l x s lXSI-l (7. 3a) and [from Eq. (3.3b)] 15M~"O. (7.3b) E (h.) - E(O) kT /0 20 On integration of Dlll = - a15Mf the layer displacement associated with (7.3) is Il(x u xs) = (a/2d)(x s l XSI-l)Xl . 10 (7.4) Thus for a sufficiently thick sample the far-field pattern of the strain field [Eq. (7.1)] is compatibl~ with rigid wedged boundaries. In this sense Eq. (7.1) is an asymptotic solution to the wedge problem. Note that in addition to satisfying the wedge boundary conditions the 15Mu 15Ms of Eqs. (7.3) also satisfy the homogeneous differential equations [Eqs. (3.1)] for the strain field away from the dislocations. Equations (7.3) fail as an exact solution to the wedge problem because they do not describe the strain field J. Appl. Phys., Vol. 45, No.4, April 1974 -20 FIG. 10. Elastic energy per unit length as a function of a dislocation located as in Fig. 9. See Eq. (6.12). ha for Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1599 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1599 o .2 FIG. 11. Elastic energy per unit length as a function of h3 for the problem discussed in Sec. VIC; see Eq. (6.13). boundary conditions, dislocations packed this densely will interact strongly with each other. The problem then becomes a many-body one and the correct choice for the configuration of minimum energy is not obvious. For example, one possibility is that as 8 increases, causing the dislocation density d;1:::: a- 1 8 to exceed (Da)-1/2, the single dislocation wall splits into two, each of which has density (2a)-18. Further increase in 8 might further proliferate the number of walls while essentially keeping the horizontal spacing between dislocations of the order of (Da)-1/2. Delaye, Ribotta, and Durand 34 observed light-scattering effects in smectic samples under conditions where wedgelike behavior is possible. Although they interpreted their data by a different (although perfectly plausible) explanation, their principal observations were that by touching the sample on one edge they established a "periodic deformation" with spacing that is essentially given by rP ~ Da [see Eq. (3) of Ref. 34]. Although there are ambiguities as to factors of 1T (i.e., X:::: a or X::::a/21T), the order of magnitude of the low-density criterion is that for a 100-jJ.mthick sample 8;;; 20 • VIII. DISLOCATION LOOP The force 11 on a dislocation is determined by Eq. (2.10), 11 = - acp 3' where - cp 3 is the average of the 0'33 to the right and to the left of the dislocation. Assuming some equilibrium value of x 1 (n), a dislocation moving to the right by OX1 increases the average strain Oa/ a by 80x/D(x1) and since CP3::::-ab6m3::::b6a/a::::b8D[x1(n)]-1 X 6xu the restoring force is (7.5) Consider a circular dislocation loop of radius R in the plane of the smectic layers. Using a right-handed set of cylindrical coordinates (p, 8,xJ we assume the dislocation density appropriate to the loop is (8.1) where the sign of the loop is such that the extra layer is exterior to the loop. If the Fourier transform of the Cartesian components of (8.1) are expressed in terms of cylindrical coordinates appropriate to reciprocal space: The dislocation if wm be stabilized at some position x (n) 1 «6x/ iJ)2):::: kTD(x1)8(a4b)-1(a/~)« 1 . (7.6a) 1M, cp, q3) == J lj(p, 8,xJexp{- i[qp cos(8 xpdpd8dx 3 , l1(q)== -21TiRsincpJ1(qR), l2(q) == + 21TiR coscpJ1(qR) , l3(q)::::0, cp) + qsX3]} (8.2) For typical numerical values «6x/ tf)2) "" 2.5 x 10 3 1) 8(a/ ~). D(x (7.6b) Taking D(x1) ~ 1 jJ.m and 8 ~ 0.10 we obtain «6x/ d2) ~4X10-4a/~ with d~10-4 cm. Although wedges with larger angles would have larger restoring forces [Eq. (7.5)], they would also have larger dislocation densities, i. e., smaller d. The net effect is an increase in the relative fluctuations with decreasing d. de Gennes 13 estimated the optical cross section for scattering by the strain pattern surrounding an individual dislocation pinned to one boundary in terms of the Fourier transform of the optical inhomog.eneity induced by Om 1 • The result is readily extended to both real or image dislocations located near flat boundaries. The scattering amplitude from an array of dislocations is just a linear combination of these individual amplitudes. Although the location of the boundaries relative to the dislocation effects the precise dependence of the scattering amplitude on scattering wave vector K::::kseat _kine, the scattering jntensity is largely confined to the region of (K)3"" 0. 10,11,34 For small d, that is in the high-density region with d 2 $Da, the approximation under discussion is inadequate and we have not been able to construct an acceptable solution. Even if we found some way to satisfy the J. Appl. Phys., Vol. 45, No.4, April 1974 (8.3) where J 1 is the Bessel function of first order. The same methods as were used in Sec. m can be employed here: i[qx6m(q)]::::I(q) , q3 0m 3(q) (8.4a) (8.4b) + ,\2qa[6m 1 (q)coscp + 6m 2(q)sincp]:::: 0. From Eqs. (8.3) and (8.4) Om3(q):::: - 21TX 2Rq3Jl(qR)(q~ + X2 q4)-1 , (8.5) 6m 3(p, 8, x 3 ) = - (21T)-2X R 2 x J Jl(qR)(q~ +X2q4t1 exp{i[qp cos(8 x q4dqdq3dcp. cp) +qsX3]} (8.6) The integrals over q3 and cp are straightforward; 6m3:::: - ~AR f" J (qp)J1(qR) exp( - Xq21 x o O 3 1 dq • )q2 (8.7) Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1600 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1600 In place of Eq. (3. 3b) we have Ii m3= Ax S X3 I 1-1 ( p -1 -ap- . o(pm p ») (8.13) (8.8) which is equal to 35 Light-scattering cross sections for dislocation loops can be estimaged by de Gennes's method 13 using the identities Iim1(q) lims = + R(41 x 3 1t1 o~ {exp[- (p2 + R2)(4A Ixsl )-1] (8.9) or =[qrn3(q)COS -il2(q)]q;1 = 21TR cosqSJ1(qR)(q~ + A2 q4)-1 (8. 14a) Xlo[pR(2A x 31)-1]}, I where Io(z) is the modified Bessel function of zero order: ~ Io(z)=Ea Iim 1(q) (Z )2m [ml]-2=Jo(iz). 2" and similarly Iim2(q) =21TR sinq J1(qR)(~ + A2q4)-1 • S (8. 14b) In the limit that pR»2Alx 31, the asymptotic form for Io(z)- (21TZ)-1/2 exp( + z) gives, as I x31- 0, lims- [8(1TA)l/2]-1(R/ p)1/2(p - R) Ix 31- 3/ 2 Xexp[ - (R - p)2(4A IX31 )-1]. (8.10) For R <:J P and Ix31 « R 2 /A this is exactly equal to the resuit obtained in Sec. III for the strain field Iim3(XUxS) with Xl =p -R [see Eqs. (3.3)]. To obtain the strain field limp(p, away from the dislocation omp oX 3 The strain fields of a dislocation loop midway between the type of flat rigid boundaries shown in Fig. 6 will be given by an infinite set of image dislocation loops directlyover (and under) the original at X3=nD. At any point (p, B, xs) the stress 0"3S will be the sum O"33=ab n="OO t Iim S (p,x 3 -nD) , (8.15) where lims is given by either Eq. (8.7), (8.8), or (8.9) using Eq. (8.7) at X3=0: 0"33 B,x s), observe that = - ~abAR J '" Jo(qp)J1(qR) coth(Aq 2D/2)q2 dq. o (8.16) = ap . om 3 (8.11) For thin samples, (AD)l/2« R, this is approximately 0"33 <:J - abRD-1 From Eq. (8.8) and the properties of Bessel functions o o(pm p oxs \f-1 1'" J (qp)J1(qR) dq, o O (8.17) I. ») R 0 (~ ap = - 2" oR Jo Jo(qp)Jo(qR) which can be integrated to obtain 36 0"3S= -ab/D if O 0. The ratio y /b has dimensions of length squared and a reasonable guess would be (8.19) The combination of line tension and stress induces a force per unit length on the dislocation, jp= - as _yR-1, . 3 ([x(n),8x(n)] __ "-.. --1---1-11 - - - 1- I-~ , J [x(n),oj; Do (8.20) FIG. 12. An image distribution to approximate the solution to the wedge problem discussed in the text. J. Appl. Phys., Vol. 45, No.4, April 1974 where s = - 2- 1 [O"ss(p > R) + 0"33(P < R)]. In the absence of inhomogeneities that might contribute further forces, Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1601 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1601 the condition. fp = 0 implies a state of unstable equilibrium for uniform boundary conditions. For example, under conditions of constant thickness I1gg(p > R), I1gg(p < R), and CPs are constant, independent of R. If CPg = - (aRl)-ly is the equilibrium stress at some radius Ru 'iJfp = Y2 > 0, 'iJR R~Rl Rl I (8.21) implying instability. Similarly the dislocation is in unstable equilibrium under conditions of constant external force F s applied to the upper boundary, Fs = - 'lTR~l1ss(p > R) - 'lTR2[I1SS(p < R) - I1ss(p > R)], (8.22) where 7'-1 = 2va2b(RoD)-1. If one could produce homogeneous samples free of dislocations, bulk measurements to determine b would have to be done at small strain OD/ D < a/2D or at high frequency w» either 7'-1 or the reciprocal of the time to nucleate dislocations, whichever is larger. Although we have not been able to produce a theoretical estimate of v, relaxation of the applied strain by dislocation motion requires that a definite number of molecules move into the sample from the sides. There must be an upper limit on R of the order of the speed of sound, R7'-l < lOS cm sec. Alternatively, one could imagine a dislocation loop of the type shown in Fig. 13(b) being nucleated in the interior of the sample. In place of Eq. (8.23) (8.29) F s = - 'lTR~l1sg(p > R) - 'lTR 2abD- l and in place of (8.24) fp= a('lTR~)"lFs + a2bD-l[(R/Ro)2 -~] _YR"1. (8.30) where 'lTR~ is the area of the flat surfaces. From Eq. (8.18) this can also be written (8.23) From Eqs. (8.18) and (8.23), Eq. (8.20) can be written fp= -a('lTR~)-lFg +a2bD-l [(R/R o)2 -~] _yR- l In this case growth of the dislocation loop is opposed by (8.24) the line tension and in order for there to be unstable growth F g"?- 'lTRHy(Ra)-l + abD-1[~ - (R/Ro)2]} , where R is the initial value for the loop. USing (8.19) and taking R «Ro Fs"?- 'lTR~ab[R-1 + (2D)-1]. Fs"''lTR~ab/D and 'iJ{fp)/'iJR> O. Consider, however, the case of a perfect smectic sample, free of any dislocations, and with initial thickness D that is an integral multiple of aO, where aO is the equilibrium thickness of an individual layer at zero stress. According to Clark and Meyer,9 the application of a tensil force F s > 0 would induce an instability against undulations of the smectic layers. For argument's sake we assume here that this does not occur and discuss relaxation of the stress by dislocations. Whether the undulation instability occurs before dislocation effects can relax the stress depends on the times for generating the undulations in comparison with the times for dislocation nucleation and motion. In many cases the undulation effect will dominate. On the other hand, for compressional forces F s the Clark-Meyer instability does not occur while the dislocation effects are unchanged except for obvious reversals of signs and the sense of the dislocations. The tensil force induces a strain OD/ D > O. So long as OD < ~a dislocations are energetically unfavorable and Fs= 'lTR~b(OD/D). If on> ~a, the elastic energy of the system would be lowered by the introduction of an extra layer. For example, imagine that a dislocation loop like that shown in Fig. 13(a) was to nucleate at the boundary. According to Eq. (8.24) with Fs=('lTR~)[ab(2D)-1_y(aR)-1] (8.31) This is now a typical nucleation and growth process, if only those loops that spontaneously appear with radii of the order of D will be able to grow. Others that are spontaneously generated will decay. Assume that in the absence of an external stress, a loop of radius R has energy 2'/TYR '" 2'ITba3R. If this is set equal to kT, the mean radius R - kT(2'ITba3)-1 < 10- 10 cm. It appears that generation of dislocation loops by thermal fluctuations alone will not be a significant factor unless b is much smaller than the values observed by Liao et al. 28 This might occur for some materials near the smectic-to-nematic phase transition. As a practical matter is is probably very difficult to attain the ideal smectic sample described above. In all likelihood typical samples will not be perfectly uniform (8. 25a) or using Eq. (8.19) Fs'" ('lTR~ab[(2D)-1 _R-l], fp = a 2bD-l[(R/Ro)2 -1] . (8.25b) (8.26) (8.27) If we assume a mobility per unit length of dislocation v, R= va bD-l[(R/Ro)2 -1]. 2 (a) (b) Assuming an initial value of 1-R/Ro"'D/Ro«1, R/Ro'" [1 - (D/2R o)exp(t/7')][1 + (D/2R o)exp(t/7')]-1 , (8.28) J. Appl. Phys., Vol. 45, No.4, April 1974 FIG. 13. Geometry of two possible dislocation loops that could spontaneously appear and relax an applied tensile stress. (a) Dislocation loop surrounding the sample with extra layer external to the loop. (h) Dislocation loop interior to the sample with extra layer within the loop. Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1602 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1602 and dislocation loops will be pinned inside by a variety of defects. Consider, for example, a sample with a single dislocation loop of radius R trapped inside. Equation (8.22) applies except that /]33(P> R) -/]ss(p < R) are no longer necessarily given by Eq. (8.18). External forces F s will cause these loops to grow or shrink and the preceding analysis is probably applicable as a crude approximation. For the example of Fig. 13(a) we should replace the sample radii Ro by a length characteristic of the sample inhomogeneity and for the analysis applicable to Fig. 13(b) that length probably corresponds to the initial value of the dislocation loop. The most important difference, however, will be that the inhomogeneous sample will be highly strained and values of BD < ia will cause the local strain in different regions to exceed a/2D. For inhomogeneous samples the simple relation F s == bBD/ D probably fails for all values of BD. On the other hand, at high frequencies it can still be applicable. not the most general that could be written. For example, consider addition of a term like iK*a2[(01 m 1 - 02m 2)2 or with m;==a- 1oi lJ. iK* {[(all - (22)1J.]2 + (0 1m2 + 02m J2] (9.4a) (9.4b) + 4(0121J.)2}. Since 0lm2 == 02m1 away from the dislocation core, Eqs. (9.4) differ from Eqs. (2.6) only by total derivatives and physical consequences of this difference can be expressed as a surface energy. The addition of Eqs. (9.4) to Eqs. (2. 6) would be equivalent to assuming some specific form for the energy of the surface surrounding the core and that assumption mayor may not be correct. Nevertheless, whether one adds Eqs. (9.4) to the energy density or not, the equations for the strain field are baSically unchanged. It is worth noting that substitution of the strains, Eq. (9.3), into Eqs. (9.4) gives positive definite energies. Further, if Eq. (9.4a) is included in the expression for the energy denSity, a nonzero term K*a(ou IX. SCREW DISLOCATIONS Although Eqs. (2.3b) and (2.8) can be used to solve for the strain fields surrounding a screw dislocation, the formalism described earlier cannot be used to calculate the forces on or between SCrew dislocations. Equation (2.10), which described the force on an edge dislocation, is derived from the total. elastic energy of the smectic using the principle of virtual work. For a given virtual displacement of a dislocation, one assumes that the change in elastic energy is dominated by the strain energy far from the dislocation core. For this to be true the core of the dislocation must be sufficiently "stiff" that it is not changed by small changes in the stresses far from the core. The energy of interaction between the core and the far field must also be negligible compared to the elastic energy far from the core. Unfortunately, the elastic energy of a screw dislocation is dominated by either the core region or the boundary forces that maintain the dislocation. To illustrate this consider the screw dislocation depicted in Fig. 14. According to the Sign convention discussed in Sec. II, the dislocation can be described by a vector (9.1) 1== - 3B(x1)B(x 2) and the nonvanishing parts of Eqs. (2.3b) and (2.8) are 22Bm1 - 01Bm2 == B(x 1)B(x2) , V2(V om) ==0 + (22)Bm 1 (9.5) must be added to the expression for tP1 [Eq. (2.7)]. Although it is tempting to simply add this term to the energy density and proceed to calculate forces, etc., we do not believe it can be justified without a better understanding of the core. The spatial Fourier transforms of Eqs. (9.5) are m 1(q) == iq2q-2 , m 2(Q)==-iQ1q-2. (9.6) Inhomogeneities in the optical dielectric constants BE1S and BE 23 which are proportional to Bm1 and 5m 2 are expected to be large and light-scattering cross sections can be estimated from Eqs. (9.6). Since the screw dislocation is uniform in the Xs direction, the scattering wave vector K must obey the selection rule Ks == (k· cat -k1nc)s == O. (9.2a) (9.2b) where the operators V 2 and V are two dimensional. The only solution to Eqs. (9.2) consistent with the cylindrical symmetry of the screw dislocation is Bm1 == (21T)"lX2(X~ + X~)-l , (9.3) It is straightforward to demonstrate that for this solution V' Bm == 0 and the stress fields tP O. The energy = density given by Eqs. (2.6) is zero everywhere away from the dislocation core and the force exerted by one screw dislocation on another [Eq. (2.10)] is predicted to be zero. Consider now Eq. (2.6) for the elastic energy density. As we remarked earlier these expressions are J. Appl. Phys., Vol. 45, No.4, April 1974 ~; FIG. 14. Illustration of the distortion of smectic layers around a screw dislocation. According to the sign convention discussed in the text, the vector I associated with this dislocation is down. . Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1603 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1603 This selection rule was approximately satisfied for scattering from edge dislocations. In contrast for the screw dislocation 1iE13 ~ q2 and 1iE23 ~ qu while for the edge dislocation 1iE13 ~ q1 and 1ie 2S ~ q2. Optical scattering effects that have been observed do not reflect the properties expected for screw dislocations. 10,11 Finally, note that for a regular array of screw dislocations at Xl = 0, x 2 = nd, where n = 0, ± 1, ± 2, ... , the strain pattern can be calculated from Eqs. (9.5), gives a pair of equations analogous to the superconductor problem: vxv&=-(e/mc)B, V XB = (41Tn&e/ c)v&. With v& analogous to 1im and e/mc analogous to B, the macroscopic mechanical properties of a smectic with dislocations may be similar to the macroscopic electrical properties of a superconductor with vortices. In the absence of more experimental information on the actual mechanical properties of smectics, it is impossible to know whether either of the two analogies discussed here has practical application. We hope that this paper and a future one 37 will help to define the experiments that should be undertaken. The alternative analogies mentioned above aside the discussion in Sec. vm have direct application to problems of isolated focal conic regions in the smectic samples. A circular focal conic defect in an otherwise uniform ideal smectic sample has the same symmetry as the dislocation loop discUssed above and stability does not seem to be a practical difficulty. The far -field stress-strain fields arising from the focal conic can be expressed in terms of the solutions to the dislocation loop problem if the latter are regarded as Green's functions. One has only to assume boundary conditions connecting the core of the focal conic to the region outside the core in which the linear elastic theory is applicable. In this way the mutual forces between isolated focal conic defects can be estimated. ACKNOWLEDGMENTS (9.7) and for Xl »d M 2'" (21Td)-1( - Xl) 1 1: (x~ + s2)ds = - (2d)"lX1! X1 !-1. (9.8) Since M2 = - a- 02 M, the smectic layers on one side of the array (Xl> 0) make an angle 1ilJ2 = a/ d with respect to those on the other side. An array of screw dislocations describe a twist wall between two domains. It would be interesting to try and observe light scattering from such a domain wall. X. FINAL REMARKS For the most part this article has dealt with idealized dislocations and boundary effects. In real samples edge dislocations will be curved and even become screwlike as one moves along the dislocation. Dislocation loops may be intertwined, and will rarely be flat and pure edgelike . One question to pose is whether or not it will be practical to say anything quantitative about the observable macroscopic properties of real samples with large densities of randomly distributed dislocations. Although we have not been able to answer the question, the analogy, discussed in Sec, n, between the magnetic problem and the smectic problem suggests one possibility for further research, If cf> and m are analogous to Band H respectively, it would be tempting to describe the elastic response relating cf> to m by methods similar to those used to treat the magnetic response of B to H. We have already alluded to this type of analogy by speaking of the "diaelasticity" of smectics, In the magnetic problem one begins with a microscopic B field and a set of microscopic currents. The macroscopic B and H are obtained by averaging the microscopic equations over suitable volume elements, thereby eliminating the details of the microscopic currents. A similar approach might be profitable in the dislocation problem. If macroscopic stresses and strains can be obtained by suitable averaging of the microscopic quantities, it may be possible to relate these quantities by a response function which depends on only statistical averages of the dislocations. We discussed the Simplest result of this type in relation to samples in which dislocations are easily created and free to move, Alternatively the mechanical properties of smectics can be considered in terms of the de Gennes analogy, 7 From the discussion in Sec. YIn regarding circular dislocation loops, one might hope to derive an expression relating V X 1 to the average stress in the vicinity of the loop. 37 The combination of Eqs. (2,3b) and this equation J. Appl. Phys., Vol. 45, No.4, April 1974 The authors are grateful to Professor M. Ashby, Professor R, Meyer, and Professor D, Turnbull for helpful conversations relevant to this work as well as criticisms of the manuscript. APPENDIX: INTERACTION ENERGY OF A PAIR OF ISOLATED EDGE DISLOCATIONS On substitution of Eq. (4.1) into the expreSSion for the energy density (2, 6b), the interaction terms are only those terms containing products 1imf1imT, i = 1, 3. Integration over all Xu Xs should obtain Eq, (4.5), For simplicity take h~ = - h~ = th, and ~'" = hf = O. Then using Eq. (4.2b) it is straightforward to obtain the following for the "interaction term" in the energy density: iDt E = 2a2b1im~1im~ if =0 x~ > (h/2)2 , if x~< (h/2)2. (Al) For x~ > (h/2)2, eiDt = a2bs "'s B(321TA)"1[X~ _ (hI2)2]-3/2X~ xexp(- x~ !x 3 ! {2A[X~ _ (h/2)2]}-1) which integrates to EiDt=2 f"" hs/2 dx 3 j"" dxlEiDt=ta2bs",sB(A/1Th)1/2. ."" (A2) From Eq. (4.4) Eint = I/J(O, h). Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions 1604 P.S. Pershan: Dislocation effects in smectic-A liquid crystals 1604 *Work supported in part by the Joint Services Electronics Program (U.S. Army, U.S. Navy, and U.S. Air Force) under Contract No. N00014-67-A-0298-0006 and by the National Science Foundation under Grants No. GH-33576 and GH-3440l. IFor an extensive review of the field up to 1967, see F.R.N. Nabarro, Theory of Crystal Dislocations (Clarendon, Oxford, England, 1967). 2A concise review of the early literature on liquid crystals is givenbyLG. Chistyakov, Sov. Phys.-Usp. 9, 551 (1967). 3J. Friedel and M. Kleman, J. Phys. (Paris) Supp!. C4, 30, 43 (1969). 4y. Bouligand, J. Phys. (Paris) 33, 525 (1972). 5p.C. Martin, O. Parodi, and P.S. Pershan, Phys. Rev. A 6, 2401 (1972). sH. Schmidt and J. Jahnig, Ann. Phys. (New York) 71, 129 (1972). 1p.G. de Gennes, Solid State Commun. 10, 753 (1972). 8See the comment at the conclusion of Appendix C of Ref. 5 regarding the analogy between vortices in superfluids and dislocations in smectics. 9Noel A. Clark and Robert B. Meyer, Appl. Phys. Lett. 22, 493 (1973). IONoel A. Clark and P. S. Pershan, Phys. Rev. Lett. 30, 3 (1973) . IIR. Ribotta, G. Durand, and D. Litster, Solid State Commun. 12, 27 (1973). 12Georges Durand, C. R. Acad. Sci. Paris B 275, 629 (1972). 13 p . G. de Gennes, C. R. Acad. Sci. Paris B 275, 939 (1972). 14M. Peach and J.S. Koehler, Phys. Rev. 80, 436 (1950). 15This result is reviewed by Nabarro in Sec. 2.3 of Ref. l. ISImplicit here, as in all macroscopic theories of elasticity, is the fact that fp(rl' ) lorl»a, the thickness of a single layer, so that the lack precision in how one accounts for fractional parts of a layer is unimportant. I7Again we only treat macroscopic phenomena and m(r) is really determined by suitable averages over a number of layers. The theory being developed here does not describe the region "near" to the dislocation. 18See , for example, Ref. 5, Sec. III. Here we will use the letter 1.1 to specify the variable called x there. 19See Eq. (2. 19a) of Ref. 5. A canonical change of independent variables from p, v, 1.1 (or x), and S to p, v, 1.1, T will give Eq. (2. 5a) on the assumption that v = O. 20p. G. deGennes, J. Phys. Supp!. C4, 30, 65 (1969). 21Although additional terms could possibly be added to Eq. (2. 6a), they can all be converted to surface terms and need not be included here. See, for example, Sec. IX, Eqs. (9.4). J 22 L. Landau and E. Lifshitz, Electrodynamics of Continuous Media (Addison-Wesley, Reading, Mass., 1960). 23See Sec. 2.3 of Ref. 1. 24M. Peach and J.S. Koehler, Phys. Rev. 80, 436 (1950). 25J . Weertman, Philos. Mag. 11, 1217 (1965). 2sThe pressure is defined in Ref. 5 to be the partial derivative of the chemical potential with respect to mass density at constant temperature and layer spacing. The "stress tensor" has not been explicitly defined since it is somewhat arbitrary. Any choice of stress tensor (Iii that gives the same force density, gi = - (Jj(Iii is physically acceptable. 21This potential is only applicable for an extended medium in which boundary effects are neglected. 28York Liao, Noel A. Clark, and P. S. Pershan, Phys. Rev. Lett. 30, 639 (1973). 29See Sec. 11.1.1 of Ref. 1. This is only a rough approximation since it does not take into account the line tension or core energy per unit length of the dislocation. More precisely, one must compare 1/J to kT/L where ~ is a characteristic length along the dislocation. If two points on a Single dislocation line are separated by a distance less than ~, the motions of the two points are assumed to be highly correlated. In practice ~ depends on both the core energy density and geometric factors related to boundary conditions. Thus kT/a >kT/~ and the region in which thermal effects dominate is smaller than that shown. 30See Sec. 2.3.2 of Ref. l. 31This argument was first suggested to the author by Robert B. Meyer (private communication). 32Alternatively if Eq. (2.3b) is regarded as analogous to the London equation [Y' x J. = - (n.e 2/mc)Bl, one has de Gennes's (Ref. 7) analogy to superconductors. The motion of dislocations in smectics is then analogous to the motion of vortices in superfluids. See, for example, Y. B. Kim and M. J. Stephen, in Superconductivity, edited by R.D. Parks (Marcel Dekker, New York, 1969), p. 1107. 33Noel Clark and Robert B. Meyer (private communication). See, however, Robert B. Meyer and P. S. Pershan, Solid State Commun. 13, 989 (1973). 34M. Delaye, R. Ribotta, and G. Durand, Phys. Lett. A 44, 139 (1973). 35Bateman Project, Tables of Integral Tran$/orms, Vol. II (McGraw-Hill, New York, 1954), see Table 8.11, No. 23. 36Bateman Project, Tables of Integral Transforms, Vol. II (McGraw-Hill, New York, 1954), Table 8.3, No. 18. 31This idea results from discussions with J. Prost and will be discussed at greater length in a subsequent publication. J. Appl. Phys., Vol. 45, No.4, April 1974 Downloaded 07 Nov 2010 to 128.103.149.52. Redistribution subject to AIP license or copyright; see http://jap.aip.org/about/rights_and_permissions