It is critical to improve the response speed of a liquid crystal wavefront corrector in order to increase the bandwidth of a liquid crystal adaptive optics system. The design of liquid crystal molecules with small rotational viscosity becomes a basic method of increasing the response speed of a liquid crystal wavefront corrector. Various phases of liquid crystal from molecular dynamics simulation are given in this paper, and the detailed computational methods of order parameter and rotational viscosity are also presented. Rotational viscosities of liquid crystals are compared based on the molecular dynamics of mixtures. The data fluctuation is reduced effectively through several simulations and the multiple analysis of original data. A detailed process of molecular dynamics of mixtures is given in this paper and the result is greatly satisfactory. We believe that one can perform a better molecular design using this process and obtain a better understanding of molecular interactions of LCs.