The Construction of the Dupin Cyclides in a Smectic A Polygonal Texture
We report experimental observations with optical microscopy of the usual so-called polygonal texture. We have made simulations of the Dupin cyclides in one small region of the sample.
Many experiments on smectic A phases of liquid crystals report on the existence of macroscopic defective textures since the pioneer work of Friedel and Grandjean of 1910 . Such a liquid crystal phase put in evidence an interconnection between physics and mathematics as nicely illustrated in . Because of the nature of the singular lines in a smectic phase, two confocal lines (ellipse and hyperbola or two confocal parabolas) both located in perpendicular plans, smectic layers take the shape of Dupin cyclides in order to keep their equidistance. Notice that other defective lines, the double helices, exist also in the smectic A phase of a liquid crystal. They were first reported by Williams  and the smectic layers at their vicinity take the form of ruled helicoids . Let us focus in this paper on the FCD textures that we describe here. The cover of  shows a nice polygonal domain, notion first introduced by Friedel . This texture is highly regular and presents different properties. Figure 1 shows the same kind of texture. We can distinguish two different networks of polygons while focusing on the bottom of the substrate (Figure 1(a)) or on the top (Figure 1(b)) for a 75 μm sample thickness.
Let us focus on a part of the ellipses network, which appear inside the white frame (see Figure 1(a)). Figure 2(a) shows a zoom of one part of the experimental network of the sample and Figure 2(b) shows the simulation of the only several ellipses: the biggest one and the other six, which are in contact with the biggest one. One of the ellipses being too small, we will not take this ellipse into account in what follows. For this simulation, we first assume that the biggest ellipse is in the plane of the top of the substrate. The six other ellipses are drawing after the realization of a small rotation of different angles in the plane of the bottom of the substrate.
After measuring different parameters like semimajor and semiminor axis for the seven ellipses and the distances between apices from the experimental observations, we try to give an overview of how the smectic layers in such a sample are keeping in mind that adjacent layers coming from two different FCDs should meet continuously. We use the Image J software, which permits us to measure directly from the optical crossed polarizers picture the values of these parameters and also of the different angles between the semimajor axis of each ellipse with respect to the biggest one (called ellipse 0) measured from its physical focus (in blue in Figure 2(b)). Figure 2(b) shows also that . By using one translation and a rotation matrix of angle in the plane of Figure 2 and also a size reduction, we were able to find how each ellipse can be obtained from the biggest one (numbered 0). What is interesting in this approach is the fact that we only use the measurement of the position of the physical focus of the biggest ellipse 0 and there is no need to take into account the positions of the other physical foci for all other ellipses , which are less well defined already due to their smaller size.
Before thinking about the construction of the layers, let us begin to recall the usual description of the smectic layers at the vicinity of two different focal lines made with an ellipse of semimajor and semiminor axes and , respectively; the distance from the center of the ellipse to the focus is called and with the confocal hyperbola . The two lines and are in two perpendicular planes and their equations are written as follows: where and ; is chosen in the locations where the hyperbola physically stops in the sample. Here, we assume in fact that the maximum and minimal values of are symmetric, whereas it is usually not the case, but it does not change our discussion here. The points where the hyperbolae meet in a sample have been called “poles” by Friedel . The segment, which joins a point of the ellipse (1) and a point of the hyperbola (1), is normal to a Dupin cyclide (2) (see ). The equation of Dupin cyclides surfaces is indeed given as follows: where and correspond to the surface parameters and each surface, that is, each Dupin cyclide, corresponds to a fixed value of . With derivatives calculations and the first two fundamental forms, we find from (2) the following two different main curvature radii (see the Appendix for a detailed calculation): In the particular case where these two numbers are equal (double root), the corresponding points are called umbilics. The strengths of these umbilical defects are often encountered in nematic liquid crystals with a higher velocity for the positive one .
The Gaussian curvature (also called total curvature) and the mean curvature can be, respectively, written as Note that the case of a torus is simply obtained replacing by zero in the expression of .
Figure 3 illustrates one smectic layer for a Toroidal Focal Conic Domain (TFCD) (Figure 3(a)) and for a Focal Conic Domain (Figure 3(b)). The FCD is made by an ellipse and a hyperbola whereas the TFCD corresponds to the case where the ellipse is degenerated to a circle and the hyperbola to a straight line. The equations of the two different sets of singular lines are given in (1) for the ellipse and the confocal hyperbola and in (6) for the circle and the straight line going through the center of the circle: where varies from to .
Taking any point of the surface, two different perpendicular plans in contain the normal at to the surface. These two plans intersect the surface into two different curvature lines, which are circles for both surfaces (torus and Dupin cyclides). These two circles of radii and , respectively, are represented in green in Figure 3.
Applying (3) and (4), we obtain the different main curvatures values taking and corresponding to the chosen point , for the torus: and ,and for the Dupin cyclide: and (note that the two main curvature lines being circles, which are planar curves, the main curvature radii of these lines can be positive or negative).By the way, let us mention that the mean curvature for the Dupin cyclide is equal to 0.
From (4), we see that the Gaussian curvature at is negative for the two surfaces. Also, we can remark that one of the main curvature lines is mingled with the generator circle of the torus but the second main curvature line is not mingled with the director circle represented in blue in Figure 3(a) .
Now, in order to build our layered system for the polygonal texture, we have to keep in mind that a focal domain possesses boundaries. The hyperbola stops in two poles and therefore two different cones exist even if only one cone whose apex is on the pole has been represented in Figure 4. Several layers existing inside the area defined by this cone have been drawn and one can distinguish the three different kinds of layers (type 1 when the Dupin cyclides are singular on the ellipse, type 3 when the Dupin cyclides are singular on the hyperbola, and type 2 when the Dupin cyclides do not meet the singular lines) .
For all the ellipses visible in the polygonal texture, we have to build the layers inside the different FCDs. One additional law named the law of corresponding cones is important . The different ellipses are in contact one to each other in one point and at this point the generators of the cones associated with the two adjacent FCDs are the same. In our tiling, we have made the simulation of the different ellipses making the hypothesis that they are all in contact, and after finding a semiempirical approach consisting of two steps, rotation and translation and size reduction of the bigger ellipse and then the direct measurement of the sizes of the ellipses via the software ImageJ, we obtain the tiling of ellipses with their confocal hyperbolae of Figure 5. Because we know how each ellipse is built from ellipse 0 using a translation, a rotation, and a size reduction, we can build the Dupin cyclides for all the ellipses using exactly the same transformations but in the 3D space. The hyperbolae have to meet all in one point but the simulations do not give enough precision to be able to find a single point; they converge towards a zone larger than a point even if it is not visible with the top view of Figure 5.
Now, we represent the different FCD associated with six ellipses and their confocal hyperbolae: one of the ellipses being too small (the yellow one). This yields to Figure 6 where a tiling has been made using Mathematica 8.0. The first remark is that the quasi-continuity of the layers is realized so that we are able to go continuously from one layer to the adjacent one. To ensure the continuity of the layers, we have to write (3) for two adjacent layers: one belonging to the layers of the FCD0 (relative to ellipse 0) and the other layer belonging to the layers of the FCD5 (relative to ellipse 5). One of the two radii of curvature is positive whereas the adjacent one is negative so that we get the following:
We therefore obtain the value of sec then the value of for all the different locations where cyclides of adjacent FCDs meet. For example, if we take the upper drawn layer of ellipse 0 and the upper corresponding one of ellipse 5, we obtain a value of ; then, the layer drawn just below corresponds to ; then, for the layer drawn just below again, .
Among Friedel’s association law, two FCDs cannot recover themselves; this argument shows that we can slightly change parameters in order that this law is better obeyed for the layers above the median plane containing approximately all the ellipses. This is clearly visible in Figure 4(b) where the layers below this median plane do not join correctly whereas they meet correctly above this same plane. In particular, our simulation permits us to slightly tilt each ellipse of a small angle (less than 2°) to permit the layers to continuously meet. There are therefore locations where Friedel’s laws are no more completely obeyed. In , a model suggests a geometrical way of extending Friedel’s law to fill the entire space of the sample with portion of Dupin cyclides and spheres, joining together not only the FCDs within a polygonal domain but also between two neighboring polygonal domains. This model has been nicely illustrated in Figure 14 of  in the case of a SmA droplet suspended in a glycerin matrix with the family of FCDs and interstices filled with spherical layers. In fact, inside a spherical droplet, we can see several circles of various sizes of radii on the surface of this spherical droplet. These circles form an Apollonius spherical packing. Two cases are to be considered.(i)When the radius of the circles is a certain value called , then these regions are filled by layers taking the form of Dupin cyclides.(ii)When the radius of the circles , then these interstitial regions are filled by layers taking the form of concentric spheres. The presence of Dupin cyclides layers permits us to reduce the free energy density because of the negative Gaussian curvature instead of a positive one for spherical layers. In , the critical value of the circle size is where denotes the radius of the droplet. In our situation, if denotes the typical size of a polygonal domain, we are in good agreement with this estimation: the size microns and the can be evaluated to be close to about 10 microns.
Figure 6 represents the case of only two ellipses belonging to the polygonal domain and their corresponding cones. Figure 6(a) corresponds to the case where so that inside the green cone FCDs have been drawn, whereas Figure 6(b) corresponds to the case where so that inside the green cone part of spheres has been drawn.
The presence of such concentric spheres has been emphasized in the construction of the defects and demonstrated using topology in . Hence, two different regions exist: one made with FCD and another one made with parts of spheres for the smallest circles.
To fill entirely the space, we have therefore to consider not only Dupin cyclides inside FCDs but also part of spheres in the interstitial regions. Figure 7 constitutes therefore just a picture to imagine how the layers could be in this tiling in first approximation. We have to keep in mind that the parts of spheres for the smallest interstitial regions have not been represented in Figure 7. In particular, the layers corresponding to the smallest ellipse (ellipse 4) should be parts of spheres because the cone of apex which lies on the ellipse 4 intercepts the spheres of center in very small circles.
We have presented a picture of the smectic layers inside a sample made with a polygonal texture. Before, for instance, the addition of nanoparticles inside a smectic phase, we should know as precisely as we can how the 3D layers are arranged inside a sample. Recently, a pillar assisted Epitaxial Assembly of TFCD in a SmA liquid crystal has been reported and presents a very regular network of TFCD . In another paper, the authors imagine a construction of an indented square pyramid around each pillar (in the case of four FCDs around a circular pillar) and the corresponding hyperbolas intersect at the pyramid’s apex along with portion of concentric spheres . This kind of questions is very interesting and in the future we can think of using the polygonal texture for including nanoparticles in a texture which is less regular than a square set of FCD but nevertheless well ordered. Notice that an important and nice work in this field has been already done in .
Calculation of the two curvature radii of a Dupin cyclide.
Making the calculation of the first derivatives of relating to and , respectively, we get the following equation: These equations yield to the expression of the three coefficients existing in the expression of the first fundamental quadratic form: where , , and are written as (see the following equation) The first quadratic fundamental form reads now as Because of the nullity of , the calculation of is easier and yields to where Now, the calculation of the different three other coefficients of the second fundamental quadratic form is needed: The calculation of the three components of the vector is needed for , , and so that we obtain at first that where . Then, the second derivatives relating to and and both of them read asUsing (A.8) and (A.9), we obtain Using (A.8) and (A.11), we obtain Using (A.8) and (A.10), we obtain The second quadratic fundamental form reads now as From , the two main curvatures of our Dupin cyclides surfaces at some point are the solutions of the following equation: Here, the introduction of the different expressions of , , , and would be long and clumsy. We therefore calculate the discriminant of (A.14). When , there are two different values for the curvatures of the surface: and . We therefore obtain the following two different main curvature radii (3) and (4) expressed in the core of the paper.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
G. Friedel and F. Grandjean, “Observation geometriques sur les liquidesa coniques fo-cales,” Bulletin De La Société Française De Minéralogie, vol. 33, pp. 409–465, 1910.View at: Google Scholar
O. D. Lavrentovich, “Designing Dupin cyclides in micro and macro worlds,” Proceedings of the National Academy of Sciences, vol. 110, no. 1, pp. 5–6, 2013.View at: Google Scholar
C. E. Williams, “Helical disclination lines in smectics A,” Philosophical Magazine, vol. 32, no. 2, pp. 313–321, 1975.View at: Google Scholar
C. Meyer, Y. A. Nastyshyn, and M. Kleman, “Helical defects in smectic-A and smectic-A* phases,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 82, Article ID 031704, 2010.View at: Google Scholar
S. Elston and R. Sambles, The Optics of Thermotropic Liquid Crystals, CRC press, 1998.
G. Friedel, “The mesomorphic states of matter,” Annals of Physics, vol. 18, pp. 273–474, 1922.View at: Google Scholar
M. Kleman and O. D. Lavrentovich, “Grain boundaries and the law of corresponding cones,” European Physical Journal E, vol. 2, pp. 47–57, 2000.View at: Google Scholar
I. Dierking, M. Ravnik, E. Lark, J. Healey, G. P. Alexander, and J. M. Yeomans, “Anisotropy in the annihilation dynamics of umbilic defects in nematic liquid crystals,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 85, no. 2, Article ID 021703, 2012.View at: Publisher Site | Google Scholar
C. Meyer, L. Le Cunff, M. Belloul, and G. Foyart, “Focal conic stacking in smectic a liquid crystals: smectic flower and apollonius tiling,” Materials, vol. 2, pp. 499–513, 2009.View at: Google Scholar
O. D. Lavrentovich, “Hierarchy of defect structures in space filling byflexible smectic A layers,” Journal of Experimental and Theoretical Physics, vol. 64, pp. 984–990, 1986.View at: Google Scholar
M. Kleman and O. D. Lavrentovich, “Liquids with conics,” Liquid Crystals, vol. 36, no. 10-11, pp. 1085–1099, 2009.View at: Publisher Site | Google Scholar
J. P. Sethna and M. Kleman, “Spheric domains in smectic liquid crystals,” Physical Review A, vol. 26, no. 5, Article ID 3037, 1982.View at: Google Scholar
A. Honglawan, D. A. Beller, M. Cavallaro, R. D. Kamien, K. J. Stebe, and S. Yang, “Pillar-assisted epitaxial assembly of toric focal conic domains of smectic-A liquid crystals,” Advanced Materials, vol. 23, no. 46, pp. 5519–5523, 2011.View at: Publisher Site | Google Scholar
A. Honglawan, D. A. Beller, M. Cavallaro Jr., R. D. Kamien, and K. J. Stebe, “Topographically-induced hierarchical assembly and geometrical transformation of focal conic domain arrays in smectic liquid crystals,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, pp. 34–39, 2013.View at: Google Scholar
Y. H. Kim, D. K. Yoon, H. S. Jeong, O. D. Lavrentovich, and H.-T. Jung, “Smectic liquid crystal defects for self-assembling of building blocks and their lithographic applications,” Advanced Functional Materials, vol. 21, no. 4, pp. 610–627, 2011.View at: Publisher Site | Google Scholar
J. Lelong-Ferrand and J. M. Arnaudiès, Cours de Mathématiques, Géométrie Et Cinématique, tome 3, relation 6’, Dunod, 2nd edition, 1977.