Monte Carlo simulations for the free energies of C60 and C70 fullerene crystals by acceptance ratio method and expanded ensemble method

J Chem Phys. 2014 Feb 28;140(8):084110. doi: 10.1063/1.4866451.

Abstract

Accurate values of the free energies of C60 and C70 fullerene crystals are obtained using expanded ensemble method and acceptance ratio method combined with the Einstein-molecule approach. Both simulation methods, when tested for Lennard-Jones crystals, give accurate results of the free energy differing from each other in the fifth significant digit. The solid-solid phase transition temperature of C60 crystal is determined from free energy profiles, and found to be 260 K, which is in good agreement with experiment. For C70 crystal, using the potential model of Sprik et al. [Phys. Rev. Lett. 69, 1660 (1992)], low-temperature solid-solid phase transition temperature is found to be 160 K determined from the free energy profiles. Whereas this is somewhat lower than the experimental value, it is in agreement with conventional molecular simulations, which validates the methodological consistency of the present simulation method. From the calculations of the free energies of C60 and C70 crystals, we note the significance of symmetry number for crystal phase needed to properly account for the indistinguishability of orientationally disordered states.