Molecular dynamics and Monte Carlo simulations are used to determine the stable states of the kaolinite-methanol intercalation system. Several evaluation scenarios are presented to describe the stability of the system using single mu SolventNClaypT ensemble simulations and a series of NpT ensemble simulations coupled with thermodynamic integration. (C) 2015 Elsevier B.V. All rights reserved.