Role of Lifshitz Invariants in Liquid Crystals

The interaction between an external action and the order parameter, via a dependence described by a so-called Lifshitz invariant, is very important to determine the final configuration of liquid crystal cells. The external action can be an electric field applied to the bulk or the confinement due to free surfaces or cell walls. The Lifshitz invariant includes the order parameter in the form of an elastic strain. This coupling between elastic strains and fields, inserted in a Landau-Ginzburg formalism, is well known and gives rise to striction effects causing undulations in the director configuration. We want to discuss here the role of Lifshitz coupling terms, following an approach similar to that introduced by Dzyaloshinskii for magnetic materials. Case studies on nematics in planar and cylindrical cells are also proposed.


Introduction
The contributions to the free-energy density of terms in the derivatives of order parameter are of great importance and recognised to be fundamental in governing the appearance of spatially modulated structures in magnetic materials and of periodic patterns in liquid crystals. It is possible to identify the same structure in the free energy, when it is represented by a Landau-Lifshitz phenomenological theory of phase transitions; this structure has the form of an invariant term, so-called Lifshitz invariant, which is linear with respect to the gradient of order parameter. As shown by Landau and Lifshitz [1], a system near its phase transition point may be unstable with respect to distortions of the appropriate OPEN ACCESS order parameter. This instability may develop, when the irreducible representation allows a quadratic antisymmetric combination, linear in the order parameter components and in their gradients.
Phases with large-scale space fluctuations of the order parameter were discovered experimentally in the 1960's [2]. Using the approach proposed by Lifshitz, Dzyaloshinskii [3] showed that these configurations are associated with the development of instabilities and found the corresponding approximate solutions of the phase equations. Presently, the family of experimentally observed modulated states has grown both in magnetic and liquid crystal systems [4][5][6][7].
The aim of this paper is to discuss those properties and behaviours of liquid crystal materials originated by Lifshitz contributes to the free energy. Before the discussion of Lifshitz invariants in liquid crystals, we prefer to devote a section of the paper to a brief remark on the use of these invariants in magnetic systems. After this remark, we show how the flexoelectric effect, the chiral elastic term and the saddle-splay surface contribution can be described as Lifshitz contributions.
A liquid-crystal material, the free-energy of which contains a Lifshitz term coupling elastic strains and external fields, can exhibit undulations in the director configuration. Periodic structures in liquid crystal materials can be achieved either in cholesteric or ferroelectric liquid crystals, which possess a natural periodic helicoidal distribution of the molecular orientation. In a nematic liquid crystal cell, a periodic structure can appear spontaneously too, with period that can be controlled by external factors such as applied fields and asymmetric anchoring conditions. The electric field controls the instability produced by the flexoelectric effect. Flexoelectric domains were first observed by Vistin and theoretically studied by Bobylev and Pikin [8,9]. More recently, Lavrentovich and Pergamenshchick discovered another interesting instability in nematics, controlled by the saddle-splay surface contribution to the free energy [10,11].
In the final part of this paper, we propose a detailed discussion of some case studies involving flexoelectricity. In particular, we discuss the hybrid nematic cell in planar geometry and show its complete phase diagrams. The same we shall do for the cylindrical confinement of nematics. To the author's knowledge, the problem of flexoelectricity in cylindrical confinement has not been discussed before: in the framework of its approximate solution, the corresponding phase diagram shows instabilities. The discussion on saddle-splay instabilities concludes the paper.

The Dzyaloshinskii-Moriya Coupling
Some magnetic structures are characterised by a modulation of the spin arrangements over periods, which are long compared to the size of the lattice cell and usually not commensurate with it. The existence of such magnetic structures can be due to competition between exchange interactions or to relativistic effects like spin-orbit coupling. Relativistic interactions were first considered by Dzyaloshinskii [3] and received a microscopic description by Moriya [12]. The Dzyaloshinskii-Moriya (DM) interaction can be written as a product of three vectors, whereas the strength of the spin-orbit coupling gives its intensity [12,13]. The macroscopic manifestation of the antisymmetric DM couplings takes place in noncentrosymmetric magnetic crystals. Dzyaloshinskii showed that, in this case, the DM interaction stabilises long-periodic spatially modulated structures of the vectors i S r , structures with a fixed sense of rotation. In antiferromagnets, the DM-interaction favours arrangements of the magnetic moments, which result in a weak spontaneous magnetisation. Within a continuum approximation for magnetic properties, the interactions responsible for these modulations are expressed by inhomogeneous invariants. In Ref. [14], these contributions to the free magnetic energy, involving first derivatives of magnetisation with respect to spatial coordinates, are defined as the inhomogeneous Dzyaloshinskii-Moriya interactions. These interactions are linear with respect to first spatial derivatives of magnetisation M r in an antisymmetric mathematical form, firstly studied in the theory of phase transitions by E. M. Lifshitz and known as Lifshitz invariants. Spiral structures arise in magnetic systems from the presence of the Lifshitz invariant in the free energy [15].
The structure of the Lifshitz invariant is, in the case of the inhomogeneous Dzyaloshinskii-Moriya interaction, a product of three vectors: a vector D r representing an internal or external field or a fixed direction in the space, a vector M r representing the local order parameter and the ∇ r operator on the order parameter components. The product has the following form: In the case of the liquid crystals, we shall see that vector D r can be an external electric field or the direction perpendicular to the sample surface. It is better to remark that in Ref. 4, we can find another choice for the DM interaction, as the pseudoscalar We will discuss this form in the Sect.5, concerning the chiral nematics. We used the DM interactions in 1996 to study the field-induced phase transition of BiFeO 3 [16]. More recently, the coupling of spin waves with the optical phonons has been discussed in the framework of Lifshitz invariant, for the same material [17]. An antiferromagnetic vector L r characterises the BiFeO 3 spin structure. The Landau-Ginzburg energy density [3] of the spin structure is the following sum of four terms: The first term L f in Equation 2 is the magneto-electric coupling as a Lifshitz invariant, where S , z P is the z-component of the spontaneous polarization vector, and α is the inhomogeneous relativistic exchange constant (inhomogeneous magneto-electric constant). The Lifshitz invariant is the responsible for the spatially modulated spin structure in BiFeO 3 , as shown in Ref. 16 in the form: as a scalar product of two fields. In our paper [16], we investigated the influence of an electric field on the spatially modulated spin structure (SDW state). The electric field has a tendency to prefer a homogeneous state and to induce a phase transition to this state. In that paper, we used the analogy with nematic liquid crystals to study magnetic materials. Here, we want to enhance the analogy of liquid crystal interactions with the two form of the Dzyaloshinskii-Moriya DM interaction.

The Flexoelectricity in Liquid Crystals
Let us consider a nematic liquid crystal and assume as order parameter the director field n r , describing the local mean orientation of molecules. This is usually a unit vector. Vector A r can be used in nematics too, rewritten in the following form: Vector A r in Equation 5 is well known in the physics of liquid crystals. A r is encountered in the structure of flexoelectric contribution to bulk free energy as Flexoelectricity is a property of liquid crystals similar to the piezoelectric effect. In certain anisotropic materials, which contain molecular asymmetry or quadrupolar ordering with permanent molecular dipoles, an applied electric field may induce an orientational distortion. Conversely any distortion will induce a macroscopic polarization within the material. The polarization vector P r in the flexoelectric term is then described with a distortion in the nematic director field: The two terms in the polarization vector are due to the splay and the bend contribution. The coupling of the polarization P r with an external electric field results in the appearance of a periodic distortion of an initial planar orientation of the nematic cell [18]. Meyer showed that the infinite liquid crystal must be disturbed, the perturbation is periodic along the director orientation and the period is inversely proportional to electric field strength [19]. This is not surprising because the polarization vector P r has the same structure of vector A r in Equation 5. In flexoelectricity, the polarization is induced by a deformation of the director field. Let us remember that in the piezoelectric materials, an applied uniform strain can induce an electric polarization or vice versa. Crystallographic considerations restrict this property to noncentrosymmetric systems. A strain, which is not uniform, can potentially break the inversion symmetry and induce polarization in non-piezoelectric materials. While the conventional piezoelectric property is different from zero only for certain select materials, the non-local coupling of strain and polarization could be potentially found in all dielectrics [20]. As a result, we find that the coupling with an external field gives the Lifshitz invariant as a DM non homogenous coupling for the electric field with the Lifshitz vector.

Periodic Distortions in Nematics
Let us discuss more deeply the Meyer result [18,19] of a periodic distortion in the infinite medium. The free energy density is given by: in the uniform elastic approximation, with K elastic constant, and with the dielectric anisotropy negligible. Moreover we assume e e e S B ≈ ≈ . Let us consider the director n r in a uniform configuration, as a vector parallel to x -axis and the electric field E r parallel to z -axis as  The free energy density, including the flexoelectric term, is: Then, a periodic distortion in a non-confined nematic is possible because it has a free energy density lower than that possessed by the uniform configuration. There is not a threshold for the electric field, since the existence of a threshold is a consequence of the medium confinement.
Let us imagine a nematic material confined in a cell composed by two plane walls, parallel to [x,y] plane, at a distance d. The anchoring conditions must be included in the energy balance. We can assume a surface energy density of the Rapini-Papoular form , for a surface treatment favouring a molecular alignment parallel with the x -axis. If the director field n r is uniform in the planar alignment, Let us choose, as in Ref. [18], the behaviour of the tilt angle in the form . The free energy density is given as in Equation 8.
Integrating the free energy density on the cell volume , where d is the cell thickness, L a fixed length in y-direction and Λ the director distortion wavelength along the x-direction, we obtain: The last term in (9) is the surface energy contribution. In the case of a uniform director field, we have a total energy as The behaviour of the two free energies uniform distorted F , F is given in Figure 2: we can see the existence of a threshold field * E . If the electric field has a value * E E < , the stable configuration of the director field is that with lower energy, in this case, when the director field is uniform. When * E E > , the stable configuration is the distorted one.
Comparing the two values of the total energy, that is: we can approximately find the threshold electric field as: where ( ) 0 , to have a real electric field: The threshold field has a value: and assuming the parameter values 2 4 10 , m / C a 11 10 − = we find a threshold voltage of Volt 10 ≈ .

The Chiral Nematic and the Smectic Phase
Much research has taken place in the field of liquid crystals to find ferroelectric materials, from the earlier studies on the smectic phases till the more recent banana-like materials [20][21][22][23]. The smectic phases are organised in layers. There are three main smectic phases: A, C and C*. In the smectic C (SmC) phase the director n r is tilted by a fixed angle, with respect to the layer normal ν r . The chiral smectic (SmC*) phase shows in addition an intrinsic twist of the director from layer to layer. The symmetry breaking 2 2 C C h → allows molecular electric dipoles to form a spontaneous electric polarization P r , which lies in the smectic planes. The macroscopic polarization vanishes in the SmC* phase, but an electric field parallel to the layers can distort the helicoidal structure, disfavouring SmC* and leading to a phase with a macroscopic polarization. In the Landau theory of smectic liquid crystals, the free energy is expanded in two order parameters: the projection n r of the director onto the smectic layer plane and the layer polarization P r [24]. The chiral term, responsible of the SmC* phase has in [20] the structure: ( ) This term has in fact a Lifshitz-like structure, if we consider the layer normal ν r , parallel to the z-axis: We can identify this expression as a pseudoscalar inhomogeneous Dzyaloshinskii-Moriya interaction, which does not involve an external field but a fixed direction in the space, that is the vector ν r normal to the smectic layer.
Chiral molecules can also form nematic phases called chiral nematic phases or cholesteric phases. The phase shows a nematic order, with the director rotating throughout the sample. The axis of this screw is normal to the director. The distance over which the director rotates by π 2 is the chiral pitch, generally of the order of the wavelength of visible light.
If the nematic phase is composed of chiral molecules, all of the same chirality, the material does not have symmetry planes and then the free energy has, according to Landau and Lifshitz, a pseudoscalar term: This is the pseudoscalar of the DM interaction as in Ref. 4, and introduced in Sect.2. If we consider the vector ν r as the direction of pitch, then the director n r lies in a plane perpendicular to it and then Equation 16 can be rewritten as: with the same structure that we encountered in the smectic term originating the helix.

The Saddle-Splay Elasticity at Surfaces
In nematics, a more general form of the distortion free-energy density, in the framework of the usual first-order continuum theory, is given as: where K is the bulk elastic constant in the case of elastic isotropy. The last term is the contribution of the saddle-splay elasticity. This contribution is not usually inserted in the bulk free energy, because it becomes a surface contribution when integration is performed on the cell thickness [25,26]. The saddle-splay contribution is then a Lifshitz invariant of the surface energy: (19) if ν r is the unit vector of direction perpendicular to the surface containing the nematic material. This term has the same form of Lifshitz scalar product in Equation 4. In addition to the anchoring energy, which is the anisotropic part of surface tension, there is an elastic contribution, which has been originally indicated as a part of the bulk elastic energy in the form of a divergence [27][28][29]. This contribution can be viewed as the elastic part of surface energy depending on the tangential gradient of director. The 24 K term may induce spontaneous twist deformations in hybrid nematic films with azimuthally degenerate anchoring conditions. Such deformations are manifested in the formation of periodic stripe domains observed in sufficiently thin hybrid NLC cells [25,30]. If the anchoring energy is sufficiently small, the Lifshitz term can produce a modulated-tilt state has recently shown by Lelidis and Barbero [31]. The saddle-splay contribution is necessary, when we have to evaluate the elastic contribution of thin films or membranes. In 1973, Helfrich studied the energetic cost of a generic sheet in a threedimensional space: we can determine, in each point of the sheet, the radii of curvature 1 r and 2 r , and local curvatures In this expression, C k is the bending (or curvature) modulus and G k is the saddle-splay (or Gaussian curvature) modulus. These two modules are set by interactions among membrane molecules. The spontaneous curvature is denoted by o c . As reported in [33], biological membranes are sheets that can be modelled with a continuum elastic approach. These membranes are two-dimensional fluids within which proteins diffuse and interact. Membranes can bend and curve, with deformations controlled by proteins and lipids; the converse is also true, it is the structure created by membrane curvature can guide the spatial organisation of membrane molecules. Then the membrane can display spatial patterning at length-scales far greater than the scale of individual molecules [33].

The Hybrid Cell and the Flexoelectricity
Let us start the discussion of some case studies. The first is on the role of flexoelectricity in hybrid nematic cells. A hybrid cell is a nematic cell where a sample is confined between two parallel walls with different anchoring conditions. One surface is treated to favour a planar alignment; the opposite one is favouring a homeotropic alignment. The cell is then named HAN, that is Hybrid Aligned Nematic cell. The hybrid cell we discuss has the y-axis perpendicular to cell walls (see the upper part of Figure 3).
An electric field can be applied parallel to y-axis: we have then where j r is the unit vector of y-axis. j r is the homeotropic direction too. The unit vector i r , parallel to the cell walls, gives the easy planar direction. The bulk free energy density is given, in the elastic isotropic approximation, by: where the last term is due to the dielectric anisotropy ε Δ of the nematic.   Figure 3) the energies of the homeotropic and planar cells, we see the possibility of an electric threshold field E* : under this value of the electric field, it is favoured the planar configuration, over the threshold value, it is the homeotropic configuration that has a lower energy. In a hybrid cell, the director changes from a planar configuration at one of the cell wall, to a homeotropic configuration at the other cell wall. The tilt angle is then depending on y, as a function ( ) The director field is given by: If the anchoring is strong, the tilt angle is 2 π θ = at 0 = y -homeotropic wall, and 0 = θ at d y = -planar wall. In the one elastic constant approximation, we have the bulk free energy density in the form: and the surface energy density To represent a hybrid configuration in a very rough approach, let us simply choose a linear function of the tilt angle with y. Then: and the total bulk energy is: Including the surface energy, the total energy is:   What is shown in Figure 4(a) is surely possible, because we can adjust the anchoring parameters. We observe then two threshold fields: when the field is lower than E' , the nematic is planar, if the field is comprised between E' and ' E' , the cell is hybrid. Over ' E' , the cell is homeotropic. As previously discussed, the electric field can be coupled with a polarization arising from an elastic deformation in the flexoelectric effect. In planar and homeotropic configurations, because there are not deformations of director, the flexoelectric effect is absent, but in the hybrid cell the deformation gives a flexoelectric polarization If θ is given by Equation (24), after integrating on the cell volume, we have the contribution of flexoelectricity to the total free energy as: In principle, the coefficient ( ) are changed from the contribution of the flexoelectricity. They could be lowered or raised by the induced polarization (see Figure 4(b)). The thresholds change according to the shape of the molecules. Comparing the thresholds we can estimate the values of the coefficients. The two electric field contributions in the HAN cell are: If they were of the same order of magnitude, we could obtain: In the case of a cell with a thickness of m μ 10 , a field of m V μ 10 , and an electric anisotropy as we obtain: in agreement with Ref. 35 and with other experimental values [36][37][38][39][40]. Recently a giant flexoelectricity has been found with bent-core nematics: a peak of m nC 35 was measured in these materials then more than 3 orders of magnitude larger than in calamitics [41]. In the next section we will study the alignment transitions in the nematic cells; such a problem was studied also in Ref. 42.

The Phase Diagrams of the Hybrid Cell
Let us consider the hybrid cell as in the previous section. We use the same notation here but we solve the Euler-Lagrange equation with the proper boundary conditions, by means of an iterative procedure previously used in Ref. 43, to investigate the ion densities in corona plasma. The Euler-Lagrange equation is: and: that is: (36) and the surface energy density: Equation 36 can be rigorously solved with elliptic functions: in Ref. 44, the existence of a critical thickness for a hybrid aligned nematic cell was predicted in the framework of a rigours solution. Nevertheless, numerical solutions and an approximate analytic theory have been already used, when studying the existence of flexoelectric instabilities in the case of asymmetric boundary conditions [45]. Here, we use an approximate solution to the non-linear problem of the form: and then Equation 36 can be written as two equations:  The boundary conditions are: (43) in which we used the dimensionless parameters . From the first equation in the boundary conditions (43): Once we chose the value of o α , from Equation (44), we have the value of 1 α , and then, after To determine the value of o α we could use the other boundary condition, the second in (43); but, in this case, we are facing a strongly oscillating function. It is better to determine the value of parameter o α , minimizing the total free energy. Adding the flexoelectricity, the term to include in the free energy density is: and, after integration on the cell thickness, we have a further contribution to the surface energy density of the form: This term can be easily inserted in the numerical calculation, to minimize the total free energy. Let us introduce the following dimensionless variables and parameter: to illustrate the results of calculations. In Figure 5 we can see the phase diagrams of the HAN cell, for a fixed choice of the surface . We can change the value of parameter 1 b and find the value of the threshold field (the electric field is dimensionless represented by ξ ) ). There are three regions in the diagrams where planar, homeotropic and hybrid alignments are allowed according to the values of the electric field. The phase diagram is depending on the values of flexoelectric parameter Π (see diagrams (a),(b) and (c) in Figure 5). The last diagram (d) shows the behaviour of a cell when we change the flexoelectric parameter Π . Note that the hybrid configuration disappears when flexoelectric parameter is higher than value 1.3.
In Figure 6, we see the behaviour of θ cos as a function of the dimensionless variable d y y = ) in the case of positive and negative flexoelectric coefficients, for different values of the electric field. Note that, as the field increases, the role of surface is suppressed and the angle at the planar surface increases. As the electric field is higher than the threshold value, the cell becomes homeotropic and 0 = θ cos . Note that, as the field increases, the role of surface is suppressed and the angle at the planar surface increases. As the electric field is higher that the threshold value, the cell becomes homeotropic and then 0 = θ cos .
To conclude this section on HAN cells, let us remember that we have another Lifshitz invariant, that gives the saddle-splay contribution to the surface free energy density, in the form: . In the hybrid cell alignment, where only the tilt angle is displayed by the elastic distortion, this contribution is zero. We will see in the last section of the paper, how this term produces a periodic distortion and how the PHAN -the Periodic HAN -texture appears. The fact that the saddle-splay contribution is zero in the HAN configuration, is in agreement with the conclusion that in the same configuration the flexoelectric contribution

Nematics in Cylindrical Geometry
Let us consider a cylinder with radius R . In this cylindrical cell we imagine to insert a nematic. We use the frame as in Figure 7 and solve the Euler-Lagrange equation in cylindrical coordinates.
The bulk density energy is given by: and then the Euler-Lagrange equation turns out to be: The surface energy density is: . For an anchoring, which favours an homeotropic alignment of the nematic perpendicular at the wall of the cylinder, we use: If we want to avoid the presence of a defect at the axis of cylinder ( − z axis), the director must escape in the − z direction. The solution, if the applied electric field is zero, is given by an inverse tangent: , for a strong anchoring at the cylinder wall. This is a well-know solution due to Belavin and Polyakov [46]. To solve the equation in the case of electric field different from zero, we choose a solution as: The use of an approximate solution could be questionable. Nevertheless, linearization and approximation of Belavin-Polyakov equation are reported in the literature [45,47]. Using (54), we have two equations to solve: The second equation can be solved with iterations. At the fourth step of iteration the solution is within 0.1%. In the following way, we have: Actually, we arrive at the following solutions: and then at final solution ( ) ( ) ( ) To determine the value of parameter β we choose the solution minimizing a reduced total free energy: where L is an arbitrary length of the cylindrical cell. Let us then consider the contribution of flexoelectricity to Euler-Lagrange equations: We use again ( ) and K WR b = as parameters. The equation to solve is: instead of Equation (56). Figures 8 and 9 show the results of calculations for different values of anchoring and flexoelectric coefficients. In the Figure 8 we can see the angle θ as a function of the reduced radial distance R r , for two values of the flexoelectric coefficient, 0 = Π and 1. The figure shows the behaviour in the case of different values of anchoring parameter b and of dimensionless electric field parameter R ξ = ξ . As the electric field is higher that a threshold value, angle θ goes to zero and the director field is parallel to cylinder axis in all the cell. The following Figure 9 shows θ cos as a function of reduced radial distance R r , for different values of Π and ξ . In this case, the value of the anchoring strength is fixed. Note that a negative value of the flexoelectric parameter is strongly favouring the alignment of the director parallel to cylinder axis, and then we find a low value of the threshold electric field. If flexoelectric parameter Π is positive and large, the distorted configuration is favoured, and the threshold field required for suppressing this configuration is increased. Moreover, if the flexoelectric parameter is large, as in the lower image in Figure 9, angle θ starts to oscillate as the field increases. We must have a huge electric field to suppress the oscillating distortion and have 0 = θ cos , with all the nematic aligned parallel to the field, in a uniform configuration.
In Figure 10, the phase diagrams are shown, when anchoring parameter b is fixed and equal to 6. We see three regions, denoted by: U for uniform alignment of director parallel to z-axis, D if the director has a deformed configuration, and O when the director is oscillating and cosine becomes negative too. Angle θ turns more than 2 π on the distance R. As told before, giant flexoelectric coefficients are possible and then the oscillation could be experimentally tested in cylindrical cells. A last note on the flexoelectric term. The flexoelectric vector is a sum of two contributions: In the previous assumptions,  As the electric field is higher that a threshold value, angle θ goes to zero, that is the director field is parallel to the cylinder axis. in all the figures. Note that a negative value of the flexoelectric parameter is strongly favouring the alignment of director parallel to cylinder axis, and then the threshold electric field is very low. If the flexoelectric parameter is positive and large, the distorted configuration is favoured, and the threshold field, needed for suppressing this configuration, is increased. As shown in the lower part of the figure, when Π is very large, θ cos is oscillating as the electric field increases. A very large field is required to suppress the distortion and have 0 = θ . Figure 10. Phase diagram of the cylindrical confinement, when the anchoring parameter b is fixed and equal to 6. The three regions are denoted by U for the uniform alignment of director parallel to the cylinder axis, D when the director has a deformed configuration, and O if director is oscillating and cosine becomes negative too.

The Saddle-Splay Contribution and the PHAN Cell
Sometimes, it is possible to note a periodicity in the HAN cells observed by the polarised light microscope. Because of this periodic configuration, the cell is in the PHAN configuration, that is a nematic cell with a period hybrid alignment. Two angles describe the PHAN configuration: θ and φ. The last angle is formed by the projection of the director in the plane of the cell with the x-axis.
The free energy density is that by Nehring and Saupe, and given by Equation (19). The frame of reference is [ ] xyz , with [ ] xy the cell plane and − z axis perpendicular to the cell plane. The homeotropic wall is at z 0 = 0, where z is the axis perpendicular to the cell plane. The planar wall is at z 1 = d, where d is the thickness of the cell. The easy-axis of the planar alignment is chosen coincident with the x-axis. The director n r is described as: The Euler-Lagrange equations are non-linear. They were solved in Ref. 25, with a numerical approach to determine the threshold thickness of the cell between the planar and the PHAN. Here, we want to grasp the role of the saddle-splay contribution, with just simple calculations. Let us then consider the tilt angle θ depending on z, and the φ angle depending on x, in the following way: The tilt is zero if 0 = z , and it is 2 π at z = d. With Λ we denote the wavelength along − x axis.
The free energy density is: where P W W = for planar anchoring with 0 = θ , and H W W = at the homeotropic anchoring 2 π θ = .
After integrating on surfaces of the cell: and then the total free energy is: Neglecting the last term, we find a threshold value for the cell thickness:

Conclusions
This paper is divided in two parts. In the first we have discussed the analogies among Lifshitz invariants in magnetic materials and liquid crystals. We saw that the structure of these invariants is the same, and that they are producing periodic instabilities in both cases. In the Lifshitz invariant, the interaction is between an external action and the order parameter, in a form that contains the gradient of order parameter. The external action can be an electric field applied to the bulk, and in this case the relevant effect is the flexoelectricity, or the confinement due to free surfaces or cell walls. The Lifshitz invariant related to surfaces gives the saddle-splay contribution to surface energy.
In the second part of the paper we discussed in depth the role of flexoelectricity in the case of confined nematics. We performed detailed calculations in the case of planar and cylindrical geometry. Phase diagrams are also shown, to see the alignment phase transitions due to electric field and the role of flexoelectric parameter.