Molecular dynamics simulations of a Gay–Berne nematic liquid crystal: Elastic properties from direct correlation functions

We report molecular dynamics simulations of a Gay–Berne nematic liquid crystal at constant temperature and density/pressure using the generalization of an algorithm recently proposed by Toxvaerd [Phys. Rev. E 47, 343 (1993)]. On the basis of these simulations the absolute values of the Oseen–Zocher–Frank elastic constants K11, K22, and K33 as well as the surface constants K13 and K24 have been calculated ab initio within the framework of the direct correlation function approach of Lipkin et al. [J. Chem. Phys. 82, 472 (1985)]. The angular coefficients of the direct pair correlation function, which enter the final equations, have been determined from the computer simulation data for the pair correlation function of the nematic by combining the Ornstein–Zernike relation and the Wiener–Hopf factorization scheme. The unoriented nematic approximation has been assumed when constructing the reference state of Lipkin et al. By an extensive study of the model over a wide range of temperatures, densities and pressu...