Calculation of excess free energy of molecular solids comprised of flexible molecules using Einstein molecule method

Abstract

Knowledge of free energies of crystalline solid phases is pivotal in determination of phase equilibria; however, their computations using molecular simulations are quite challenging. In this paper, we apply the Einstein molecule method to compute the excess free energies of crystalline molecular solids comprising of molecules described by flexible models. This work can be viewed as an extension of the method described in Aragones et al. [J Chem Phys. 2013;139:a034104]. In this article, we present our calculations of the free energies of ice-Ih described by q-TIP4P/F model and Orcinol Form I described by the OPLS-AA forcefield. The free energies computed for ice-Ih are in agreement with previous reported results, and those of Orcinol are shown to be thermodynamically consistent. We demonstrate that the efficiency of our calculations compares favourably with existing methods. The method described here can be easily implemented in popular publicly available molecular dynamics packages without any modification to their existing source codes.