Microscopic shear flow simulations of a biaxial smectic A liquid crystal based on the soft ellipsoid string-fluid

We have studied the behaviour of a biaxial smectic A liquid crystal based on the soft ellipsoid string-fluid in shear flow by molecular dynamics simulation using the SLLOD equation of motion. This is facilitated by the fact that the biaxial symmetry allows linear relations between the pressure and the velocity gradient. This means that linear irreversible thermodynamics can be applied independently of the simulations to obtain the torques determining the orientations of the system and that the predictions of this theory can be cross-checked by the simulations. It turns out that there is a torque turning the smectic layers to the orientation parallel to the vorticity Wane if the simulation is started in another orientation. In the orientation parallel to the vorticity Wane where the director formed by the long axes of the molecules, n(w), is perpendicular to the vorticity Wane there is another torque keeping the director formed by the normals of the broadsides of the molecules, n(u), parallel to this Wane at a constant alignment angle, psi relative to the streamlines independently of the strain rate. Moreover, this alignment angle seems to be the one where the irreversible energy dissipation rate (w) over dot, is minimal. This is in agreement with a recently proven theorem according to which (w) over dot is minimal in the linear regime of a nonequilibrium steady state. Finally, we studied the orientation of n(u) when the smectic layers are parallel to the shear Wane. In a simulation this orientation is stabilised by the periodic boundary conditions. Then we found that there was a nonlinear torque turning n(u) to the orientation perpendicular to the streamlines thus minimising the value of (w) over dot even though this value is larger than the value of (w) over dot in the orientation parallel to the vorticity Wane. This means that (w) over dot is minimized given the external boundary conditions.