We derive exact expressions for the forces and torques between biaxial molecules interacting via the RE-squared potential, a recent variant of the Gay-Berne potential. Moreover, efficient routines have been provided for rigid body MD simulations, resulting in 1.6 times speedup compared to the two-point finite difference approach. It has also been shown that the time cost of a MD simulation will be almost equal to a similar MC simulation, making use of the provided routines.
@article{arxiv.cond-mat/0605193,
title = {Analytical First Derivatives of the RE-squared Interaction Potential},
author = {M. Babadi and M. R. Ejtehadi and R. Everaers},
journal= {arXiv preprint arXiv:cond-mat/0605193},
year = {2009}
}