Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Electromagnetism wikipedia , lookup
Yang–Mills theory wikipedia , lookup
Equation of state wikipedia , lookup
Aharonov–Bohm effect wikipedia , lookup
Circular dichroism wikipedia , lookup
Introduction to gauge theory wikipedia , lookup
Renormalization wikipedia , lookup
Superconductivity wikipedia , lookup
Condensed matter physics wikipedia , lookup
Vol. 120 (2011) ACTA PHYSICA POLONICA A No. 3 Statistical Theory of Biaxial Nematic and Cholesteric Phases A. Kapanowski∗ Marian Smoluchowski Institute of Physics, Jagiellonian University, W.S. Reymonta 4, PL-30-059 Kraków, Poland (Received May 24, 2011) A statistical theory of biaxial nematic and cholesteric phases is presented. It is derived in the thermodynamic limit at small density and small distortions. The considered phases are composed of short rigid biaxial or less symmetric molecules. The expressions for various macroscopic parameters involve the one-particle distribution function and the potential energy of two-body short-range interactions. The cholesteric phase is regarded as a distorted form of the nematic phase. Exemplary calculations are shown for several systems of molecules interacting via Corner-type potential based on the Lennard–Jones 12-6 functional dependence. The temperature dependence of the order parameters, the elastic constants, the phase twists, the dielectric susceptibilities, and the flexoelectric coefficients are obtained. The flow properties and defects in nematics are also discussed. PACS: 61.30.Cz, 77.84.Nh 1. Introduction The liquid crystalline state of matter is found between solids and liquids. There are numerous different types of liquid crystalline phases, but here we will focus on the nematic state. Nematic liquid crystals are characterized by anisotropic positional short-range order and orientational long-range order [1]. In the biaxial nematic (NB ) phase the mesogenic molecules tend to be parallel to selected axes, labelled by the unit orthogonal vectors (L, M , N ). In the uniaxial nematic (NU ) phase only the N vector survives. Generally, two uniaxial phases are distinguished: prolate (NU + ) and oblate (NU − ), formed by rod-like and disk-like molecules, respectively. The uniaxial nematic phase has the D∞h point group symmetry and any second-rank tensorial property has two different principal components. The biaxial nematic phase has the D2h point group symmetry and the corresponding second-rank tensorial property has three different principal components. Most nematogenic molecules are intrinsically biaxial but they give rise to uniaxial phases as a consequence of the rotational disorder around the long molecular axis. However, this rotational disorder can be hampered by interactions favouring the molecules to stick parallel to one another. The first idea to assess the macroscopic effects of molecular biaxiality was of Freiser [2, 3]. Later, liquid crystals phases formed by biaxial molecules have been studied using molecular field treatments [4–9], the Landau–de Gennes theory [10–13], computer simulations of lattice models [14–17], an SU(3) representation [18, 19]. It was shown that single-component models consisting of biaxial molecules and interacting by properly chosen continuous potentials, can produce a biaxial phase. The transition between biaxial and uniaxial ∗ e-mail: [email protected] phases is usually second order, and a direct transition between isotropic and biaxial nematic phases is also predicted. Similar results were obtained from the analytical study of single-component systems consisting of biaxial molecules interacting via hard-core potentials [20–25], and confirmed by simulation results [26–28]. On the experimental side, a biaxial nematic phase was discovered in a lyotropic system in 1980 [29]. In 1991 phase biaxiality was reported for side-chain polymers [30]. In 1992 biaxial defects were observed in a biaxial thermotropic polymer [31]. The existence of stable biaxial phases in thermotropic low molecular weight systems was not certain for many years [32–34]. The first convincing thermotropic biaxial nematic phase was jointly reported in 2004 by Kumar, Samulski and coworkers [35–37] from bent-core mesogens based on a mesogenic oxadiazole core with lateral substituents. In the same year 2004, a novel class of biaxial phases bases on organosiloxane tetrapodes was found by Mehl and coworkers [38–40]. We would like to note that bent-core molecules have C2v symmetry and their phases and phase transitions were studied by Lubensky and Radzihovsky [41]. Among many possible phases, the orientationally ordered but optically isotropic tetrahedral nematic T phase is described (Td symmetry). Tetrahedral symmetry in nematogenic molecules was studied by Fel [42, 43] and a third-rank order parameter was introduced. Recently it was shown within the extended Straley model that the following transition sequence is possible on decreasing temperature: the second-order transition from the isotropic to the tetrahedral phase and the first-order transition from the tetrahedral to the biaxial nematic phase [44]. A detailed symmetry classification of unconventional nematic phases has been carried out by Mettout [45]. Structures and properties of liquid crystalline phases formed by bent-core molecules were reviewed by Takezoe and Takanishi [46]. (351) 352 A. Kapanowski There is another system in which a biaxial nematic phase could occur: a mixture of rod-like and plate-like molecules of comparable size and in comparable amounts [47]. Unfortunately, it has been found by theory and experiment that these mixtures tend to phase separate [48–50]. In order to avoid demixing the rods and discs must have an attractive potential, such that they would be attracted more to one another than to each other. Alternatively, the rod-like and the disk-like components should be linked together covalently, either at the terminal end or in a lateral position of the rod [51]. Recent synthetic [52] and computer simulation [53] results have renewed optimism for the achievement of the first biaxial rod-plate thermotropic nematic mixture. Biaxial materials could offer significant advantages over conventional nematic liquid crystals in display applications by leading to improved response times and to better viewing characteristics. First confirmations of the speed advantage of biaxial nematics have been obtained both experimentally [54] and from computer simulations [55]. However, a number of fundamental material-related issues have to be successfully addressed before biaxial nematics can be applied in competitive display applications: engineering of room-temperature low mass materials, controlling of optical axes alignment, selective addressability of optical axes by applied electric field, and many others [56]. In this paper we present a statistical theory of biaxial nematic and cholesteric phases. The cholesteric phase is regarded as a distorted form of the nematic phase. In Sect. 2 the description of the orientational order is presented. In Sect. 3 different potentials used for biaxial systems are described. In Sect. 4 a statistical theory of the nematic phase is showed. In Sect. 5 static distortions of the nematic phase are analysed. Sections 6 and 7 deal with electric field effects, namely dielectric anisotropy and the flexoelectric effect. Section 8 is devoted to chiral biaxial nematics. In Sect. 9 viscous properties and defects of the nematic phases are briefly described. 2. Orientational order The long-range orientational order is the defining characteristic of a liquid crystal. It is responsible for the anisotropy in the properties of the nematic phase. It is important to understand the orientational order and tools describing it. Here we introduce the definitions of the orientational distribution function and order parameters. 2.1. Orientational distribution function Let us assume that molecules are rigid blocks of D2h symmetry. The orientation of a molecule is described by the three Euler angles R = (φ, θ, ψ) or by the three orthonormal vectors (l, m, n). In a space-fixed reference frame (ex , ey , ez ) they can be written as l = lα eα = R1α eα , n = nα eα = R3α eα . m = mα eα = R2α eα , (1) The most complete description of the orientational order in the homogeneous phase is provided by the orientational distribution function, f (R), which gives the probability density of finding a molecule at a given orientation. By definition the distribution function is normalized, that is, Z Z dRf (R) = dφdθ sin θ dψf (φ, θ, ψ) = 1 . (2) The distribution function can be expanded in a complete basis set of functions spanning the relevant orientational space X (j) (j) f (R) = fµν Fµν (R) , (3) jµν (j) Fµν where are defined in Ref. [57]. The definition and (j) the main properties of the functions Fµν are recalled in Appendix A. Four different cases can be distinguished, depending on the molecule and phase symmetry. (j) (i) Uniaxial molecules in uniaxial phases. F00 (R) = Pj (cos θ) are used, j is even, Pj are the Legendre polynomials. (j) (ii) Biaxial molecules in uniaxial phases. F0ν are used, j is even. (j) (iii) Uniaxial molecules in biaxial phases. Fµ0 are used, j is even. Here the second macroscopic axis must have an external cause, an external field for example, because no molecular interaction in a system of uniaxial molecules can produce a macroscopic ordering that is less symmetric than the molecules themselves. (j) (iv) Biaxial molecules in biaxial phases. Fµν are used, j can be even or odd. (j) The expansion coefficients fµν can be identified by us(j) ing the orthogonality of the functions Fµν 2j + 1 , (4) 8π 2 where the angular brackets denote an ensemble average over the orientations of the molecules. For any function A(R), this average is defined by Z hAi = dRf (R)A(R) . (5) (j) (j) fµν = hFµν i The expansion (3) does not converge rapidly partly because of the growth of the numerical factors. An alternative expansion of the orientational distribution function has generally a better convergence X (j) (j) ln f (R) = Sµν Fµν (R) . (6) jµν In the case of the Corner potential energy of interaction the expansion (6) can be finite. We note that for non-rigid molecules the orientational distribution function should depend also on a set of torsional angles which describe the conformational state. Statistical Theory of Biaxial Nematic and Cholesteric Phases In practice the alignment is characterised not through the complete distribution function f but by some numerical parameters — order parameters. They should be zero in the less ordered phase and nonzero in the more ordered phase. In the case of the nematic phase the most impor(2) (2) tant are the averages hFµν i. The parameters hF00 i and (2) hF02 i are nonzero in both uniaxial and biaxial nematic (2) (2) phases. The parameters hF20 i and hF22 i are zero in the uniaxial phase and nonzero in the biaxial phase. In the (2) (2) limit of the perfect uniaxial order (hF00 i = 1, hF22 i = 0) (2) (2) or biaxial order (hF00 i = hF22 i = 1) the parameters (2) (2) hF02 i and hF20 i vanish. Note that in the case of the (3) tetrahedral phase the order parameter hF22 i should be (4) taken into account. The order parameters hF00 i and (4) hF44 i were computed in the Monte Carlo simulations of hard tetragonal parallelepipeds [58] in order to identify the cubatic (parquet) phase. 2.2. Order parameter tensor 3 2 1 3 αβ 1 2 Qαβ = x(Nα Nβ − δ ) − y(Lα Lβ − Mα Mβ ), (7) where L, M , and N are the orthogonal eigenvectors of Q corresponding to the eigenvalues −(x + y)/2, −(x − y)/2, and x, respectively. Three different phases can occur: isotropic (x = y = 0), uniaxial nematic (x 6= 0, y = 0), and biaxial nematic (x 6= 0, y 6= 0). The order parameter tensor can be constructed by means of the anisotropic part of the diamagnetic susceptibility or, in some cases, by means of other macroscopic response functions [1]. In the absence of electric and magnetic fields the bulk free energy for the isotropic and the nematic phases in the Landau–de Gennes theory has the form [13] £ ¤ F (Q) = F Tr(Q2 ), Tr(Q3 ) . (8) It is possible to construct order parameters in relation to specific molecular models. These order parameters will be called microscopic. For rigid biaxial molecules the order tensors can be defined [27]: ¡ ¢ 1 Qll (9) αβ = 2 3hlα lβ i − δαβ , ¡ ¢ 1 Qmm (10) αβ = 2 3hmα mβ i − δαβ , ¡ ¢ nn 1 Qαβ = 2 3hnα nβ i − δαβ . (11) The order tensors are used to determine the principal axes of the system and the order parameters in computer simulations [16]. We can express the eigenvalues of the order tensors by means of the order parameters (see Table I). On the other hand, the order parameters can be calculated from the eigenvalues of the order tensors in several ways (2) ll mm hF00 i = Qnn zz = −Qzz − Qzz , TABLE I Expansion coefficients of the order tensors eigenvalues (2) over the four order parameters hFµν i. (12) (2) Qij αβ hF00 i Qll xx Qll yy Qll zz Qmm xx Qmm yy Qmm zz Qnn xx Qnn yy Qnn zz 1/4 1/4 −1/2 1/4 1/4 −1/2 −1/2 −1/2 1 √ (2) hF02 i √ − 3/4 √ − 3/4 √ 3/2 √ 3/4 √ 3/4 √ − 3/2 0 0 0 (2) hF20 i √ − 3/4 √ 3/4 0 √ − 3/4 √ 3/4 0 √ 3/2 √ − 3/2 0 (2) hF22 i 3/4 −3/4 0 −3/4 3/4 0 0 0 0 (2) mm 3hF02 i = Qll zz − Qzz mm ll ll = Qmm xx + Qyy − Qxx − Qyy , √ The ordering of the phase can be also described by a symmetric traceless second-rank tensor Q of the form [12] 353 (13) (2) nn 3hF20 i = Qnn xx − Qyy mm ll mm = Qll yy + Qyy − Qxx − Qxx , ¢ ¡ (2) mm mm ll . hF22 i = 13 Qll xx − Qyy − Qxx + Qyy (14) (15) The procedure of assigning the eigenvalues should satisfy (2) the following conditions: (a) hF00 i > 0; (b) the same order parameters must have the same values in all the ways they are computed; (c) for each configuration at one temperature the order parameters must be as close as possible to the mean value of the order parameters of the previous temperature [16]. 3. Model potentials Statistical theories of nematic liquid crystals span between two extreme approaches. The first approach starts from Onsager’s paper [59], where the ordering isotropic– nematic transition was found by assuming that nematogenic molecules can be represented as elongated hard cylinders only interacting by excluded volume. The Onsager theory is fully entropic and athermal, the number density is a control parameter. Anisotropic, short-range, repulsive interactions are taken into account only. The theory is exact in the limit of infinitely long molecules at low concentrations. Although the theory gave rather poor predictions for the I–N transition, the statistical mechanics of simple fluids suggest that anisotropic hard-body interactions need to be understood before theories of softer interactions become predictive. A review of early computer simulations of hard-core molecules can be found in the paper by Frenkel [60]. The second approach is a mean-field description of thermotropic liquid crystals by Maier and Saupe (MS) [61]. Here the temperature is the control parameter and the theory is entirely formulated in terms of anisotropic long-range attractive interactions. The MS 354 A. Kapanowski theory is exact only when the molecular interactions are infinitely long-range. Let us present extensions of the two approaches to biaxial molecules. 3.1. Hard core interactions In hard-core interactions, the excluded volume between two molecules is proportional to the second virial coefficient in the functional expansion of the configurational free energy. That is why many efforts have been made to find the excluded volume for bodies with biaxial symmetry. Two technical approaches have been followed to achieve this. The former relies on convex-body coordinate systems [62] and takes biaxial ellipsoids as prototypes for biaxial molecules. The excluded volume between biaxial ellipsoids was first computed by Tijpto-Margo and Evans [63], and then recast in a different form. In a second approach, biaxial bodies are obtained by Minkowski addition of a sphere and another, suitably chosen body. Mulder [21] first applied this technique to compute the excluded volume between two identical spheroplatelets and his results were generalized by Taylor [64] to the case of unequal spheroplatelets. In 2005 Mulder [65] computed the excluded volume for an important class of convex bodies, the spherozonotopes. By applying Mulder’s results, Rosso and Virga [66] computed the quadrupolar projection of the excluded volume between two spherocuboids. The equation of state of gases may be expressed in the form of the virial expansion pV N N2 = 1 + B2 + B3 2 + . . . , N kB T V V (16) where p is the pressure, V is the volume, kB is the Boltzmann constant, and T is the temperature. The coefficients Bn , which are functions of the temperature, are called the virial coefficients. These coefficients have been given in terms of the interactions between molecules. In the case of hard ellipsoids it is possible to calculate B2 exactly (elliptic integrals) by the formula B2 = Re Se +Ve , where Re , Se , and Ve are the mean radius of curvature, surface and volume of the hard ellipsoid, respectively [67]. The first five virial coefficients were evaluated numerically by Vega [68] and the proposed isotropic phase equation of state was checked by means of the Monte Carlo simulations [69]. The first seven virial coefficients for isotropic fluids of hard uniaxial convex bodies were calculated in Ref. [70], where the radius of convergence was also estimated. 3.2. Anisotropic attractive interactions The extension of the MS theory for biaxial molecules was first proposed by Freiser [2, 3] and next by Straley [5]. The Straley model builds a mean-field theory upon the most general quadrupolar effective interaction compatible with the D2h symmetry X (0) (2) (2) V (R1 , R2 ) = v00 + vµν Fµν (R2−1 R1 ) , (17) µν where (2) vµν = (2) vνµ due to particle interchange symmetry. The phase diagram obtained for a fixed choice of length-to-width ratio of the platelets showed a direct isotropic-to-biaxial transition at an isolated Landau point. However, it was realized [71] that the Straley potential could lead to a different phase diagram, also exhibiting a tricritical point in the uniaxial-to-biaxial transition line and a direct isotropic-to-biaxial transition. The proper form of the potential energy is (2) (2) (2) (2) V (R1 , R2 ) = v00 F00 (R2−1 R1 ) + v22 F22 (R2−1 R1 ) , (18) and it was called the minimum coupling model capable of producing biaxial ordering [72]. Actually, there is another tricritical point in the same model, but along the line of the direct isotropic-to-biaxial transition [73]. A special case of the Straley potential was introduced by Luckhurst et al. [6] h (2) (2) V (R1 , R2 ) = −² F00 (R2−1 R1 ) + λF02 (R2−1 R1 ) i (2) (2) + λF20 (R2−1 R1 ) + λ2 F22 (R2−1 R1 ) . (19) The model was compared with the Monte Carlo simulations [14, 16]. Simulation results suggest that the potential of the form (2) (2) V (R1 , R2 ) = v00 F00 (R2−1 R1 ) i h (2) (2) (2) + v02 F02 (R2−1 R1 ) + F20 (R2−1 R1 ) (20) entails the absence of biaxial order [14]. 3.3. Gay–Berne (GB) models In 1972 Berne and Pechukas proposed simple analytical forms for the orientational dependence of the potential between two molecules which they derived from a Gaussian overlap model [74]. Next the model was modified to mimic a linear site–site potential [75]. The first generalization of a Gaussian overlap model for biaxial molecules was that of Ayton and Patey [76]. By reducing the thickness of a hard uniaxial 1:1:3 ellipsoid (lacking in a nematic phase), first a uniaxial nematic phase and then biaxial nematic phase were found. The interactions between unlike biaxial particles were described by Berardi et al. [77, 78] Φ(R1 , R2 , u) = 4²0 ²(R1 , R2 , u) £ ¤ × w12 (R1 , R2 , u) − w6 (R1 , R2 , u) , (21) where w(R1 , R2 , u) = σ0 /[u − σ(R1 , R2 , u) + σ0 ], σ(R1 , R2 , u) is the anisotropic contact distance, ²(R1 , R2 , u) is the anisotropic interaction term, σ0 and ²0 are empirical parameters. A similar generalization of the GB potential to non-homogeneous biaxial interactions is that of Cleaver et al. [79]. 3.4. Corner potential energy The Corner potential energy [80] has the form Φ12 (u/σ), where u is the distance between molecules and σ depends on the molecule orientation R1 and R2 , and on the vector ∆. For σ it is possible to write the general Statistical Theory of Biaxial Nematic and Cholesteric Phases expansion proposed by Blum and Torruella [81] which is evidently invariant under rotations and translations. It involves the 3–j Wigner symbols and the standard rotation matrix elements. In the case of the chiral biaxial molecules the lowest order terms of the expansion are £ ¤ σ = σ0 + σ1 (∆ · n1 )2 + (∆ · n2 )2 + σ2 (n1 · n2 )2 £ ¤ + σ3 (∆ · l1 )2 + (∆ · l2 )2 + σ4 (l1 · l2 )2 £ ¤ + σ5 (l1 · n2 )2 + (l2 · n1 )2 £ ¤ + σ6 (n1 · n2 )(n1 × n2 ) · ∆ £ + σ7 (l1 · n2 )(l1 × n2 ) · ∆ − (m1 · n2 )(m1 × n2 ) · ∆ + (n1 · l2 )(n1 × l2 ) · ∆ ¤ − (n1 · m2 )(n1 × m2 ) · ∆ £ + σ8 (l1 · l2 )(l1 × l2 ) · ∆ + (m1 · m2 )(m1 × m2 ) · ∆ − (l1 · m2 )(l1 × m2 ) · ∆ ¤ − (m1 · l2 )(m1 × l2 ) · ∆ . (22) The meaning of the expansion coefficients σi is explained in Table II. The potential energy of the form Φ12 (u/σ) allows us to close the temperature dependence of many physical quantities inside the functions Bs (T ) defined in Appendix B. In the case of wedge-shaped molecules or banana-shaped molecules other terms in σ are necessary (see the description of the flexoelectric effect). Let us note that σ0 determines the length scale, whereas ² determines the energy scale. It is convenient to define the kernel Z K(R1 , R2 ) = d∆(σ/σ0 )3 = XX j (j) (j) Fµν (R2−1 R1 ) , Kµν (23) µν where the following symmetries can be identified: (j) (j) (i) K(R2 , R1 ) = K(R1 , R2 ) or Kµν = Kνµ ; (ii) K(RR1 , RR2 ) = K(R1 , R2 ) for any rotation R. The molecule volume is approximated by (0) Vmol = K00 σ03 /24 . (24) The Corner potential is connected with the so-called Parsons–Lee decoupling approximation which leads to TABLE II 355 the good description of hard and soft repulsive nematogens [82, 83]. In the case of soft interactions an effective molecule volume should be taken into account due to the interpenetration of the compressible molecular core [84]. 4. Statistical theory of the phase In this section the free energy of the system, based on the statistical mechanics, is presented. The free energy depends on the one particle distribution function which describes the state of the molecule alignment. 4.1. The free energy of the phase We assume that the molecules interact via two-body short-range forces that depend on the distance between the molecules and their orientations; Φ12 gives the potential energy of interactions. The free energy for our system can be derived in the thermodynamic limit (N → ∞, V → ∞, N/V = const) from the Born–Bogoliubov– Green–Kirkwood–Yvon hierarchy [57] or as the cluster expansion [85]. The total free energy F consists of the entropy term and the interaction term, namely, F = Fent + Fint , where (25) Z βFent = βFint dr dRG(r, R)[ln(G(r, R)Λ) − 1], Z 1 dr 1 dR1 dr 2 dR2 =− 2 × G(r 1 , R1 )G(r 2 , R2 )f12 . (26) (27) Here G is the one particle distribution function with the normalization Z dr dRG(r, R) = N , (28) f12 = exp(−βΦ12 ) − 1 is the Mayer function, β = 1/(kB T ), Λ is related to the ideal gas properties. The equilibrium distribution G minimizing the free energy (25) satisfies the following equation: Z ln(G(r 1 , R1 )Λ) − dr 2 dR2 G(r 2 , R2 )f12 = const . (29) In the homogeneous phase the distribution G does not depend on the position of a molecule and we can write G(r, R) = G0 (R). It is convenient to use a dimensionless function f (R) = G0 (R)V /N with the normalization Z dRf (R) = 1 . (30) The meaning of the expansion coefficients σi . σi Description σ0 σ1 , σ2 σ3 , σ4 , σ5 σ6 σ7 , σ8 isotropic achiral uniaxial nematic achiral biaxial nematic chiral uniaxial nematic chiral biaxial nematic 4.2. Exemplary calculations Let us define a dimensionless parameter λ = B2 (T )σ03 (N/V ) . (31) (j) The set of coefficients Sµν (6) is derived from the equations 356 A. Kapanowski (j) Sµν = X (j) (j) λKρν hFµρ i for j > 0 , (32) ρ (0) hF00 i = 1 the normalization condition . (33) Fig. 1. Temperature dependence of the order parame(2) ters hFµν i. T is the dimensionless temperature. Fig. 2. Temperature dependence of the order parame(3) (4) ters hF22 i and hFµν i. T is the dimensionless temperature. We used the Lennard–Jones 12-6 potential energy, σ1 /σ0 = 0.5, σ2 /σ0 = 0.1, and σ3 /σ0 = −0.4, packing N Vmol /V = 0.1. We were looking for the orientationally stable, the lowest energy solution. On decreasing the temperature we meet the first-order transition to the uniaxial nematic phase at T = 0.374 and the second-order transition to the biaxial nematic phase at T = 0.220. The temperature dependence of the order parameters is presented in Figs. 1 and 2. Temperatures are expressed in ²/kB . Let us note that in the uniaxial nematic phase the (2) (4) most significant are hF00 i and hF00 i, whereas in the bi(2) (3) (4) axial nematic phase we should add hF22 i, hF22 i, hF22 i, (4) and hF44 i. 5. Static distortions of the nematic phase In the ideal biaxial nematic phase, the molecules are (on average) aligned along common directions (L, M , N ). In the absence of external fields the orientation of the directions is arbitrary, which is a consequence of the fact that the nematic phase arise from the spontaneous breaking of the continuous rotational symmetry. In most practical circumstances the ideal conformation will not be compatible with the constraints that are imposed by limiting surfaces of the sample and by external fields acting on the molecules [1]. There will be some deformation of the alignment and the order parameter tensor Q(r) will vary from point to point. Generally, both eigenvectors and eigenvalues of Q depend on r. But in the hydrodynamic limit the eigenvalues are nearly constant and only the directors depend on position [86]. The inhomogeneity of the director field results in a distortion (elastic) free energy. In a continuum approach this energy is obtained as an expansion about an undistorted reference state with respect to gradients of the tensor order parameter Q [87–89] or gradients of the directors [90]. There are many equivalent forms of the biaxial elastic free energy and it is possible to transform one to another after making use of relations that follow from constraints imposed by the orthonormality of the directors. The elastic constants are material constants that determine the change in the elastic free energy of a nematic phase. The phenomenological theory of the elastic constants and viscosity coefficients of biaxial liquid crystals was developed for the first time by Saupe [91]. The density-functional formalism for the elastic constants was used by Singh et al. [92]. Kini and Chandrasekhar [93] discussed the feasibility of determining some of the biaxial elastic constants by studying the effects of external magnetic and electric fields applied in different sample geometries. They showed that at certain conditions a Freedericksz transition should occur above a critical value of the field. The elastic constants influence almost all phenomena observed in liquid crystals: light scattering, defect shapes, long-wavelength fluctuations, etc. They play a crucial role in almost all applications of liquid crystalline materials, especially in display devices. On the other hand, they provide valuable information regarding the nature and importance of various anisotropies of the intermolecular potentials and of the spatial and angular correlation functions. 5.1. Distortion free energy Let us call Fd the free energy due to the distortion of the local frame (L, M , N ). A general form of its density fd (r) was derived in [90] in the case of small distortions. When a considered phase has a D2 symmetry group (the biaxial cholesteric phase) we get fd = K11 D11 + K22 D22 + K33 D33 + 12 K1111 (D11 )2 + 12 K1212 (D12 )2 + 12 K1313 (D13 )2 + 12 K2121 (D21 )2 + 21 K2222 (D22 )2 + 12 K2323 (D23 )2 + 12 K3131 (D31 )2 + 21 K3232 (D32 )2 + 12 K3333 (D33 )2 + K1122 D11 D22 Statistical Theory of Biaxial Nematic and Cholesteric Phases + K1133 D11 D33 + K2233 D22 D33 + K1221 D12 D21 + K1331 D13 D31 + K2332 D23 D32 + L123 ∂α (Lα D23 + Mα D13 ) + L231 ∂α (Mα D31 + Nα D21 ) + L312 ∂α (Nα D12 + Lα D32 ) , (34) where Kij are the chiral elastic constants, Kijkl = Kklij , Lijk = Ljik are the achiral elastic constants, D11 = Lα Mβ ∂α Nβ , D12 = Lα Nβ ∂α Lβ , D13 = Lα Lβ ∂α Mβ , D21 = Mα Mβ ∂α Nβ , D22 = Mα Nβ ∂α Lβ , D23 = Mα Lβ ∂α Mβ , D31 = Nα Mβ ∂α Nβ , D32 = Nα Nβ ∂α Lβ , The terms with Kii give 3 bulk terms, the terms with Lijk give 3 surface terms, the terms with Kijkl give 12 bulk and 3 surface terms of the form ∂α (Lβ ∂β Lα − Lα ∂β Lβ ) = 2(D23 D32 − D22 D33 ) , (36) ∂α (Mβ ∂β Mα − Mα ∂β Mβ ) = 2(D31 D13 − D11 D33 ) , (37) ∂α (Nβ ∂β Nα − Nα ∂β Nβ ) = 2(D12 D21 − D11 D22 ) . (38) The total numbers of bulk and surface terms are 15 and 6, respectively. In the case of the uniaxial cholesteric phase the distortion free-energy density has the form fd = K0 N · (∇ × N ) + 12 K1 (∇ · N )2 2 + 12 K2 [N · (∇ × N )] + 12 K3 [N × (∇ × N )] + 21 K4 ∇ · [(N · ∇)N − N (∇ · N )] (39) Therefore, we have 4 bulk (from K0 to K3 ) and 2 surface terms (K4 and K5 ). 5.2. Basic deformations of the biaxial nematic phase Three basic types of deformations appear in the continuum theory of uniaxial nematic liquid crystals: splay, twist, and bend. They extract from the distortion free energy terms with K1 , K2 , and K3 , respectively. Basic deformations of biaxial nematic liquid crystals were defined by Kapanowski [94] and only nine deformations extract one constant from the distortion free energy. These deformations are collected in Table III together with deformation parameters qi (i = 1, . . . , 9). For the sake of future reference, we will write the spatial dependence of the vectors (L, M , N ). In the case of twists, which are important for the theory of chiral phases, we have L(r) = [1, q3 z, −q2 y] , Parameter Deformation Elastic constant FE coefficient q1 M twist, N twist K1111 (K2 ) none q2 L twist, N twist K2222 (K2 ) none q3 L twist, M twist K3333 (0) none q4 N splay, L bend K1212 (K1 ) a4 (e1 ) q5 M splay, L bend K1313 (0) a5 (0) q6 N splay, M bend K2121 (K1 ) a6 (e1 ) q7 L splay, M bend K2323 (0) a7 (0) q8 M splay, N bend K3131 (K3 ) a8 (−e3 ) q9 L splay, N bend K3232 (K3 ) a9 (−e3 ) N (r) = [q2 y, −q1 x, 1] . (35) + 12 K5 ∇ · [(N · ∇)N + N (∇ · N )]. TABLE III Basic deformations of the biaxial nematic phase. The corresponding elastic constants and the FE coefficients are given, the values for the uniaxial nematic phase are in parentheses. M (r) = [−q3 z, 1, q1 x] , D33 = Nα Lβ ∂α Mβ . 2 357 (40) Splays and bends will be used in the description of the flexoelectric effect L(r) = [1, q7 y − q5 x, q9 z − q4 x] , M (r) = [q5 x − q7 y, 1, q8 z − q6 y] , N (r) = [q4 x − q9 z, q6 y − q8 z, 1] . (41) 5.3. Elastic constants In the homogeneous biaxial nematic phase composed of C2v or D2h molecules the distribution function has the form [94] G0 (R) = G0 (l · L, l · N , n · L, n · N ) . (42) In order to derive expressions for the elastic constants it is enough to assume that, in the deformed phase, the phase orientation depends on the position but the magnitude of the alignment is constant, £ ¤ G0 (r, R) = G0 l · L(r), l · N (r), n · L(r), n · N (r) . (43) Let us denote the derivatives Uα = ∂1 G0 lα + ∂3 G0 nα , Wα = ∂2 G0 lα + ∂4 G0 nα . (44) The twist elastic constants are Z 1 βK1111 = dR1 dR2 duf12 u2x W1y W2y , 2 Z 1 dR1 dR2 duf12 u2y (U1z − W1x ) βK2222 = 2 × (U2z − W2x ) , Z 1 dR1 dR2 duf12 u2z U1y U2y . βK3333 = 2 The splay-bend elastic constants have the form Z 1 βK1212 = dR1 dR2 duf12 u2x (U1z − W1x ) 2 × (U2z − W2x ) , (45) (46) (47) (48) 358 A. Kapanowski Z 1 dR1 dR2 duf12 u2x U1y U2y , 2 Z 1 βK2121 = dR1 dR2 duf12 u2y W1y W2y , 2 Z 1 dR1 dR2 duf12 u2y U1y U2y , βK2323 = 2 Z 1 βK3131 = dR1 dR2 duf12 u2z W1y W2y , 2 Z 1 βK3232 = dR1 dR2 duf12 u2z (U1z − W1x ) 2 × (U2z − W2x ) . βK1313 = (49) (50) (51) (52) (53) Other constants are connected with more complicated deformations Z 1 βK1122 = dR1 dR2 duf12 ux uy 4 £ ¤ × (U1z − W1x )W2y + W1y (U2z − W2x ) , (54) Z 1 βK1133 = dR1 dR2 duf12 ux uz 4 £ ¤ × − U1y W2y − W1y U2y , (55) Z 1 βK2233 = dR1 dR2 duf12 uy uz 4 £ ¤ × − U1y (U2z − W2x ) − (U1z − W1x )U2y , (56) K1221 = K1122 , K1331 = K1133 , K2332 = K2233 , (57) L123 = L231 = L312 = 0 . (58) In the case of the uniaxial nematic phase Uα = 0 and we get several equivalent expressions for elastic constants K1 , K2 , and K3 , K1 = K1212 = K2121 , K2 = K1111 = K2222 , K3 = K3131 = K3232 , 1 2 K4 = (K1 + K2 ), K5 = 0 . (59) (60) 5.4. Exemplary calculations Let us define η = kB T B4 (T )σ05 (N/V )2 . (61) The elastic constants can be written as finite sums of the form X [I][J] Kijkl = η hF [I] ihF [J] iAijkl , (62) [I][J] [I][J] Fig. 3. Temperature dependence of the splay-bend elastic constants expressed in ²/σ0 . T is the dimensionless temperature. where the coefficients Aijkl depend on σi only. The parameters of the model system from Sect. 4 were used to calculate the temperature dependence of the elastic constants shown in Figs. 3 and 4. The elastic constants are expressed in ²/σ0 . If we assume the molecular length σ0 = 1 nm and the interaction energy ² = 0.1 eV, then we can estimate the values of the elastic constants ²/σ0 = 16 pN and the temperature ²/kB = 1160 K. In Fig. 4. Temperature dependence of the twist elastic constants expressed in ²/σ0 . T is the dimensionless temperature. the uniaxial nematic phase three independent bulk constants are present: K1 , K2 , and K3 . In the biaxial nematic phase K1 splits into K1212 > K2121 , K2 splits into K2222 > K1111 , and K3 splits into K3232 > K3131 . The new twist elastic constant is K3333 > 0 and the new splay-bend elastic constants are K1313 > K2323 > 0. Other new constants are positive K1133 , K2233 and negative K1122 . The elastic constants in display applications provide the restoring mechanism when the applied electric field is switched off. The small biaxial elastic constants signal possible difficulties with the restoring mechanism for biaxial modes. 6. Dielectric anisotropy A static electric field imposed on a nematic is connected with at least two different processes: namely, the dielectric constant anisotropy and the flexoelectric effect. In this section we will focus on the dielectric constant and in the next section we will describe the flexoelectric effect. 6.1. Theory of the dielectric susceptibility The dielectric permittivity (²) and the susceptibility (χ = ² − 1) are tensors with three different values along Statistical Theory of Biaxial Nematic and Cholesteric Phases the different axes, at least in the most general case. The purpose of the dielectric theory is to relate them to the molecular properties, i.e., to the polarizability α and the dipole moment µ. The dielectric susceptibility is defined as: ²0 χλσ = ∂Pλ , ∂Eσ (63) where E is the electric field and P = P ind + P dip is the total electric polarization which can be divided into the part due to the molecular polarizability and the part due to the dipole moment. According to the theory of polar nematic liquid crystals by Maier and Meier (MM) [95], the susceptibility can be calculated from the equation X N 2 χλλ = hF hRiλ i(αi + µ2i βF/²0 ) , (64) V i where h and F are the factors introduced by Onsager and can be written as 3χ̃ + 3 h= , (65) 2χ̃ + 3 1 F = , (66) 1 − α̃²0 f µ ¶ 2χ̃ 3N ²0 f = , (67) 2χ̃ + 3 V αi and µi are the components of the molecule polarizability and the molecule dipole moment, respectively. In the expressions for h, F , and f the averaged susceptibility χ̃ and the polarizability α̃ were used, i.e. the anisotropy was neglected. Apart from that, the molecule under consideration was contained in a spheroidal cavity. It is unrealistic for liquid crystals and that is why the generalized Claussius–Mossotti (GCM) approach was proposed by Kapanowski and Wietecha [96, 97]. We note that Maier and Meier took into account the polarization of the surroundings by the permanent and induced dipole moment. This leads to the presence of the reaction field. In the case of small density, small dipoles or non-polar molecules, the reaction field is small and can be neglected. Let us assume that the selected molecule is in a non-spherical cavity of molecular dimensions and surrounded by the anisotropic polarized continuum. Then the internal electric field acting on the molecule depends on the external electric field, the polarization and the shape factors Ωi . Exemplary formulae for different cavities are given in Appendix B. The susceptibilities have the form P 2 (N/V ) i hRiλ i(α + µ2i β/²0 ) P 2 i χλλ = . (68) 1 − (N/V ) i hRiλ iΩi (αi + µ2i β/²0 ) We have to add that both expressions for the susceptibility (64) and (68) have some limitations. Firstly, the polarizability of a molecule is represented by an anisotropic point polarizability, whereas in reality the polarizability of a molecule will be distributed over its whole volume. Secondly, the local Lorentz field approximation is not 359 valid for dense fluids because the field from the molecules inside the cavity (the near field) cannot be neglected (the Mossotti catastrophe). In the case of isotropic polar liquids the Kirkwood–Fröhlich theory provides in principle a framework in which short ranged dipole–dipole interactions can be explicitly evaluated. The theory has been extended to liquid crystals by Bordewijk [98]. Corresponding expressions for the biaxial nematic phase can be found in Ref. [99]. 6.2. Exemplary calculations In order to calculate the susceptibilities, we have to provide data on the polarizability and the dipole moment. The ordering of the phase is described by the model from Sect. 4, where the molecules are similar to ellipsoids with three different axes 2a × 2b × 2c and a : b : c is 1 : 2 : 5. The elements of the polarizability tensor α are estimated from the model described in [96] α1 = 3Vmol , α2 = 4Vmol , α3 = 5.567Vmol . (69) The temperature dependence of the susceptibilities calculated from Eq. (64) for non-polar molecules is shown in Fig. 5. In the uniaxial nematic phase we have χzz > χyy = χxx . At the transition to the biaxial nematic phase the lower susceptibility branch splits into two independent χyy > χxx . For the considered system the reaction field is important as the susceptibilities calculated from Eq. (68) are twice smaller. Fig. 5. Temperature dependence of the susceptibilities. T is the dimensionless temperature. 7. Flexoelectric effect In a deformed uniaxial nematic liquid crystal, there should appear in many cases a spontaneous dielectric polarization described by Meyer [100]: Pα = e1 Nα ∂β Nβ + e3 Nβ ∂β Nα , (70) where e1 and e3 are the splay and the bend flexoelectric (FE) coefficients, respectively. The appearance of spontaneous polarization in liquid crystals as a result of orientational deformations is referred to as the flexoelectric effect. The inverse effect is also possible: when an 360 A. Kapanowski electric field is applied on the nematic phase, the alignment may become distorted, since a suitable distortion will imply a polarization parallel to the electric field. The FE effect can influence electro-optical properties, defect formation, and structural instability. That is why different techniques have been suggested to observe possible mechanisms producing the FE polarization: measurements of the polarisations or the surface charges induced by an imposed distortion, pyroelectric-effect-based techniques, a technique inspired by the flexoelectric-optic effect [101]. 7.1. Flexoelectric polarization in the biaxial nematic phase A statistical theory for the dipole FE polarization in the uniaxial nematic phase was derived by Kapanowski [102] in the thermodynamic limit at small distortions and a small density. The theory was applied to two systems with polar wedge-shaped and banana-shaped molecules [103]. In the case of the biaxial nematic phase the FE polarization has the form [104] Pα = s11 Lα ∂β Lβ + b11 Lβ ∂β Lα + s22 Mα ∂β Mβ + b22 Mβ ∂β Mα + s33 Nα ∂β Nβ + b33 Nβ ∂β Nα , (71) where the parameters sii and bii (i = 1, 2, 3) are not unique because if we add any constant to all of them, the polarization will not change. The physical FE coefficients ai (i = 4, . . . , 9) are connected with the proper basic deformations of the phase (see Table III). The FE polarization calculated for splay and bend deformations (41) gives Px = q7 (s11 − b22 ) + q9 (s11 − b33 ) = q7 a7 + q9 a9 , (72) Py = q5 (s22 − b11 ) + q8 (s22 − b33 ) = q5 a5 + q8 a8 , (73) Pz = q4 (s33 − b11 ) + q6 (s33 − b22 ) = q4 a4 + q6 a6 . (74) The physical coefficients satisfy the identity a4 − a5 − a6 + a7 + a8 − a9 = 0 . (75) du dR1 dR2 f12 G0 (R1 )µx (R1 )uy U2y , (80) Z a8 = du dR1 dR2 f12 G0 (R1 )µy (R1 )(−uz )W2y , (81) Z a9 = du dR1 dR2 f12 G0 (R1 )µx (R1 )uz (U2z − W2x ) . (82) In the case of the uniaxial nematic phase composed of wedge-shaped molecules (n is the long molecule axis) Uα = 0 and the Meyer two FE coefficients are recovered e1 = a4 = a6 , e3 = −a8 = −a9 . (83) The expressions (83) are valid for C2v molecules, whereas very often simpler expressions for C∞v molecules are present in the literature. We also note that the polarization in a deformed liquid crystal can be produced as a result of a gradient in the average density of the molecule quadrupole moments [105]. Such a polarization does not need asymmetry in molecular shape. 7.3. Exemplary calculations In order to describe C2 molecules we write σ of the form (n is the symmetry axis) £ ¤ σ = σ0 + σ1 (∆ · n1 )2 + (∆ · n2 )2 + σ2 (n1 · n2 )2 £ ¤ + σ3 (∆ · l1 )2 + (∆ · l2 )2 + σ4 (l1 · l2 )2 £ ¤ + σ5 (l1 · n2 )2 + (l2 · n1 )2 + σ6 (n1 · n2 ) + σ7 (∆ · n1 − ∆ · n2 ) + σ8 [(n1 × n2 ) · ∆] , (84) where the terms with σ6 , σ7 (achiral) and σ8 (chiral) are the lowest order terms compatible with the C2 symmetry. Analytical calculations suggest that the most important for the FE effect is the σ7 term. Let us define ζ = µ3 B3 (T )σ04 (N/V )2 . (85) The FE coefficients can be written as finite sums of the form X [I][J] ai = ζ hF [I] ihF [J] iAi , (86) [I][J] The molecule electric dipole moment is defined as α = x, y, z . a7 = [I][J] 7.2. Flexoelectric coefficients µα (R) = µ1 lα + µ2 mα + µ3 nα , Z (76) The microscopic expressions for the physical FE coefficients are Z a4 = du dR1 dR2 f12 G0 (R1 )µz (R1 )(−ux ) × (U2z − W2x ) , (77) Z a5 = du dR1 dR2 f12 G0 (R1 )µy (R1 )(−ux )U2y , (78) Z a6 = du dR1 dR2 f12 G0 (R1 )µz (R1 )uy W2y , (79) where the coefficients Ai depend on σi only. We used the Lennard–Jones 12-6 potential energy, σ1 /σ0 = 0.5, σ2 /σ0 = 0.1, σ3 /σ0 = −0.4, σ7 /σ0 = 0.1, σ4 = σ5 = σ6 = σ8 = 0, packing N Vmol /V = 0.1. The molecules are similar to cones pointing down. The temperature dependence of the FE coefficients is shown in Fig. 6. The FE coefficients are expressed in µ3 /σ02 , temperatures in ²/kB . On decreasing the temperature we meet the first-order transition from the isotropic phase to the uniaxial nematic phase at T = 0.373 and the second-order transition to the biaxial nematic phase at T = 0.220. The splay FE coefficient e1 splits into a4 > a6 > 0 whereas the bend FE coefficient −e3 splits into 0 > a8 > a9 . The coefficients a5 and a7 are small and negative. If we assume the molecular length σ0 = 1 nm Statistical Theory of Biaxial Nematic and Cholesteric Phases Fig. 6. Temperature dependence of the FE coefficients expressed in µ3 /σ02 . T is the dimensionless temperature. and the electric dipole moment µ = 1 D, then we can estimate the values of the FE coefficients µ/σ02 = 3.3 pC/m. 8. Chiral biaxial nematics An object is chiral if its image in a plane mirror cannot be brought to coincide with itself. A chiral molecule induces chiral interactions that produce intermolecular torques of a given sign and can give rise to equilibrium chiral structures [106]. The effective chiral interactions can arise from classical or quantum interactions [107]. In liquid crystal physics, chiral interactions lead to many different phases, such as cholesteric phases, blue phases or some smectic phases. The (uniaxial) cholesteric phase is characterized by a helical arrangement of the director field N along a common pitch axis £ ¤ N (r) = 0, − sin(q0 x), cos(q0 x) , (87) where q0 = K0 /K2 , P = 2π/|q0 | is the cholesteric pitch and the distortion free-energy density is fd = −K02 /2K2 . The sign of K0 distinguishes between right- and left-handed helices. It is interesting that dimensional analysis does not predict the pitch correctly; the cholesteric phase possesses an additional mesoscopic length scale. Chiral thermotropic mesogens are used in optoelectronic applications due to the unique rheological, electrical, and optical properties. Chiral lyotropic mesogens are also common; for example, DNA, the rod-like fd-virus, or polypeptides. In these systems, the cholesteric pitch is very sensitive to concentration, temperature, pressure, as well as to the solvent conditions. The temperature dependence types for the pitch of the helix can be falling, rising, and inverting. The first of these types is most often found in pure cholesterics, the second is rather rare [108, 109], and the third is observed in mixtures of right-handed and left-handed cholesterics [110]. Theoretical attempts to predict the behavior of the cholesteric pitch are challenging owing to the complexity of the chiral interaction and the inhomogeneous and anisotropic nature of the phase [111, 112]. Early models of the uniaxial cholesteric phase were discussed in [113], where a gener- 361 alized van der Waals theory of lyotropic cholesterics was presented. The influence of fluctuations on the phase diagram of chiral nematic liquid crystals was studied by Englert et al. [114, 115]. The analysis of all five helicity modes was done by Longa et al. [116]. There is a connection between chirality and biaxiality. In 1974 Priest and Lubensky [117] recognized that cholesteric liquid crystals must have some slight biaxial order because of the difference between the directions along and perpendicular to the helical axis. A connection between the phase biaxiality and the stability of various chiral liquid crystalline phases was studied in [118] using the extended de Gennes–Ginzburg–Landau theory. Issaenko and Harris [119] showed that an accurate evaluation of the intermolecular dispersion forces contributing to chiral ordering requires consideration of biaxial correlations between molecules. In 1989 Krolin et al. [120] measured the chirality-induced biaxiality in lyotropic liquid crystals. Recently, Dhakal and Selinger [121] investigated the statistical mechanics of chirality and biaxiality in liquid crystals through a variety of theoretical approaches. Their calculations show that the twist acts as a field on the biaxial order, and conversely, the biaxial order increases the twist (reduces the pitch). In this paper, only one twist is taken into account in the biaxial cholesteric phase. Some experiments with lyotropic liquid crystals suggest [122, 120] that the biaxial cholesteric phase chooses its chiral axis along the axis associated to the smallest eigenvalue of the order parameter tensor Qαβ . It seems that in real experiments boundary conditions can be an important factor determining the phase configuration. We will show how three twists can minimize the distortion free energy. 8.1. Chiral elastic constants The statistical theory of biaxial cholesteric liquid crystals was derived by Kapanowski [123] in the limit of small twists. The cholesteric phase was regarded as a distorted form of the nematic phase. Microscopic expressions for the chiral elastic constants of the biaxial cholesteric phase are as follows: Z 1 βK11 = − dR1 dR2 duf12 ux G0 (R1 )W2y , (88) 2 Z 1 βK22 = − dR1 dR2 duf12 uy G0 (R1 )(U2z − W2x ) , 2 (89) Z 1 βK33 = dR1 dR2 duf12 uz G0 (R1 )U2y . (90) 2 In the case of the uniaxial cholesteric phase Uα = 0 and we get K11 = K22 = K0 , K33 = 0 . (91) 8.2. The equilibrium distortion in the chiral biaxial phase The equilibrium distortion of the biaxial cholesteric phase is composed of three twists with respect to the 362 A. Kapanowski orthogonal axes. In the limit of small distortions the vectors (L, M , N ) have the form (40). Minimizing the distortion free-energy density fd we get the equilibrium twists q1 , q2 , and q3 [123]. It is important that a global minimum of fd exists only if the matrix K1111 K1122 K1133 (92) K1122 K2222 K2233 K1133 K2233 K3333 is positive definite. In the opposite case there is no stable biaxial cholesteric (and nematic) phase. In the limit of the uniaxial cholesteric phase the stability condition has a simple form K1 < 3K2 . (93) As far as we know the condition (93) is satisfied in almost all theories and experiments. The equality K1 = 3K2 appears in some theories with hard molecules [124, 125]. Recently, the elastic constants and orientational viscosities of bent-core nematic liquid crystals have been determined [126] and the unusual anisotropies of these parameters have been discussed in terms of short-range, smectic-C correlations among molecules. The relation (93) was not satisfied and the viscosity ratios also contrasted with ordinary calamitics. Now we would like to explain a certain ambiguous situation concerning the pitch at the transition from the biaxial nematic to the uniaxial nematic phase. If the twists q1 , q2 , and q3 (40) are applied to the biaxial nematic phase and we move to the uniaxial nematic phase, the distortion free energy will be transformed to the form fd = −K0 (q1 + q2 ) + 12 K2 (q1 + q2 )2 − K4 q1 q2 . (94) The q3 twist is not present because it is undefined in the uniaxial nematic phase but there is a surface term with K4 in Eq. (94). The K4 term is present because K1122 = K2 − K4 at the transition. In the thermodynamic limit surface terms are negligible and as the result the equilibrium twist q0 = q1 + q2 = K0 /K2 is obtained. As a side effect, there is a discontinuity in the pitch at the transition. 8.3. Exemplary calculations Let us define ξ = kB T B3 (T )σ04 (N/V )2 . (95) The chiral elastic constants can be written as finite sums of the form X [I][J] Kii = ξ hF [I] ihF [J] iAii , (96) [I][J] [I][J] where the coefficients Aii depend on σi only. We used the Lennard–Jones 12-6 potential energy, σ of the form (22) with σ1 /σ0 = 0.5, σ2 /σ0 = 0.1, σ3 /σ0 = −0.4, σ6 /σ0 = 0.01, σ7 /σ0 = σ8 /σ0 = 0.001, packing N Vmol /V = 0.1. The achiral elastic constants are expressed in ²/σ0 , the chiral elastic constants in ²/σ02 , the twist parameters qi in 1/σ0 , temperatures in ²/kB . Fig. 7. Temperature dependence of the chiral elastic constants expressed in ²/σ02 . T is the dimensionless temperature. Fig. 8. Temperature dependence of the pitches expressed in σ0 . T is the dimensionless temperature. The temperature dependence of the achiral elastic constants is similar to shown in Figs. 3 and 4. The temperature dependence of the chiral elastic constants is shown in Fig. 7. The temperature dependence of the pitches is shown in Fig. 8. In the uniaxial cholesteric phase there is one chiral elastic constant K0 . On decreasing the temperature, K0 splits into K11 > K22 and K33 becomes nonzero. The twist parameters satisfy the inequalities q1 < q2 < q3 and all are greater than q0 . We would like to add that relatively small changes of chiral parameters σ6 , σ7 , and σ8 can change the described picture qualitatively. Generally, the pitches in the biaxial cholesteric phase are smaller than in the uniaxial cholesteric phase. It is in accordance with other approaches to the description of the biaxial cholesteric phase [121]. If we assume the molecular length σ0 = 1 nm and the interaction energy ² = 0.1 eV, then we can estimate the values of the chiral elastic constants ²/σ02 = 0.16 N/m. 9. Other physical properties of nematics In this section we would like to described briefly two physical properties of biaxial nematics where the theory is well developed and interesting. Statistical Theory of Biaxial Nematic and Cholesteric Phases 9.1. Viscous properties Let us consider a nematic film between two parallel glass plates. If the plates slide parallel to each other, the molecules will flow, just as an ordinary fluid. However, due to the coupling between orientation and flow the nematodynamics is more difficult to study than the hydrodynamics of isotropic liquids. The theory relates the viscous stress tensor σαβ to viscosity coefficients, the velocity gradient tensor and the directors. A complete set of viscosity coefficients for the uniaxial nematic consists of five elements [127]: three Mięsowicz coefficients η1 , η2 , and η3 , the coefficient η12 , and the rotational viscosity coefficient γ1 . Theoretical treatments of the viscous properties often use different sets of coefficients and the most common is a set of five αi introduced by Leslie [128, 129]. The most complete formulation of the viscosity theory for biaxial liquid crystals was given by Leslie et al. in Ref. [130], where twelve independent viscosity coefficient were identified (three rotational viscosities and nine effective shearing viscosities). The method of calculation of the Leslie viscosity coefficients was given by Fiałkowski for the case of ellipsoids with three different principal axes [131, 132]. The derived formulae were expressed in terms of order parameters, temperature, number density, and diffusion constants. 9.2. Defects A nematic liquid crystal can be regarded as an ordered medium, where a region of space is described by a function ψ that assigns to every point of the region of an order parameter [133]. The possible values of the order parameter constitute a space known as the order parameter space M , manifold of internal states or a space of degeneracy [134]. The region of space is usually three-dimensional space, though two-dimensional regions can be used for films. Distortions of ψ can contain singularities, where ψ is not defined. For a three-dimensional medium, the singular regions might be either zero-dimensional (points), one-dimensional (lines), or two-dimensional (walls). These are defects. Whenever an inhomogeneous state cannot be eliminated by continuous variations of the order parameter, it is called a topological defect. If the inhomogeneous state does not contain singularities, but nevertheless is not deformable continuously into a homogeneous state, one says that the system contains a topological configuration or soliton [135]. Defects are local breakings of symmetry in an ordered medium. Dislocations break translational symmetries. Disclinations break rotational symmetries and are basic defects in liquid crystals [136]. An adequate description of defects in ordered media requires the theory of homotopy, which is part of algebraic topology [137]. Then it is possible to classify the defects types, to find the laws of decay, merger and crossing of defect, to trace out their behaviour during phase transitions, and so on. The key point is the concept of topological invariant or topological charge which is inherent 363 in every defect. The stability of the defect is guaranteed by the conservation of its topological charge. Topologically stable wall, line, and point defects are described by elements of the homotopy groups π0 (M ), π1 (M ), and π2 (M ), respectively. The situation becomes more complex when we wish to study surface defects, namely defects that arise on a surface bonding the body [138]. The order parameter space for uniaxial nematics is the surface of unit sphere, where opposite points are identified. This space is known as the real projective plane M = RP 2 . M can be expressed as the coset space M = SO(3)/D∞ . The important results are π0 (M ) = Z2 , π1 (M ) = Z2 , π2 (M ) = Z. The order parameters space for biaxial nematics is the full proper three-dimensional rotation group SO(3), provided the appropriate discrete sets of points are identified (proper point group D2 ). Thus M = SO(3)/D2 but it turns out that the natural representation is M = SU(2)/H, where H is the quaternion group. Here π1 (M ) = H and there are five types of line defects because there are five classes of conjugated elements in H. The groups π0 (M ) and π2 (M ) are trivial. The structures, energies, and interactions of defects in many liquid crystalline phases were presented by Chandrasekhar and Ranganath [139]. In 1997 Biscari and Peroli studied the topological properties of manifolds which describe uniaxial and biaxial nematics in order to compare and classify all types of bulk and surface defects in nematics [138]. The manifold of the traceless order tensor Q consists of three disjoint subsets with different symmetries (isotropic, uniaxial, and biaxial distributions). The surfaces induce in the order tensor the additional symmetry which is important for surface defects. Information about the nature of surface defects come from the structure of the relative homotopy groups. 10. Conclusions In this paper the statistical theory of biaxial nematic and cholesteric phases was presented. It was derived in the thermodynamic limit at small density and small distortions. The microscopic expressions for many macroscopic parameters were shown: for the elastic constants, the dielectric susceptibilities, the FE coefficients, and the phase twists. The expressions involve the one-particle distribution function and the potential energy of two-body short-range interactions. In addition flow properties and defects of biaxial nematics were briefly discussed. Exemplary calculations were shown for several systems of molecules interacting via Corner-type potential based on the Lennard–Jones 12-6 functional dependence. Estimates for all parameters were given as the direct comparison between the theory and the experiment is in many cases difficult at present. The presented theory can be generalized in many directions. It is desirable to go beyond the low density limit where the Mayer function f12 is replaced with a better approximation of the direct correlation function. It is 364 A. Kapanowski known that to understand the thermodynamic properties of dense fluids, the repulsive forces must be treated to all orders whereas attractive forces can be treated by a perturbation theory [140]. Let us consider the convex peg model [140, 141], where molecules are envisioned to have a hard biaxial core embedded in a spherically symmetric square well. Then the interaction term in the free energy (27) can be written as Z N2 βFint = dR1 dR2 f (R1 )f (R2 ) 2V ¸ ½ · 4 − 3ρ β² + e −1 × Vexcl (R1 , R2 ) 4(1 − ρ)2 ¾ + (1 − eβ² )Vsw , (97) where ρ = N Vmol /V , Vexcl (R1 , R2 ) is the excluded vol3 is the exume of the hard core, and Vsw = (4/3)πσsw cluded volume of the spherical square well. Hard core interactions are mapped into the dense hard sphere system using a Carnahan–Starling approximation. If we consider the Corner square well potential energy (104), then by means of similar assumptions we derive Z N2 βFint = dR1 dR2 f (R1 )f (R2 )K(R1 , R2 )(σ03 /3) 2V · ¸ 4 − 3ρ × − 3B2 (T ) − 1 . (98) 4(1 − ρ)2 The theory can be generalized to the case of many component liquid crystalline mixtures. In fact, mixtures are essential for technological applications of liquid crystals because various combinations of physical properties are required [51]. The design and synthesis of materials exhibiting the biaxial nematic phase with desired properties constitute one of the major challenges in liquid crystal science. Appendix A Invariants (j) We list the main properties of the functions Fµν , which describe biaxial and uniaxial molecules in biaxial and uniaxial phases [94]. (j) (i) The definition of the functions Fµν : µ ¶2+δµ0 +δ0ν 1 (j) Fµν (R) = √ 2 X (j) × (−1)j(σ−ρ)/2 Dρ×µ,σ×ν (R) , (99) ρ,σ=±1 where j is a non-negative integer. If j is even, then 0 ≤ µ ≤ j and 0 ≤ ν ≤ j. If j is odd, then 2 ≤ µ ≤ j (j) and 2 ≤ ν ≤ j. Functions Dµν (R) are standard rotation matrix elements. (j) (ii) The functions Fµν are real. (j) (j) (iii) For the inverse rotation Fµν (R−1 ) = Fνµ (R). For (j) the zero rotation Fµν (0) = δµν . (iv) The orthogonality relations have the form Z (j) (k) dRFµν (R)Fρσ (R) = δjk δµρ δνσ 8π 2 /(2j + 1) . (100) (v) The invariants can be labelled by the indicator [I], where I runs from 1 by 1 to infinity. For two indicators [I] = (j, µ, ν) and [J] = (k, ρ, σ) we can write that [I] < [J] only when j < k, or when j = k and µ < ρ, or when j = k, µ = ρ, and ν < σ. Appendix B Corner potential energy We list some examples for the dependence of the Corner potential energy Φ12 on u/σ, together with the functions Bs (T ) defined as [β = 1/(kB T )]: Z ∞ Z ∞ Bs (T ) = dxxs f12 (x) = dxxs [exp(−βΦ12 (x)) − 1] . 0 0 (101) (i) The soft-core potential energy (m ≥ s + 2): Φ12 (u/σ) = ²(σ/u)m , µ ¶ m−s−1 −1 Γ (β²)(s+1)/m . Bs (T ) = s+1 m (ii) The square-well potential energy +∞ for (u/σ) < 1, Φ12 (u/σ) = −² for 1 < (u/σ) < Rsw , 0 for (u/σ) > Rsw , ª 1 © s+1 Bs (T ) = [exp(β²) − 1](Rsw − 1) − 1 . s+1 (iii) The soft-well potential energy (m ≥ s + 2) ( +∞ for (u/σ) < 1, Φ12 (u/σ) = −²(σ/u)m for (u/σ) > 1, µ ¶ s+1 s+1 −1 ,1 − , β² . Bs (T ) = 1 F1 − s+1 m m (102) (103) (104) (105) (106) (107) (iv) The Lennard–Jones m–n potential energy (m > n > s + 1): £ ¤ Φ12 (u/σ) = 4² (σ/u)m − (σ/u)n , (108) ∞ Bs (T ) = µ −1 X 1 (4β²)[p(m−n)+s+1]/m s + 1 p=0 p! pn + m − s − 1 ×Γ m + ¶ ∞ 1 X 1 n (4β²)[p(m−n)+m−n+s+1]/m s + 1 p=0 p! m µ ×Γ ¶ pn + n − s − 1 . m (109) Statistical Theory of Biaxial Nematic and Cholesteric Phases Fs (0, 0, 0) = 1/3 Appendix C Shape factors for biaxial molecules We list shape factors for model biaxial molecules with the D2h symmetry. The reference point, where we calculate the electric field, is always placed in the molecule mass center. (i) Ellipsoids with the axes 2a × 2b × 2c: Ω1 = J(1 − a2 /b2 , 1 − a2 /c2 ) , (110) 2 2 2 2 (111) 2 2 2 2 (112) Ω2 = J(1 − b /a , 1 − b /c ) , Ω3 = J(1 − c /a , 1 − c /b ) , where the function J(x, y) is defined for x < 1 and y < 1 as Z 1 dt t2 p J(x, y) = . (113) (1 − xt2 )(1 − yt2 ) 0 (ii) Platelets of dimensions a × b × c: Ω1 = Fp (a, b, c) , Ω2 = Fp (b, a, c) , Ω3 = Fp (c, a, b) , Fp (x, y, z) = Fp (x, z, y) ! Ã 2 yz = arcsin p . π (x2 + y 2 )(x2 + z 2 ) (114) (115) Let us note that Ω1 + Ω2 + Ω3 = 1. (iii) Spherocuboids [65] obtained by the Minkowski addition of a parallelepiped with sides of length a, b, and c, and a sphere of radius R. As a special case one can obtain platelets, spheroplatelets or spherocylinders. A spherocuboid can be bounded by the box of dimensions (a + 2R) × (b + 2R) × (c + 2R). Ω1 = Fp (a + 2R, b, c) + Fs (a, b, c) + Fc (a, b, c) + Fc (a, c, b) , (116) Ω2 = Fp (b + 2R, a, c) + Fs (b, a, c) + Fc (b, a, c) + Fc (b, c, a) , (117) Ω3 = Fp (c + 2R, a, b) + Fs (c, a, b) + Fc (c, a, b) + Fc (c, b, a) , Z π/2 Fs (x, y, z) = dφ (118) 0 Z π/2 × 0 8R2 (x + 2R cos θ) sin θ cos θ dθ , πξs3 (119) ξs2 = (x + 2R cos θ)2 + (y + 2R sin θ sin φ)2 + (z + 2R sin θ cos φ)2 , (120) Z π/2 4R cos θ(x + 2R cos θ)z p , (121) Fc (x, y, z) = dθ πξc2 z 2 + ξc2 0 ξc2 = (x + 2R cos θ)2 + (y + 2R sin θ)2 , (122) like for spheres , z Fc (0, 0, z) = √ like for cylinders . 2 z 2 + 4R2 365 (123) (124) References [1] P.G. de Gennes, J. Prost, The Physics of Liquid Crystals, Clarendon Press, Oxford 1993. [2] M.J. Freiser, Phys. Rev. Lett. 24, 1041 (1970). [3] M.J. Freiser, Mol. Cryst. Liq. Cryst. 14, 165 (1971). [4] R. Alben, J. Chem. Phys. 59, 4299 (1973). [5] J.P. Straley, Phys. Rev. A 10, 1881 (1974). [6] G.R. Luckhurst, C. Zannoni, P.L. Nordio, U. Segre, Mol. Phys. 30, 1345 (1975). [7] D.K. Remler, A.D.J. Haymet, J. Phys. Chem. 90, 5426 (1986). [8] B. Bergersen, P. Palffy-Muhoray, D.A. Dunmur, Liq. Cryst. 3, 347 (1988). [9] Z. Zhi-Dong, H. Xi-Min, Acta Phys. Sin. (Overseas Ed.) 6, 671 (1997). [10] D.W. Allender, M.A. Lee, Mol. Cryst. Liq. Cryst. 110, 331 (1984). [11] D.W. Allender, M.A. Lee, N. Hafiz, Mol. Cryst. Liq. Cryst. 124, 45 (1985). [12] E.F. Gramsbergen, L. Longa, W.H. de Jeu, Phys. Rep. 135, 195 (1986). [13] D. Allender, L. Longa, Phys. Rev. E 78, 011704 (2008). [14] G.R. Luckhurst, S. Romano, Mol. Phys. 40, 129 (1980). [15] C.D. Mukherjee, N. Chatterjee, Phys. Lett. A 189, 86 (1994). [16] F. Biscarini, C. Chiccoli, P. Pasini, F. Semeria, C. Zannoni, Phys. Rev. Lett. 75, 1803 (1995). [17] C. Chiccoli, P. Pasini, F. Semeria, C. Zannoni, Int. J. Mod. Phys. C 10, 469 (1999). [18] B.M. Mulder, Th.W. Ruijgrok, Physica A 113, 145 (1982). [19] J. Alfaro, O. Cubero, L.F. Urrutia, Physica A 272, 314 (1999). [20] C.-S. Shih, R. Alben, J. Chem. Phys. 57, 3055 (1972). [21] B.M. Mulder, Liq. Cryst. 1, 539 (1986). [22] B. Mulder, Phys. Rev. A 39, 360 (1989). [23] R. Holyst, A. Poniewierski, Mol. Phys. 69, 193 (1990). [24] M.P. Taylor, J. Herzfeld, Phys. Rev. A 44, 3742 (1991). [25] A.N. Zakhlevnykh, P.A. Sosnin, Mol. Cryst. Liq. Cryst. 293, 135 (1997). [26] M.P. Allen, Liq. Cryst. 8, 499 (1990). [27] P.J. Camp, M.P. Allen, J. Chem. Phys. 106, 6681 (1997). [28] R. Berardi, L. Muccioli, S. Orlandi, M. Ricci, C. Zannoni, J. Phys., Condens. Matter 20, 463101 (2008). [29] L.J. Yu, A. Saupe, Phys. Rev. Lett. 45, 1000 (1980). [30] H.F. Leube, H. Finkelmann, Makromol. Chem. 192, 1317 (1991). 366 A. Kapanowski [31] T. De’Neve, M. Kleman, P. Navard, J. Phys. II 2, 187 (1992). [32] R. Berardi, C. Zannoni, J. Chem. Phys. 113, 5971 (2000). [33] G.R. Luckhurst, Thin Solid Films 393, 40 (2001). [34] G.R. Luckhurst, Nature (London) 430, 413 (2004). [35] L.A. Madsen, T.J. Dingemans, E.T. Samulski, Phys. Rev. Lett. (2004). M. Nakata, 92, 145505 [36] B.R. Acharya, A. Primak, S. Kumar, Phys. Rev. Lett. 92, 145506 (2004). [37] L.A. Madsen, T.J. Dingemans, E.T. Samulski, Phys. Rev. Lett. (2006). M. Nakata, 96, 219804 [38] K. Merkel, A. Kocot, J.K. Vij, R. Korlacki, G.H. Mehl, T. Meyer, Phys. Rev. Lett. 93, 237801 (2004). [39] J.L. Figueirinhas, C. Cruz, D. Filip, G. Feio, A.C. Ribeiro, Y. Frere, T. Meyer, G.H. Mehl, Phys. Rev. Lett. 94, 107802 (2005). [40] K. Neupane, S.W. Kang, S. Sharma, D. Carney, T. Meyer, G.H. Mehl, D.W. Allender, S. Kumar, S. Sprunt, Phys. Rev. Lett. 97, 207802 (2006). [41] T.C. Lubensky, L. Radzihovsky, Phys. Rev. E 66, 031704 (2002). [42] L.G. Fel, Phys. Rev. E 52, 702 (1995). [43] L.G. Fel, Phys. Rev. E 52, 2692 (1995). [44] A. Kapanowski, to be published in Mol. Cryst. Liq. Cryst. [45] B. Mettout, Phys. Rev. E 74, 041701 (2006). [46] H. Takezoe, Y. Takanishi, Jpn. J. Appl. Phys. 45, 597 (2006). [47] R. Alben, J. Chem. Phys. 59, 4299 (1973). [48] R. Hashim, G.R. Luckhurst, S. Romano, Mol. Phys. 53, 1535 (1984). [49] P.J. Camp, M.P. Allen, P.G. Bolhuis, D. Frenkel, J. Chem. Phys. 106, 9270 (1997). [50] F.M. van der Kooij, H.N.W. Lekkerkerker, Phys. Rev. Lett. 84, 781 (2000). [51] C. Tschierske, D.J. Photinos, J. Mater. Chem. 20, 4263 (2010). [52] D. Apreutesei, G. Mehl, Chem. Commun., 609 (2006). [53] A. Cuetos, A. Galindo, G. Jackson, Phys. Rev. Lett. 101, 237802 (2008). [54] J.-H. Lee, T.-K. Lim, W.-T. Kim, J.-I. Jin, J. Appl. Phys. 101, 034105 (2007). [55] R. Berardi, L. Muccioli, C. Zannoni, J. Chem. Phys. 128, 024905 (2008). [56] S.D. Peroukidis, P.K. Karahaliou, A.G. Vanakaras, D.J. Photinos, Lig. Cryst. 36, 727 (2009). [57] M. Fiałkowski, A. Kapanowski, K. Sokalski, Mol. Cryst. Liq. Cryst. 265, 371 (1995). [58] B.S. John, C. Juhlin, F.A. Escobedo, J. Chem. Phys. 128, 044909 (2008). [59] L. Onsager, Ann. New York Acad. Sci. (1949). [60] D. Frenkel, Mol. Phys. 60, 1 (1987). 51, 627 [61] W. Maier, A. Saupe, Z. Naturforsch. A 13, 564 (1958); 14, 882 (1959); 15, 287 (1960). [62] G.S. Singh, B. Kumar, Ann. Phys. (New York) 294, 24 (2001). [63] B. Tjipto-Margo, G.T. Evans, J. Chem. Phys. 94, 4546 (1991). [64] M.P. Taylor, Liq. Cryst. 9, 141 (1991). [65] B.M. Mulder, Mol. Phys. 103, 1411 (2005). [66] R. Rosso, E.G. Virga, Phys. Rev. E 74, 021712 (2006). [67] G.S. Singh, B. Kumar, J. Chem. Phys. 105, 2429 (1996). [68] C. Vega, Mol. Phys. 92, 651 (1997). [69] C. McBride, E. Lomba, Fluid Phase Equilibria 255, 37 (2007). [70] X.-M. You, A.Yu. Vlasov, A.J. Masters, J. Chem. Phys. 123, 034510 (2005). [71] A.M. Sonnet, E.G. Virga, G.E. Durand, Phys. Rev. E 67, 061701 (2003). [72] L. Longa, P. Grzybowski, S. Romano, E. Virga, Phys. Rev. E 71, 051714 (2005). [73] G. De Matteis, E.G. Virga, Phys. Rev. E 71, 061703 (2005). [74] B.J. Berne, P. Pechukas, J. Chem. Phys. 56, 4213 (1972). [75] J.G. Gay, B.J. Berne, J. Chem. Phys. 74, 3316 (1981). [76] G. Ayton, G.N. Patey, J. Chem. Phys. 102, 9040 (1995). [77] R. Berardi, C. Fava, C. Zannoni, Chem. Phys. Lett. 236, 462 (1995). [78] R. Berardi, C. Fava, C. Zannoni, Chem. Phys. Lett. 297, 8 (1998). [79] D.J. Cleaver, C.M. Care, M.P. Allen, M.P. Neal, Phys. Rev. E 54, 559 (1996). [80] J. Corner, Proc. R. Soc. Lond. A 192, 275 (1948). [81] L. Blum, A.J. Torruella, J. Chem. Phys. 56, 303 (1972). [82] J.D. Parsons, Phys. Rev. A 19, 1225 (1979). [83] Sin-Doo Lee, J. Chem. Phys. 87, 4972 (1987). [84] A. Cuetos, B. Martínez-Haya, S. Lago, L.F. Rull, J. Phys. Chem. B 109, 13729 (2005). [85] J. Stecki, A. Kloczkowski, J. Phys. Paris 40, C3-360 (1979). [86] L. Longa, J. Stelzer, D. Dunmur, J. Chem. Phys. 109, 1555 (1998). [87] E. Govers, G. Vertogen, Phys. Rev. A 30, 1998 (1984). [88] D. Monselesan, H.-R. Trebin, Phys. Status Solidi B 155, 349 (1989). [89] L. Longa, H.-R. Trebin, Phys. Rev. A 39, 2160 (1989). [90] S. Stallinga, G. Vertogen, Phys. Rev. E 49, 1483 (1994). [91] A. Saupe, J. Chem. Phys. 75, 5118 (1981). [92] Y. Singh, K. Rajesh, V.J. Menon, S. Singh, Phys. Rev. E 49, 501 (1994). Statistical Theory of Biaxial Nematic and Cholesteric Phases [93] U.D. Kini, S. Chandrasekhar, Physica A 156, 364 (1989). [94] A. Kapanowski, Phys. Rev. E 55, 7090 (1997). [95] W. Maier, G.Z. Meier, Z. Naturforsch. A 16, 1200 (1961). [96] A. Kapanowski, T. Wietecha, Phys. Rev. E 71, 021710 (2005). [97] A. Kapanowski, Mol. Cryst. Liq. Cryst. 436, 1091 (2005). [98] P. Bordewijk, Physica 75, 146 (1974). [99] K. Toriyama, D.A. Dunmur, S.E. Hunt, Liq. Cryst. 5, 1001 (1989). [100] R.B. Meyer, Phys. Rev. Lett. 22, 918 (1969). [101] A.G. Petrov, in: Physical Properties of Liquid Crystals, EMIS Datareview Series 25, Eds. D.A. Dunmur, A. Fukuda, G.L. Luckhurst, IEE, London 2001. [102] A. Kapanowski, Phys. Rev. E 75, 031709 (2007). [103] A. Kapanowski, Opto-Electron. Rev. 16, 9 (2008). [104] A. Kapanowski, Phys. Rev. E 77, 052702 (2008). [105] J. Prost, J.P. Marcerou, J. Phys. (Paris) 38, 315 (1977). [106] A.B. Harris, R.D. Kamien, T.C. Lubensky, Phys. Rev. Lett. 78, 1476 (1997). [107] S.A. Issaenko, A.B. Harris, T.C. Lubensky, Phys. Rev. E 60, 578 (1999). [108] T. Harada, P. Crooker, Mol. Cryst. Liq. Cryst. 30, 79 (1975). [109] P.M. Agics, L. Bata, Mol. Cryst. Liq. Cryst. 107, 391 (1984). [110] A.N. Zakhlevnykh, M.I. Shliomis, Sov. Phys. JETP 59, 764 (1984). [111] A.B. Harris, R.D. Kamien, T.C. Lubensky, Rev. Mod. Phys. 71, 1745 (1999). [112] G.S. Chilaya, L.N. Lisetskii, Sov. Phys. Usp. 24, 6 (1981). [113] H.H. Wensink, G. Jackson, J. Chem. Phys. 130, 234911 (2009). [114] J. Englert, L. Longa, H. Stark, H.-R. Trebin, Phys. Rev. Lett. 81, 1457 (1998). [115] J. Englert, H. Stark, L. Longa, H.-R. Trebin, Phys. Rev. E 61, 2759 (2000). [116] L. Longa, M. Cieśla, H.-R. Trebin, Phys. Rev. E 67, 061705 (2003). [117] R.G. Priest, T.C. Lubensky, Phys. Rev. A 9, 893 (1974). 367 [118] L. Longa, W. Fink, H.-R. Trebin, Phys. Rev. E 50, 3841 (1994). [119] S.A. Issaenko, A.B. Harris, Phys. Rev. E 61, 2777 (2000). [120] T. Krolin, A.M. Figueiredo Neto, L. Liebert, Y. Galerne, Phys. Rev. A 40, 4647 (1989). [121] S. Dhakal, J.V. Selinger, Phys. Rev. E 83, 020702(R) (2011). [122] A.M. Figueiredo Neto, Y. Galerne, L. Liebert, J. Phys. Chem. 89, 3939 (1985). [123] A. Kapanowski, Z. Naturforsch. 57a, 105 (2002). [124] J.P. Straley, Phys. Rev. A 8, 2181 (1973). [125] H. Kimura, M. Hosino, H. Nakano, J. Phys. Soc. Jpn. 51, 1584 (1982). [126] M. Majumdar, P. Salamon, A. Jákli, J.T. Gleeson, S. Sprunt, Phys. Rev. E 83, 031701 (2011). [127] Liquid Crystals: Experimental Study of Physical Properties and Phase Transitions, Ed. S. Kumar, Cambridge University Press, Cambridge (UK) 2001. [128] F.M. Leslie, Q.J. Mech. Appl. Math. 19, 367 (1966). [129] F.M. Leslie, Arch. Ration. Mech. Anal. 28, 265 (1968). [130] F.M. Leslie, J.S. Laverty, T. Carlsson, Quart. Jour. Mech. Appl. Math. 45, 595 (1992). [131] M. Fiałkowski, Phys. Rev. E 55, 2902 (1997). [132] M. Fiałkowski, Phys. Rev. E 58, 1955 (1998). [133] N.D. Mermin, Rev. Mod. Phys. 51, 591 (1979). [134] H.-R. Trebin, Adv. Phys. 31, 194 (1982). [135] M. Kleman, O.D. Lavrentovich, Soft Matter Physics: An Introduction, Series: Partially Ordered Systems, Springer-Verlag, New York 2003. [136] M. Kleman, Rep. Prog. Phys. 52, 555 (1989). [137] L. Michel, Rev. Mod. Phys. 52, 617 (1980). [138] P. Biscari, G.G. Peroli, Commun. Math. Phys. 186, 381 (1997). [139] S. Chandrasekhar, G.S. Ranganath, Adv. Phys. 35, 507 (1986). [140] B. Tjipto-Margo, G.T. Evans, Mol. Phys. 74, 85 (1991). [141] D.C. Williamson, F. del Rio, J. Chem. Phys. 109, 4675 (1998).