3.3. Liquid crystal order¶
Various standard initialisations for the tensor order parameter \(Q_{\alpha\beta}\) are available and are discussed below.
3.3.1. Nematic initialisation¶
Remember that to obtain a nematic, the cholesteric pitch wavevector \(q_0\) should be set to zero.
The following initial conditions are specifically aimed at nematics (although other initialisations discussed for cholesterics below may also be relevant).
lc_q_initialisation nematic # uniform nematic...
lc_init_nematic 1.0_0.0_0.0 # director n_alpha
Here, based on the input director \(n_\alpha\), the tensor order parameter will be initialised uniformly in uni-axial fashion:
The value of the amplitude \(A_0\) is taken from the corresponding current free energy parameter. If the input director is not a unit vector, it is normalised to be so before the initialisation.
lc_q_initialisation random_xy # random in (x,) plane
Will select a random unit director in the (x,y) plane and initialise the \(Q_{\alpha\beta}\) tensor as above.
3.3.2. Cholesteric initialisation¶
The following simple cholesteric initialisations are available:
lc_q_initialisation cholesteric_x
lc_q_initialisation cholesteric_y
lc_q_initialisation cholesteric_z
These give, respectively, a helical axis along the \(x\)-direction, \(y\)-direction or \(z\)-direction. For example, if the \(x\)-direction is selected, then a director \([0,\cos(q_0 x), \sin(q_0 x)]\) is computed and the tensor order parameter initialised from the unixial approximation
with \(A_0\) taken form the corresponding free energy parameter.
For the \(y\)-direction the uniaxial director is \([\cos(q_0 y),0,-\sin(q_0 y)]\), and for the \(z\)-direction the director is \([\cos(q_0 z), \sin(q_0 z),0]\).
3.3.3. Blue phase initialisations¶
The reader is referred to Wright and Mermin [WrightAndMermin1989] for the theoretical background on blue phases.
3.3.3.1. Blue phase I (“o8m”)¶
An approximation to the order parameter for BPI (referred to as \(O^{8-}\) by Wright and Mermin) is available in the high-chirality limit, and can be used as the basis for initialisation.
A suitable initialisation might be:
lc_q_initialisation o8m # Initialisation is BPI
lc_q_init_amplitude -0.2 # amplitude A
If we write \(C_x = \cos(\sqrt{2} q_0 x)\), \(S_x = \sin(\sqrt{2} q_0 x)\) etc, then the order parameter components are:
and so on. Note that a negative amplitude \(A\) is associated with oblate anisoptropy (+ve is prolate).
The wavevector is typically selected so that the system length is \(L = \sqrt{2} p\), that is, a unit cell.
3.3.3.2. Blue phase II (“o2”)¶
The approximation in the high chirality limit referred to by Wright and Mermin as \(O^2\) is appropriate for blue phase II, and is introduced via
lc_q_initialisation o2
lc_q_init_amplitude 0.3
Again using \(C_x\) and \(S_x\) the approximation is:
and so on. Note the amplitude here is positive.
3.3.3.3. Rotated blue phase initialisations¶
An additional rotation can be applied to the following blue phase initialisations to adjust their orientation:
lc_q_initialisation o8m # Available
lc_q_initialisation o2 # Available
lc_q_init_euler_angles 45_45_0 # e.g., Euler angles (a,b,c).
The rotation applied for Euler angles (\(\alpha, \beta, \gamma\)) is a standard Euler rotation
representing a sequence of rotations: first, around the \(z\)-axis by angle \(\alpha\) to obtain rotated frame \((x', y',z)\); second around the new \(x'\)-axis to obtain frame \((x', y'', z')\); and third around the \(z'\)-axis to obtain the final result. The rotation is around the centre of the system.
The angles entered should be in degrees.
Note that an arbitrary rotation of the initial conditions based on the high chirality limit may introduce frustration at the periodic boundaries if the disclination lines do not match. No such issue occurs if boundaries are not periodic.
3.3.4. Miscellaneous liquid crystal initialisations¶
D.C. Wright and N.D. Mermin, Crystalline liquids: the blue phases, Rev. Mod. Phys., 61 385-432 (1989).