A spin dynamics algorithm, combining checkerboard updating and a rotation algorithm based on the local second-rank ordering field, is developed for the Lebwohl-Lasher model of liquid crystals. The method is shown to conserve energy well and to generate simulation averages that are consistent with those obtained by Monte Carlo simulation. However, care must be taken to avoid the undesirable effects of director rotation, and a method for doing this is proposed.

• Received 21 May 2005

DOI:https://doi.org/10.1103/PhysRevE.72.036703