ANALYSIS OF A MODEL FOR BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES. Tiziana Giorgi and Feras Yousef. (Communicated by Patricia Bauman)

Transcription

1 DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS SERIES B Volume X, Number X, XX 15 doi:1.3934/dcdsb.15.xx.xx pp. X XX ANALYSIS OF A MODEL FOR BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES Tiziana Giorgi and Feras Yousef Department of Mathematical Sciences New Mexico State University Las Cruces, NM , USA Communicated by Patricia Bauman Abstract. We consider a model originally introduced to study layer-undulated structures in bent-core molecule liquid crystals. We first prove existence of minimizers, then analyze a simplified version used to study how in columnar phases the width of the column affects the type of switching, which occurs under an applied electric field. We show via Γ-convergence that as the width of the column tends to infinity, rotation around the tilt cone is favored, provided the coefficient of the coupling term, between the polar parameter, the nematic parameter, and the layer normal is large. 1. Introduction. Bent-core molecule liquid crystals BLCs are achiral molecule materials, which exhibit phases typical of chiral liquid crystals LCs. Due to this characteristic, they have been in recent years the subject of a noticeable amount of theoretical and experimental work in the physics community, see 7, 13] for a review of the available literature on BLCs. Among other approaches, phenomenological Landau models have been successfully employed to analyze and numerically approximate the smectic behavior of these amazingly phase-rich materials. The main focus of this paper is a version, presented in 8, 5], of the free energy first introduced by Vaupotič and Čopič in 16, 5]. Vaupotič and Čopič s interest in 16] lies on the polarization-modulated layer-undulated structure found in some tilted smectic phases, and the key feature of their model is the coupling term between the polarization splay and the tilt of the material molecules with respect to the normal vector of the smectic layers. In 8], instead the authors look at the type of switching under applied electric field encountered in the so-called columnar phases, where the material forms two-dimensional structures made of fragments of smectic layers, which is seen experimentally in BLCs. We first introduce a slight modification of the full energy functional proposed in 16], needed to guarantee mathematical well-posedness, and prove existence of minimizers. We then set up the mathematical framework to study the reduced functional considered in 8]. These are inherently one-dimensional problems, which nevertheless are of mathematical interest as the two unknown functions describing 1 Mathematics Subject Classification. Primary: 76A15; Secondary: 49J99. Key words and phrases. Bent-core molecules, liquid crystals, columnar phases, Γ-convergence, energy minimization. Funding to the first author was provided by the NSF grant #DMS The research of the second author was partially funded by the NSF grant #DMS

2 TIZIANA GIORGI AND FERAS YOUSEF the state of the BLC are mutually orthogonal unit vectors. We use Γ-convergence to derive a limiting energy functional for the large-column width case, and argue that there are only two candidates for minimizers of the limiting functional, which respectively represent switching by rotation on the tilt cone, and along the axis of the molecule, see Remark. We obtain some basic estimates to show that the state representing switching by rotation on the tilt cone is in this limiting case the minimizer for certain values of the physical parameters. The paper is organized as follows. In Section, we introduced the relevant terminology, and provide some background to illustrate the physical picture underlying the problem. We present in Section 3 the full version of the model, suitably modified, and we show existence of minimizers. Section 4 is devoted to the study of columnar phases in the case of large columnar width. In Section 5, we study the minimizers of the Γ-limit energy and obtain some simple energy estimates, which are functions of the physical parameters.. Background. In a very simplified illustration, we think to temperature as providing a measure of random motion, with its increase affecting the ability of attractive intermolecular forces to keep order. Solids can be then pictured as regularly stacked molecules having positional and orientational order, and liquids as molecules free to randomly move, with no positional or orientational order. For thermotropic liquid crystals, the transition from crystalline states to isotropic liquid states, which is induced by varying temperature, is characterized by in between thermodynamically anisotropic fluid-like stable phases. In these in between phases, liquid crystals may retain some of the mechanical properties of fluids, but more importantly they retain some of the optical features of crystals, hence their use and importance in the manufacturing of display devices. Figure 1 provides a standard simple microscopic Figure 1. Schematic thermotropic liquid crystals phases. In the inset, the smectic layer normal ν is taken to be parallel to the vertical axis z. picture of this behavior, which originates from the fact that many LC materials have rod-like strongly elongated in one direction molecules. The key point to keep in mind is that in a liquid crystalline phase some positional or all order may be lost,

3 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 3 but some orientational order remains, and one can define a preferred local average direction, represented by a unit vector, n, called the director, see inset of Figure 1. The liquid crystalline phases are grossly divided in two families: nematic and smectic. The nematic phase, N, is distinguished by the presence of some degree of long-range orientational order, but no positional order. The molecules can be polar but the states described by n and n are physically the same. The director in the variation chiral nematic phase, N*, rotates in a helical sense throughout the sample, see Figure. The smectic phases manifest also a small amount of orientational order, but now also some degree of positional order is present: the molecules are arranged on average in equidistant planes, and one has a layered structure with a well-defined interlayer spacing. The smectic layer normal is denoted by ν. Two common smectic phases are the smectic A SmA, and the smectic C SmC, with their chiral variants, SmA* and SmC*. In the SmA phase, the axes of orientational and positional order are parallel, while in the SmC phase orientational and positional orders are tilted one with respect to the other, see inset of Figure 1, and so the director n is tilted at a defined angle, the tilted angle θ, to the layer normal ν. In general, the tilt magnitude θ is well-defined at given values of pressure and temperature, while the azimuthal or phase angle φ is not, and n is located somewhere on the cone of axis ν and opening angle θ, called the tilt cone. The chiral SmC* phase is characterized by an helical superstructure, the helix axis is along ν while the azimuthal angle φ is the quantity changing when proceeding along the axis, see Figure. Figure. Schematic of chiral phases. The SmC* is a particularly interesting phase in soft matter, since chirality has significant consequences, which include the appearance of a mesoscopic electric polarization. In this phase the polarization is parallel to n ν, and it too changes direction as n travels along the helical axis, thus typically, the macroscopic polarization is zero. Nevertheless, under proper manipulation some LCs can exhibits ferroelectric behavior, which is of great interest for technological applications. Usually, a LC which displays this phase has chiral molecules, that is molecules that differ from their mirror image, but still have no head or tail. Recently it was discovered that macroscopically chiral smectic phases are observed in certain achiral compounds, where the bent geometry of the molecules gives rise to a polar order

4 4 TIZIANA GIORGI AND FERAS YOUSEF perpendicular to the the molecular axis, see Figure 3. Cost efficiency manufacturing considerations, and richness in the variety of phases observed, have made these materials the subject of many theoretical and experimental physical studies, 7]. The mathematical treatments are nevertheless still very scant, although validation and analysis of the available models might result in new relevant contributions to the area, see for example 3]. Figure 3. The spontaneous polarization is along the direction given by the polarization director p. Bent-core molecule liquid crystals, schematically depicted in Figure 3, have a characteristic banana shape, and can be described by the tilt angle θ, the azimuthal angle φ, the layer tilt angle, and the polar angle α, 16, 3]. We will use the following convention for these fundamental angles. We take the layer tilt angle to be defined modulo π as the angle from the z-axis toward the x-axis about the positive y-axis, with this choice the layer normal is given by ν =< sin,, cos >. The tilt angle θ is the angle between n and ν, that is n ν = cos θ, and has values θ π. We define the azimuthal angle φ as the angle from t toward s about ν, where t is the unit vector in the x-z plane perpendicular to ν, and such that t ν is in the negative y direction; while s is a unit vector in the ν-n plane, perpendicular to ν and such that t ν is in the same direction as n ν. Finally, the polar angle α is the angle obtained from n ν to p about n, with α π, π. Note that α and φ are not defined if n and ν are parallel, in this work, when using this representation, we will always assume that θ, π. We can write the layer normal, and the nematic and polar directors in terms of the molecular angles as follows: ν =< sin,, cos >; n =< sin cos θ sin θ cos φ cos, sin θ sin φ, cos θ cos sin θ sin cos φ >; p =< sin α sin θ sin sin α cos θ cos cos φ sin φ cos cos α, 1 sin α cos θ sin φ cos φ cos α, sin α sin θ cos sin α cos θ cos φ sin sin φ sin cos α >.

5 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 5 This representation contains implicitly the constraints n = p = 1 and n p =, a fact which proves to be particularly useful when deriving bounds for the energy of specific states. A popular reference book for the physics of LCs is 13] by de Gennes and Prost, while a mathematical introduction to classical results for nematic and smectic LCs are the books by Stewart 15], and Virga 17]. For an understanding of the properties and modeling of ferroelectric LCs standard references include Lagerwall 1] and Mu sevi c & Blinc & Zek s 1]. 3. Model. The Landau free energy density introduced in 16] includes some usual terms, that is the one-constant approximation Frank-Oseen nematic elastic energy, the Chen-Lubensky energy describing the smectic structure and the dipole-dipole self-interaction energy, as well as three terms pertinent to the bent-core molecules: a coupling term between the polarization splay and the tilt of the material molecules with respect to the normal vector of the smectic layers, a term proportional to the polarization splay and finally a coupling term between the polar and nematic parameters and the layer normal, which serves as penalization for deviation of the polarization from being parallel to n ν. Specifically, they propose the free energy density: 1 K n n n ] c n iq ψ c n ψ D n ψ P ɛɛ p 1 K p p n ψ 1 K p p K np p n ψ, with the constraints n = p = 1 and n p =. Here, q = π d where d is the smectic layer thickness in the non-tilted phase, p 1 is the x-component of the polarization director, and all the coefficients are positive, except c which is taken to be strictly less than zero, as it is the case in the SmC phase. In postulating the form of the energy, the authors make the assumption that the problem of layer undulation formation starting from an homogenous tilted phase is essentially one-dimensional, therefore in the above expression everything depends only on the variable x. The authors also assume that the polarization P = P p has constant magnitude P, and neglect the Maxwell s equations in the electric self-interaction energy term, hence the simplified form of the self-interaction dipole term. An additional implicit assumption is that the material is well in the smectic phase so that ψ is constant. The energy for bent-core molecules 3, 1] and the one for ferroelectricity 4, 13] used by Bauman and collaborators incorporate the full self-interaction energy term, but assume that the ground state of the system is the so-called SmC G state, resulting in a different choice for the term describing the interaction between the polarization director, the layer normal and nematic director. There are two mathematical shortcomings in the above energy. First, to show existence of energy minimizers it is well-known that the term involving the second order derivatives of the Chen-Lubensky smectic energy 6] needs to be modified, for example following the work of Luk yanchuk 11], since in its original form does not allow for the control of the second derivatives of ψ in the n direction, 9]. Second, the quadratic term involving the divergence of p is not sufficient to mathematically control its gradient.

6 6 TIZIANA GIORGI AND FERAS YOUSEF We use the Chen-Lubensky energy density introduced in 3], and consider the modifications mentioned above, so that the energy density in becomes fψ, n, p = 1 K n n n ] c D ψ c D ψ a D D ψ a D D ψ P p 1 ɛɛ K p p D ψ 1 K p p p ] K np p n D ψ, 3 where Dψ = ψ iq cos θ B n ψ, and D, D are its components parallel and perpendicular to the direction of the nematic vector, that is D ψ = n ψ iq cos θ B ψn and D ψ = Dψ D ψ. Here, q is the smectic layer periodicity along the z-axis, and q = q cos θ B, with θ B the cone angle in the homogeneous phase. The order parameter ψ describes the smectic layering, and corresponds to the distribution of the centers of mass of the molecules, more precisely: ψx, y, z = ηx, y, ze iωx,y,z, where η ηx, y, z cosωx, y, z is the molecular mass density, η the locally uniform mass density and η the mass density of the smectic layers. The smectic layers themselves are the level sets of ω. Therefore, if ψ = the molecules have no positional order, the centers of masses are distributed uniformly, and the liquid crystal is in the nematic phase. The assumption that the phenomenon under investigation is essentially one-dimensional is reflected in the choice the authors make for ψ in 16], they consider ηx, y, z = η constant and ωx, y, z = q z ux, where ux measures the displacement from the unmodulated structure. Rewriting the energy in terms of u, and the components of nx = n 1 x, n x, n 3 x and px = p 1 x, p x, p 3 x yields: fu, n, p = 1 K n n 1 n n 3 ] c η q u n 1 n 3 cos θ B c η q u 1 u n 1 n 3 ] a η q 4 u 1 u n 1 n 3 ] a η q u u n 1 n 3 n 1 ] a η q 4 u n 1 n 3 cos θ B 4 a η q u n 1 n 3 n 1 n 1 cos θ B ] P p 1 4 ɛɛ K p η q u 1 u n 1 n 3 ] p 1 1 K p p 1 p p 3 ] K np η q u p 1 p 3. For a fixed bounded interval = a, b, we let A = { u, n, p : u H, n, p H 1 ; S with n p = in }, and define the free energy: Fu, n, p = fu, n, p dx. Proposition 1. The functional F is bounded below in the class of functions u, n, p A, which verifies either periodic boundary conditions for u, n 1 and p 1, or zero boundary conditions only for p 1. More precisely, for such functions there are positive constants

7 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 7 C and C for which Fu, n, p C u u n p ] dx C. 5 Proof. We start by showing some lower bounds for the terms involving u and its derivatives. Using the simple fact that c d ɛ 1ɛ c ɛ d for any ɛ >, and real numbers c and d, we have for fixed ɛ, ɛ 1 > : a u u n 1 n 3 n 1 ] a u n 1 n 3 n 1 n 1 cos θ B ] min { } ɛ a, a u ɛ1 ɛ u n 1 n 3 n 1 ] 1 ɛ 1 ɛ 1 ] ɛ 1 cos θ B n 1, as well as, for fixed ɛ > : a u 1 u n 1 n 3 ] a ɛ 1 ɛ u 1 ] a ɛ u n 1 n 3 4. Additionally, we note that, for any ɛ 3, ɛ 4 > it holds a u n 1 n 3 cos θ B 4 a ɛ 3 1 ɛ 3 u n 1 n 3 u n 1 n 3 cos θ B a ɛ 3 cos θ B u n 1 n 3 cos θ B a ɛ 3 1 ɛ 3 ɛ 4 1 ɛ 4 u n 1 n 3 4 a ɛ 3 1 ɛ 3 ɛ 4 cos θ B u n 1 n 3 a ɛ 3 cos θ B u n 1 n 3 cos θ B, and, for any ɛ 5, ɛ 6 > : a u u n 1 n 3 n 1 ] a ɛ 5 1 ɛ 5 u n 1 n 1 n 3n 1 n 3 n 1 a ɛ 5 u 1 n 1 4 a ɛ 5 1 ɛ 5 ɛ 6 1 ɛ 6 u n 1 n 1 a ɛ 5 ɛ 6 1 ɛ 5 n 3n 1 n 3 n 1 a ɛ 5 u 1 n 1 4 a ɛ 5 1 ɛ 5 ɛ 6 1 ɛ 6 u n 1 n 1 a ɛ 5 ɛ 6 1 ɛ 5 n 3n 1 a ɛ 5 ɛ 6 1 ɛ 5 n 3 n 1 a ɛ 5 u 1 n 1. The choices ɛ = ɛ1 1ɛ 1 and ɛ j = 1 for j = 3, 4, and the condition n = 1 result in the following lower bound inequality, which holds for any ɛ k >, k =, 5, 6: Fu, n, p Kn 1 min { } a, a η q ɛ 1 cos θ B a η q ɛ 5 ɛ 6 1 ɛ 5 n 1 n n 3 ] c η q u n 1 n 3 cos θ B c η q u 1 u n 1 n 3 ] a η q 4 u 1 u n 1 n 3 ] a η q 4 u n 1 n 3 cos θ B 4

8 8 TIZIANA GIORGI AND FERAS YOUSEF a ɛ η q 4 u 1 ] 1 1 ɛ a 4 a ɛ η q 4 u n 1 n 3 4 a 4 η q 4 cos θ B u n 1 n 3 a η q 4 cos θ B u n 1 n 3 cos θ B a η q u n 1 n 3 n 1 n 1 cos θ B ] a η q ɛ 5 1 ɛ 5 ɛ 6 1 ɛ 6 u n 1 n 1 a η q ɛ 5 u 1 n 1 1 min { a, a } η q ɛ 1 1 ɛ 1 u P p 1 ɛɛ K p η q u 1 u n 1 n 3 ] p 1 ] 1 K p p 1 p p 3 ] K np η q u p 1 p 3 dx. 6 To bound the coupling term we instead need to enforce the boundary conditions, and we obtain for all ɛ j >, j = 7..11: u 1 u n 1 n 3 ] p 1 dx = u b p 1 u u p 1 dx a b p 1 u n 1 b p 1 u n 1 u n 1 p 1 dx a a u n 1 n 3 n 3 p 1 dx ɛ 7 u 1 u p 1 dx ɛ 7 u n 1 u n 1 u n 1p 1 dx ɛ 8 ɛ 9 p 1 dx 1 u n ɛ 1n 3 dx 1 n 3 4 dx, 7 8 ɛ 9 with u n 1 u n 1 u n 1p 1 dx ɛ 1 ɛ 11 u n 1n 1 1 ɛ 11 u dx 1 ɛ 1 u n 4 1p 1 dx u p 1 dx. 8 Recalling that p = n = 1, and using 7 and 8 in equation 6, leads to our next inequality: Fu, n, p Kn 1 min { } a, a η q ɛ 1 cos θ B a η q ɛ 5 ɛ 6 1 ɛ 5 n 1 n n 3 ] c η q u n 1 n 3 cos θ B c η q u 1 u n 1 n 3 ] a η q 4 u 1 u n 1 n 3 ] a η q 4 u n 1 n 3 cos θ B 4 a ɛ η q 4 u 1 ] 1 a 1 ɛ 4 a ɛ η q 4 u n 1 n 3 4 a 4 η q 4 cos θ B u n 1 n 3 a η q 4 cos θ B u n 1 n 3 cos θ B

9 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 9 a η q u n 1 n 3 n 1 n 1 cos θ B ] a η q ɛ 5 ɛ 6 1 ɛ 5 1 ɛ K p η q ɛ 11 u n 1 n min { } a, a η q ɛ 1 1 ɛ K p η q ɛ 7 ɛ 1 a ] 1 η q ɛ 5 u K 1 p η q u K p η q n 3 4 ɛ 7 ɛ 8 ɛ 1 ɛ 11 ɛ 9 P ] 1 p p 1 ɛɛ K p K p η q ɛ 8 ɛ 9 1 p p 3 ] K np η q u p 1 p 3 ] dx, 9 which is verified for any ɛ j >, j = 1, and j = To obtain the wanted lower bound, we start by choosing ɛ 1 small enough to have K n 1 min { } a, a η q ɛ 1 cos θ B > C > for some constant C ; then we select ɛ 5, ɛ 7 and ɛ 1 so to have the coefficient of the u term positive, say greater of a constant C 1 >, and after that we pick ɛ 6 small enough to verify C a η q ɛ 5 ɛ 6 > C >, for some constant C. Once, ɛ 5 and ɛ 6 are fixed, we 1 ɛ 5 can find an ɛ 11 for which the coefficient of u n 1 n 1 is a positive constant, denoted by C 3. Finally, we choose ɛ 8 and ɛ 9 so that the coefficient of p is greater than some constant C 4 >, and ɛ so that the one of u n 1 n 3 4 is greater than some constant C 5 >. If we next set C 6 = K 1 p η q 1 1, C 7 = K p η q ɛ 7 ɛ 8 a ɛ 1 1 ɛ 1 ɛ η q 4, we are left with the inequality: Fu, n, p C n 1 n n 3 ] ɛ 11 ɛ 9, and C 8 = c η q u n 1 n 3 cos θ B c η q u 1 u n 1 n 3 ] a η q 4 u 1 u n 1 n 3 ] a η q 4 u n 1 n 3 cos θ B 4 C 8 u 1 ] C5 u n 1 n 3 4 a 4 η q 4 cos θ B u n 1 n 3 a η q 4 cos θ B u n 1 n 3 cos θ B a η q u n 1 n 3 n 1 n 1 cos θ B ] C 3 u n 1 n 1 C 1 u C 6 u C 4 p 1 p p 3 ] ] P ɛɛ p 1 K np η q u p 1 p 3 dx C 7, 1 for some positive constants C j >, j = We let C = min{c 1, C, C 4, C 8 }, and use in 1 the elementary bound αt βt β α, which holds for every α, β >, to conclude that there exist positive

10 1 TIZIANA GIORGI AND FERAS YOUSEF constants C and C such that Fu, n, p K np η q and the proposition follows. C u p 1 p 3 dx u u n p ] dx C, 11 The previous proposition implies the existence of minimizers of F in the following classes of functions: { } A = u, n, p A with u = p 1 =, A p = {u, n, p A with periodic boundary conditions for u, n 1 and p 1 }. Proposition. There exist minimizers for F in A and A p. Proof. We will first prove the result for A, and then point out the modifications needed for the proof with A p. The energy functional F is bounded below, by Proposition 1, and above, since if we pick n = 1,, and p =, 1, we have that F, n, p = η q ] c cos θ B c a q a q cos 4 θ B. Therefore, we can consider a minimizing sequence {u j, n j, p j } A, namely F u j, n j, p j inf F u, n, p, and from 5 for any c > fixed, we can find u,n,p A a j c such that C u j u j n j p j ] dx C F u j, n j, p j inf F u, n, p c, u,n,p A for any j > j c. From this inequality, since n j L = p j L = 1/, and due to the zero boundary conditions on u j, we conclude that {u j, n j, p j } is uniformly bounded in H H 1 ; R 3 H 1 ; R 3, j > j c. Hence, there exists a subsequence, still denoted by u j, n j, p j, such that u j, n j, p j u, n, p weakly in H H 1 ; R 3 H 1 ; R 3, see 5, Theorem 3.18], strongly in H 1 L ; R 3 L ; R 3, by the Rellich-Kondrachov Compactness Theorem 1], and such that u j, u j, n j, p j all converge almost everywhere in, 1, pg 3]. In fact, since R and bounded, we know that H j1 C j,λ for j =, 1 and < λ < 1 always compact imbedding, so the convergence is uniform, 1, pg 168]. From this, it is straightforward to conclude that u, n, p A. By the weak lower semicontinuity of the L norm, and Proposition 3.5 in 5], since u j 1 u j n j 1 n j 3 ] converges strongly in L, we have F u, n, p lim inf j F u j, n j, p j = inf u,n,p A F u, n, p. For functions in A p, the lower and upper bounds on the energy, which are still valid, allows us to find as before a subsequence whose nematic and polarization directors converge to suitable n and p. Instead, for the u j s we can assume that the subsequence u j converges weakly in H1, strongly in L and uniformly to a function w H 1, such that w = w 1, this last condition due to the periodic boundary conditions on the u j and the uniform convergence in. Hence, the function u x = x w t dt is such that u, n, p A p. The proposition then follows if we notice that w j x = u j x u j converges weakly

11 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 11 in H, strongly in H 1 and uniformly to u, and verifies F w j, n j, p j = F u j, n j, p j. 4. Columnar phases simplified model Model and notations. The so-called B 1 phase observed in bent-core molecule liquid crystals is characterized by a two-dimensional structure, where columns of broken smectic layers result in an antiferroelectric structure. In this phase the spontaneous polarization is perpendicular to the column axis, if instead the polarization is along the column axis, the material is in a so-called reversed B 1 phase, denoted by B 1rev or B 1rev tilted, depending if the nematic director is parallel or tilted with respect to the normal of the smectic layers, 7, 8]. The experimental work presented in 8] deals with both reversed B 1 phases, since in these phases the spontaneous polarization can be reoriented using an electric field, without a complete rebuilding of the underlining columnar structure. The modeling and numerical part is dedicated to the study of the B 1rev tilted phase, where the tilt allows the molecules to switch polarization either by rotating around the molecular axis or around the cone, or more in general by a combination of both, see Figure 4. Figure 4. Schematic representation of the B 1rev tilted phase structure, and switching by an external electric field In particular, the authors are interested in the dependence on physical parameters of the type of switching, that is in knowing if the rotation is along the axis of the molecules or the cone for particular materials and sample size, since the two different rotations have different consequences on the structural chirality. They adopt an energy, inspired by the density, which includes a term that accounts for the interaction of the external electric field E and the polarization. They observe that in the B 1rev tilted the layer normal and the tilt angle θ are constant, so that θ = θ B or π, and all the terms involving only the smectic order parameter are

12 1 TIZIANA GIORGI AND FERAS YOUSEF constant. In here, to represent this situation, we consider the energy density: 1 K n n n ] P p 1 ɛɛ ˆK p p n ν 1 K p p p ] ˆK np p n ν P p E, 1 where ν is the layer normal. We non-dimensionalize the energy by rescaling distances with respect to the column width L, and multiply it by the factor ɛɛ. As in 8], we also consider a reference frame with the y-axis parallel to the columnar axis, so that ν = e 3 <,, 1 >, and take the applied electric field in the same direction as the columnar axis, that is E = E e. In conclusion, using the fact that n ν = sin θ B, we are lead to work with the following non-dimensionalized energy functional: E L n, p = e L n, p dx, 13 where D =, 1, and with D e L n, p = 1 k n n L 1 n n 3 ] 1 kp p L 1 p p 3 ] L p 1 ˆk p sin θ B p 1 1 np p 3 1 ] L k E p 1 p sgne p 3, kn = ɛɛ P K n ; ˆk p = ɛɛ P P ˆK p ; kp = ɛɛ P K p ; k np = ˆK ɛɛ np P ; k E = ɛɛ E. 14 P To mimic the switching mechanism, we consider the boundary conditions: p = p1 =<, 1, > and n = n1 =< sin θ B,, cos θ B >, 15 and take E > so that the applied electric field is in the opposite direction as the spontaneous polarization at the boundary. This boundary conditions reflect the observation made in 8] that in a sample only the molecules of every other column will rotate, which can be modeled by considering only one single column with strong anchoring of the directors at the boundary. In particular, this implies the use of zero boundary conditions for the first component of the polarization director p. Under these assumptions the energy can be written as: E L n, p = D 1 k n n L 1 n n 3 ] 1 kp L L p 1 1 np p 3 1 L k E p 1 p p 3 ] p 1 p 1 p 3] ] dx. 16 In what follows, we study how the relative size of the physical parameters in the energy density influence the minimizing states in the case of large column width. More precisely, we consider the Γ-limit of the energy for L tending to infinity, and study the minimizers of the limiting energy.

13 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 13 Remark 1. While considering the limiting behavior as the width of the column tends to infinity is consistent with the discussion in 8], one could take the point of view of keeping the width fixed, and let the elasticity constants become small relatively to the electric and coupling constants. In this interpretation, one could still work with the family of energies 16, where now L would be the mathematical parameter expressing the relative behavior of the different constants. A mathematical analysis similar to the one below could then be applied. We recall that given a metric space X and a family of functionals F L : X, ], parameterized by L, we say that F L Γ-converges to a functional F : X, ] on X as L if the following two conditions are satisfied: i Liminf inequality: For every x X, and every sequence {x L } X converging to x, we have lim inf L F Lx L F x. ii Limsup inequality: For every x X, there exists a sequence {x L } X converging to x for which lim sup F L x L F x. L In this paper, to study the behavior of the energy minimizers for the large domain case, we consider the space: X = { n, p : D S measurable functions, with n p =, n 3 = cos θ B, n = n1 =< sin θ B,, cos θ B >, p = p1 =<, 1, >}. As expected, the minimizers of the limiting energy, obtained considering the L 1 metric, belong to the the following subspace of X: Y = {n, p Y such that Sp = } X, where Sv = number of points of discontinuity of v, and with Y = {n, p such that p BV l D; { 1, 1}, p = 1 in D, n 1 BV r D; { sin θ B, sin θ B }, n 3 = cos θ B, p 1 = p 3 = n = }, BV l D; { 1, 1} = {v BV, 1 with values in { 1, 1}, and such that v = v1 = 1}. We recall that functions in BV, 1 are continuous, with possibly a finite number of jump discontinuities. Thus, in the above, if x = or x = 1 are points of discontinuity of v, the boundary conditions are intended as v = v1 = 1, in other words, BV l D; { 1, 1} can be identified with the space of functions v BV loc R with values in { 1, 1}, such that v 1 on, 1,, 14]. The space BV r D; { sin θ B, sin θ B } is defined in an analogous way, with now boundary conditions v = v1 = sin θ B.

14 14 TIZIANA GIORGI AND FERAS YOUSEF The expression of the limiting energy for functions n, p Y is given by: { } { T } T Fn, p = inf inf N Ew, z dx inf inf Ev, u dx T T inf T inf { T Ev 1, u 1 dx where N= total number of points of discontinuity of n 1 and p counted once, and the first infimum is taken over functions w, z, u, v, u 1, v 1 H 1, T ; S, such that w±t = ± sin θ B,, cos θ B ; v±t = n 1 ±,, cos θ B ; v 1 ±T = n 1 1 ±,, cos θ B ; z±t =, 1, ; u±t =, ±1, ; u 1 ±T =, 1, ; } w z = ; v u = ; v 1 u 1 = ; w 3 = cos θ B ; v 3 = cos θ B ; v 1 3 = cos θ B ; and Ea, b = 1 k n a kp b b 1 1 k np b 3 1 k E b 1 b 1 b ] 3. Finally, we extend the energy functional in a standard fashion so that it will be defined for all functions in X: E L n, p if n, p H 1 D; S H 1 D; S X F L n, p = if n, p X \ H 1 D; S H 1 D; S 17, and similarly for the limiting functional: Fn, p if n, p Y F n, p = if n, p X \ Y. 4.. Compactness and Γ-convergence. We start by deriving a compactness result. Proposition 3 Compactness. Let and {n k, p k } X such that F Lk n k, p k is uniformly bounded, then {n k, p k } is pre-compact in X, that is there exists n, p Y X and a subsequence { n kj, p kj } such that nkj, p kj n, p in the L 1 -norm. Proof. Assume there exists M > independent of k such that F Lk n k, p k M, then by definition of F Lk we deduce that n k, p k H 1 D; S H 1 D; S. Additionally, if we denote by p j k and n j k, j = 1,, 3, the components of the vectors p k and n k, we see that 1 p 1 k dx M, 1 p 3 M k dx, k np, 1 p k 1] dx M k E, hence, as, we have p 1 k, p k 1 and p 3 k in L D, and D bounded implies that they also converge in L 1 D and a.e. up to subsequences. 18

15 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 15 Looking at boundary condition, since : p k has a continuous representative, we can extend it to the whole real line as follows: { p p k = k, if x D 1, if x R \ D, so that p k L 1 loc R, and p k { 1, if x D 1, if x R \ D, pointwise a.e. up to subsequences, and in L 1 loc R. In other words, p k p BV l D; { 1, 1} in the L 1 -norm. Again from the uniform bound for the energies of the elements of the sequence, since p 3 k = p k n k ν = n 1 n n1 p n p 1, 19 we have, where we drop the subscript k for simplicity of notation, that M 1 Lk np = 1 Lk np = 1 Lk np = 1 Lk np = 1 Lk np from which we gather 1 Lk np M 1 Lk np M 1 Lk np n 1 n n1 p n p 1 ] dx sin θ B n 1p n p 1 n 1 p 1 n p dx sin θ B n 1 n 1 1 p n p ] 1 n 1 p 1 n p dx sin θ B n 1 n 1 1 p n p 1 n 1 p 1 n 3 p 3 n 1 p 1 ] dx sin θ B n 1 n 1 1 p n p 1 n 1p 1 n 3 p 3 n 1 p 1 ] dx, sin θ B n ] 1 n 1 1 p dx 1 1 n p 1dx 1 Lk np p 1dx 1 Lk np 1 1 n 1p 1dx 1 Lk np 3 n 1 p 1 n 3p 3 dx 1 n 3 p 3 n 1 p 1 dx

16 16 TIZIANA GIORGI AND FERAS YOUSEF M 1 k np L 1 p 1dx 3 k np L 1 p 1dx 1 Lk np 1 p 3dx k np M M. Equation allows us to follow the classical Modica-Mortola argument for the term involving n 1. Namely, we start by noticing that knp k n 1 then we define φ 1 t = sin θ B n 1 n 1 dx 1 Lk np 1 kn L t sin θ B 1 1 sin θ B n 1 dx n 1 dx M k np M M, 1 sin θ B s ds, so that the change of variable t = n 1 x, and the chain rule in W 1,1 D, here φ 1 C 1 sin θ B, sin θ B ], give 3M k np M k np k n = k np k n 1 1 sin θ B n 1 x n 1x dx φ 1 n 1 x] dx, that is, φ 1 n 1 x] is uniformly bounded in L 1 D. On the other hand, since n 1 x sin θ B, we have φ 1 n 1 x sin θ B, which implies φ 1 n 1 x is uniformly bounded in L 1 D. In conclusion, reintroducing the subscript k, the sequence of functions W k := φ 1 n 1 k is uniformly bounded in W 1,1 D. By the Rellich Kondrachov theorem, there exists a subsequence of {n 1 k }, which we still label with k, and a function W BV D such that W k := φ 1 n 1 k W in L 1 loc D, also by taking a further subsequence, we can assume W kx W x a.e. in D. From, the function wn 1 k x = sin θ B n 1 k x is such that 1 wn1 k x 1 dx = sin θ B n 1 kx dx 4 M k np M, and, for we have wn 1 k in L 1 D, a.e. up to a subsequence. φ 1 t is continuous and strictly increasing in sin θ B, sin θ B ], hence φ 1 1 is continuous and n 1 k x = φ 1 1 W kx φ 1 1 W x a.e. We set n 1x := φ 1 1 W x, and remark that since wn 1 k x a.e. in D it must hold wn 1 x = a.e. in D, which gives n 1 x { sin θ B, sin θ B } a.e. in D. In turn, W x {φ 1 sin θ B, φ 1 sin θ B }, but φ 1 sin θ B =, so there exists a set E D such that W x = φ 1 sin θ B χ E. Now, W BV D and D has finite measure thus χ E BV D, and n 1 = sin θ B χ E sin θ B 1 χ E. From this, we conclude that n 1 k converges a.e. in D to a BV function with values in { sin θ B, sin θ B }. Furthermore, if outside D we extend n 1 k to be sin θ B, we

17 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 17 obtain that n 1 k sin θ B a.e. in R\D, we can then conclude that n 1 k converges a.e. to a BV r D; { sin θ B, sin θ B } function. Convergence of n 1 k in L 1 loc R follows by Egoroff s and Vitali s convergence theorems if we show {n k } is equi-integrable, i.e. for all ε > there exists δ > such that n k dµ ε for all k and F D with µf δ, but since n k = 1, F n k dµ µf, we can pick δ = ε. F Finally, using the fact that n 1 k n k = sin θ B we also have that n k converges to in L 1 and a.e. We are able to prove the following standard Γ-convergence result: Theorem 4.1. Consider X as a metric space endowed with the L 1 -norm, then for L the functional F L Γ-converges to F. Proof. Liminf Inequality: Consider n, p X and {n k, p k } X such that n k, p k n, p in L 1, we need to show that lim inf F n k, p k F n, p. If lim inf F n k, p k = there is nothing to prove, otherwise there exists a subsequence such that lim F L nkj kj, p kj = lim inf F L nkj kj, p kj <, and for j j j large enough: < F Lkj nkj, p kj < lim inf j F j nkj, p kj 1 M, which implies n kj, p kj H 1 D; S H 1 D; S. Proceeding as in the proof of Proposition 3, we can find another subsequence n kjl, p kjl which converges in L 1 to an element n, p in Y. By uniqueness of limit, we conclude that n, p n, p Y and lim F jl n kjl, p kjl = lim inf F L nkj kj, p kj. l j Therefore, we have reduced our task to prove that if n, p Y and {n k, p k } H 1 D; S H 1 D; S X is such that n k, p k n, p in L 1, and < F Lk n k, p k < M < then lim inf F n k, p k F n, p. Given n, p Y, we denote by N the total number of jump discontinuities of n 1 and p counted once. Since Sp = for elements in Y, we know that N, while from n 1 BV r D; { sin θ B, sin θ B } we have N <. Let = θ 1 <... < θ j <... < θ N = 1 be the points of discontinuity, we observe that if 1 < j < N only n 1 has a jump at θ j, and at j = 1, N p 1 is discontinuous while n 1 might or might not be discontinuous. By passing to a subsequence if necessary, we can assume pointwise convergence a.e., and for any η, sin θ B we can find δj 1, δ j such that θ j δj 1, θ j δj are points of convergence for n k and p k, and the intervals I j = θ j δj 1, θ j δj are disjoint. Then, since N <, for k large enough, say k > k η, we can assume that the following relations hold for the components of n k and p k. For n k : sin θ B η < n 1 k θj δ 1 j sin θb sin θ B n 1 k θj δ j sin θb η 3

18 18 TIZIANA GIORGI AND FERAS YOUSEF or if 1 < j < N, and if j = 1, and sin θ B n 1 k θj δ 1 j sin θb η sin θ B η n 1 k θj δ j sin θb, 4 n 1 k θ1 δ 1 1 = sin θb n 1 n 1 k θ1 δ 1 n1 η, 5 n 1 k θn δ N = sin θb n 1 1 n 1 k θn δ 1 N n1 1 η, 6 for j = N. Which also imply, using the constraints n k = 1, and n 3 k = cos θ B : with η as η. For p k : if 1 < j < N, and if j = 1, and η n k θj δ 1 j η, η n k θj δ j η, 7 η < p 1 k θj δ 1 j < η, η < p1 k θj δ j < η, η < p 3 k θj δ 1 j < η, η < p3 k θj δ j < η, 1 η < p k θj δ 1 j < 1, 1 η < p k θj δ j < 1, 8 p 1 k θ1 δ 1 1 =, η < p1 k θ1 δ 1 < η, p 3 k θ1 δ 1 1 =, η < p3 k θ1 δ 1 < η, p k θ1 δ 1 1 = 1, 1 η < p k θ1 δ 1 < 1, 9 η < p 1 k θn δ 1 N < η, p1 k θn δ N =, η < p 3 k θn δ 1 N < η, p3 k θn δ N =, 1 η < p k θn δ 1 N < 1, p k θn δ N = 1, 3 if j = N. Therefore, up to subsequences, for k > k η, starting from N 1 1 F Lk n k, p k 1 L j= k k n n k 1 dx k p p k dx I j I j p 1 1dx k np p 1 3dx k E p 1 p 1 p ] 3 dx I j I j I j

19 1 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 19 1 k n θ 1δ 1 θ 1δ 1 θ 1 δ 1 1 n k dx 1 k p θ 1δ 1 θ 1 δ 1 1 p k dx p 1 1 k npp 3 1 k E p 1 p 1 p 3] dx 1 θ 1 δ k n θ N δ N n k dx 1 k p θ N δ N p k dx θ N δ N θ N δ 1 N θ N δ 1 N θ N δ 1 N p 1 1 k npp 3 1 k E p 1 p 1 p 3] dx, and denoting T j = δ j, where δ j = δ j δ1 j, the change of variable x = 1 s θ j δ j δ1 j gives: F Lk n k, p k N 1 j= T j j E w j,k, z j,k dx T N T 1 1 E w 1,k, z 1,k dx N E w N,k, z N,k dx, 31 with w j,k, z j,k H 1 j, T j ; S, w j,k z j,k =, w j,k 3 = cos θ B, and where at ±T j the functions w j,k, z j,k verify conditions 3 or 4, and 7, 8 if 1 < j < N, while w 1,k, z 1,k verify conditions 5, 7, and 9, and w N,k, z N,k verify conditions 6, 7, and 3. Taking the infimum over all functions w j, z j, v, u, v 1, u 1 verifying the assumptions above, for k > k η we derive: F Lk n k, p k inf w j, z j v, u v 1, u 1 { N 1 j= T j j E w j, z j dx T N } E v 1, u 1 dx. T 1 1 E v, u dx N Additionally, since the limit function p is constant in I j in both 3 and 4, by symmetry we can fix one of the two conditions, say 4, and conclude, for all k > k η, that: { N 1 F Lk n k, p k inf j= T j j E w, z dx T 1 1 E v, u dx

20 TIZIANA GIORGI AND FERAS YOUSEF T N N { T inf inf T } { T E v 1, u 1 dx inf inf N T } { T E v, u dx inf inf T } E w, z dx } E v 1, u 1 dx, 3 where the first infimum is taken over functions in H 1, T ; S, with w z = v u = v 1 u 1 =, w 3 = v 3 = v 1 3 = cos θ B, and where at ±T the functions w, z verify condition 4, 7, and 8, while v, u verify conditions 5, 7, and 9, and v 1, u 1 verify conditions 6, 7, and 3, all of which we remind the reader depend on η, letting η we obtain: lim inf F n k, p k F n, p. Limsup Inequality: We need to show that for every n, p X there exists {n k, p k } X such that n k, p k n, p in L 1 and lim supf Lk n k, p k F n, p. Let n, p X, if n, p X \ Y then F n, p =, and we set n k, p k = n, p for all k. If n, p Y instead, then F n, p <, that is F n, p <. Let N be the total number of points of discontinuity of n 1 and p counted once. By the definition of F n, p, for any k >, there exists T 1, T, T 3 and w, z H 1 1, T 1 ; S, v, u H 1, T ; S and v 1, u 1 H 1 3, T 3 ; S, with suitable boundary conditions, such that N T 1 1 E w, z dx T E v, u dx T 3 3 E v 1, u 1 dx F n, p We let = θ 1 < θ <... < θ N = 1 denote the points of discontinuity. Note that T 1, T, T 3 and all the functions above depend on the choice of k. We next pick > k large enough so to have: < T and define: n k x = < θ T1 < θ T1 v x T w ± x θ j v 1 x T3 1 nx < θ 3 T1 if θ j T1 <... < θ N 1 T1 1 k. < 1 T3 < 1, if x T x θ j T1, j =,, N 1 if 1 T3 x 1 otherwise. 33 We also consider an analogous definition for p k x. As in ], we remark that the choice between the plus and minus sign in 33 can be made in such a way that the resulting function is continuous, and the sequence n k, p k is an element of X. It is also a consequence of its definition that the sequence converges in L 1 to n, p.

21 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 1 The energy of an element n k, p k of the sequence is F Lk n k, p k = T k n n k 1 kp p k p k 1 1 k np p k 3 1 k E p k 1 p k 1 p k 3] dx N 1 j= θ j T 1 θ j T 1 k n n k 1 kp p k p k 1 1 k np p k 3 1 k E p k 1 p k 1 p k 3] dx 1 1 T 3 k n n k 1 kp p k p k 1 1 k np p k 3 1 k E p k 1 p k 1 p k 3] dx, where the other parts of the integral are zero since n k = n = constant, p k = p = constant, and p 1 = p 3 =, p = 1 between singularities. If we now make the changes of variables: s = x T on, T ], we obtain: s = ± x θ j s = x T3 1 F Lk n k, p k = N T 1 F n, p 1 k, 1 E w, z ds on θ j T1, θ j T1 ], ] on 1 T3, 1, T E v, u ds T 3 3 E v 1, u 1 ds and by taking the lim sup for k, recalling that we picked > k, we gather lim supf Lk n k, p k F n, p. 5. Minimizers of Γ-limit. In this section, we study the minimizers of the limiting functional F in 18, here we are interested in their behavior as functions of the physical parameters. We remark that existence of minimizers for this functional is a consequence of the Γ-limit procedure, since the Γ-limit as defined in the previous section is always lower semicontinuous in X. The energy F of an element n, p in X is unbounded unless n, p Y, therefore we can restrict our attention to elements in Y. If, as mentioned in the previous section, we think of BV r D; { sin θ B, sin θ B } as the space of functions v in BV loc R such that v sin θ B on, 1,, for n, p Y we have that

22 TIZIANA GIORGI AND FERAS YOUSEF n 1 has one of four fundamental behaviors described for simplicity via the graphs presented in Figure 5. Figure 5. Possible behaviors near the boundary of the component n 1 of n, p Y. As a consequence, we claim that there are only two candidates for minimizers of F, namely, n c, p d and n d, p d, where and n c x = sin θ B,, cos θ B if x, 1, n d x = sin θ B,, cos θ B if x, 1, p d x =, 1, if x, 1, 34 n c = n c 1 = n d = n d 1 = sin θ B,, cos θ B, p d = p d 1 =, 1,. 35 This can be justified, by noticing that if n, p n d, p d is of Type I then F n d, p d < F n, p, while if n, p n c, p d is of Type II then F n c, p d < F n, p, and, finally, that if n, p is of Type III and ñ, p is of Type IV, then min {F n c, p d, F n d, p d } < min {F n, p, F ñ, p}. These inequalities are all true since for n c, p d and n d, p d one has N =, and the last two infima in the definition of F depend only on the behavior of the functions near the boundary. To determine the minimizers of F, we then need to compare the energies of n c, p d and n d, p d, which for these states simplify to: T T inf inf E v, u dx T inf inf E v T 1, u 1 dx, 36

23 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 3 for v, u, v 1, u 1 H 1, T ; S such that v±t = n 1 ±,, cos θ B, v 1 ±T = n 1 1 ±,, cos θ B, u±t =, ±1,, u 1 ±T =, 1,, v u = = v 1 u 1, v 3 = v 1 3 = cos θ B. 37 Furthermore, if we fix either n c, p d or n d, p d, it s straightforward to see that given v, u verifying the relevant conditions in 37, if we consider v 1 s := v s and u 1 s := u s, then v 1, u 1 will verify the analogous conditions in 37, and E v, u = E v 1, u 1. In a similar way, given a v 1, u 1 we can build an admissible v, u with same energy E, and we can conclude that the two infima in equation 36 are equal. A summary of our discussion can be express in the { form of the energies } for our two candidates minimizers: F n c, p d = inf T inf H 1, T ; S with T v±t = sin θ B,, cos θ B, u±t =, ±1,, E v, u dx, for v, u v u =, v 3 = cos θ B, 38 { } T and, F n d, p d = inf inf E v, u dx, for v, u H 1, T ; S with T v±t = ± sin θ B,, cos θ B, u±t =, ±1,, v u =, v 3 = cos θ B. 39 At first sight one might think that the state with continuous n 1 component, that is n c, p d, might provide the least energy, but a careful examination of the terms involved shows that the relative size of the physical parameters will play a role, as displayed in the following result. Proposition 4. For fixed k E, k n, and k p, if k n k p there exists c > such that for any k np > c the global minimizer of F is n d, p d. Proof. We rewrite the energy in terms of the angles introduced in equation 1, where we notice that in here we are considering = and θ = θ B, so that if v = u = 1, v u =, and v 3 = cos θ B, we have: v =< sin θ B cos φ, sin θ B sin φ, cos θ B >; u =< sin α cos θ B cos φ sin φ cos α, sin α cos θ B sin φ cos φ cos α, and, we can rewrite T sin α sin θ B >, E v, u dx = T 1 k n sin θ B kp cos θ B kp sin θ B cos α φ

24 4 TIZIANA GIORGI AND FERAS YOUSEF and 1 k p α kp cos θ B α φ 1 k E 1 k np k E sin α sin θ B sin α cos θ B cos φ sin φ cos α k E sin α cos θ B sin φ cos φ cos α 1 The boundary conditions in 38 and 39 translate to α =, αt = π, ] dx. 4 φ = = φt, 41 α = = αt, respectively. We use the assumption k n k p to rewrite 4: T E v, u dx = k p T cos θ B α φ 1 k E φ =, φt = π, 4 1 k p sin θ B α 1 k np k E sin α sin θ B 1 k n kp sin θ B 1 ] k p sin θ B cos α φ sin α cos θ B cos φ sin φ cos α k E sin α cos θ B sin φ cos φ cos α 1 ] dx. 43 and by a standard procedure, we find a lower bound for the energy of n c, p d : F n c, p d T 1 k p sin θ B α 1 k np k ] E sin α sin θ B dx sin θ B k p k np k E T sin α α dx = 4 sin θ B k p k np k E, 44 where we used in an essential way the boundary condition 41. To find an upper bound for F n d, p d, we just consider the linear interpolation between the boundary values, and pick αx, φx = π T x π.

25 BENT-CORE LIQUID CRYSTALS COLUMNAR PHASES 5 This choice yields: F n d, p d inf T T > 1 k E { π 1 = inf T > T { π = inf T > 4 1 k n sin θ B kp π T sin π T x π k E k n sin θ B kp ] π cos 1 k E k n sin θ B k p 1 T 1 k E T T x π T k E }, ] 1 dx } 3T which, taking the infimum over T >, results in to the bound: π F n d, p d 4 k 4 n sin θ B kp 1 ke. 45 Comparing 44 and 45, we have that F n d, p d < F n c, p d, for k np large enough. The discussion proceeding the statement of this proposition implies that n d, p d is the global minimizer for such k np values. Remark. As a final remark, we should explain why we claim that the state n c, p d represents rotation around the molecular axis, while n d, p d suggests rotation around the tilt cone. Looking at the boundary conditions 38, expressed as 41 via the polar angle α and the azimuthal angle φ, we see that α is changing from to π, while φ is fixed at, that is the molecule is rotating around its axis while preserving its position on the tilt cone. Similarly, looking at the boundary conditions 39 in the energy of n d, p d as given in 4, we see that now is the azimuthal angle which is changing from to π, while the polar angle is fixed at, suggesting that the molecule is rotating around the tilt cone, while the relative position of n and p remains unchanged. Acknowledgments. The first author would like to acknowledge the Isaac Newton Institute for Mathematical Sciences, where this work was initiated, for its hospitality and support during her Spring 13 visit. Special thanks to Professor Martin Čopič for kindly explaining some aspects of the model related to the physics of the problem. The authors would also like to thank the referees for their useful suggestions. REFERENCES 1] R. Adams and J. Fournier, Sobolev Spaces, nd edition, Academic Press, 3. ] N. Ansini, A. Braides and V. Valente, Multiscale analysis by Γ-convergence of a onedimensional nonlocal functional related to a shell-membrane transition, SIAM J. Math. Anal., 38 6, ] P. Bauman and Phillips, Analysis and stability of bent-core liquid crystal fibers, Discrete Contin. Dyn. Syst. Ser. B, 17 1, ] P. Bauman, D. Phillips and J. Park. Existence of solutions to boundary value problems for smectic liquid crystals, Discrete Contin. Dyn. Syst. Ser. S, 8 15, ] H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Universitext, Springer, 11. 6] J.-H. Chen and T. C. Lubensky, Landau-Ginzburg mean-field theory for the nematic to smectic-c and nematic to smectic-a phase transitions, Phys. Rev. A, , ] A. Eremin and A. Jákli, Polar bent-shape liquid crystals - from molecular bend to layer splay and chirality, Soft Matter, 9 13,

26 6 TIZIANA GIORGI AND FERAS YOUSEF 8] E. Gorecka, N. Vaupotič, D. Pociecha, M. Čepič and J. Mieczkowski, Switching mechanism in polar columnar mesophases made of bent-core molecules, ChemPhysChem, 6 5, ] S. Joo and D. Phillips, Chiral nematic toward smectic liquid crystals, Comm. Math. Phys., 69 7, ] S. T. Lagerwall, Ferroelectric and antiferroelectric liquid crystals, Encyclopedia of Materials: Science and Technology, 1, ] I. Luk yanchuk, Phase transition between the cholesteric and twist grain boundary C phases, Phys. Rev. E, , ] I. Mu sevi c, R. Blinc and B. Zek s, The Physics of Ferroelectric and Antiferroelectric Liquid Crystals, World Scientific Publishing Company,. 13] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals, nd edition, Clarendon Press, Oxford, ] L. Scardia, A. Schlömerkemper and C. Zanini, Boundary layer energies for non convex discrete systems, Math. Models Methods Appl. Sci., 1 11, ] I. W. Stewart, The Static and Dynamic Continuum theory of Liquid Crystals, Taylor & Francis, 4. 16] N. Vaupotič and M. Čopič, Polarization modulation instability in liquid crystals with spontaneous chiral symmetry breaking, Phys. Rev. E, 7 5, ] E. G. Virga. Variational Theories for Liquid Crystals, Chapman & Hall, London, Received April 14; revised February 15. address: tgiorgi@nmsu.edu address: faamy@nmsu.edu