Flexoelectricity, incommensurate phases and the Lifshitz point

The solutions for the minimizers of the energy density f (q, p)  =  Aq² + Bq⁴ + p² + gA,B + β(q'p - p'q)+ |q'|² +κ|p'|²] describe the flexoelectric effect with a flexoelectric coupling coefficient β. The order parameters q and p can be visualized as strain and polarisation, respectively. The parameter κ denotes the ratio of intrinsic length scales for q and p. We show that the structural ground-states include 3 phases, namely the paraelastic state q  =  p  =  0, the ferroelastic state where polarization exists inside and near twin boundaries, and the incommensurate (modulated) phases with a very rich array of structural modulations ranging from nearly pure sine waves to kink arrays and ripple states. The phases coincide in the multicritical Lifshitz point. Linear flexoelectricity p~q' is encountered only approximately inside the ferroelastic phase and near the phase boundary between the paraelastic phase and the incommensurate phase. The relationship between the polarisation and the strain gradient is highly non-linear in all other states. The spatial profiles and energy distributions are discussed in detail.