Direct calculation of solid-liquid equilibria from density-of-states Monte Carlo simulations

A density-of-states Monte Carlo method is proposed for simulations of solid-liquid phase equilibria. A modified Wang-Landau density-of-states sampling approach is used to perform a random walk in regions of potential energy and volume relevant to solid-liquid equilibrium. The method provides a direct estimate of the relative density of states [Omega(U,V)] and thus the relative free energy within these regions, which is subsequently used to determine portions of the melting curve over wide ranges of pressure and temperature. The validity and usefulness of the method are demonstrated by performing crystallization simulations for the Lennard-Jones fluid and for NaCl.