Communication: A simple method for simulation of freezing transitions

Despite recent advances, precise simulation of freezing transitions continues to be a challenging task. In this work, a simulation method for fluid-solid transitions is developed. The method is based on a modification of the constrained cell model which was proposed by Hoover and Ree [J. Chem. Phys. 47, 4873 (1967)]. In the constrained cell model, each particle is confined in a single Wigner-Seitz cell. Hoover and Ree pointed out that the fluid and solid phases can be linked together by adding an external field of variable strength. High values of the external field favor single occupancy configurations and thus stabilize the solid phase. In the present work, the modified cell model is simulated in the constant-pressure ensemble using tempering and histogram reweighting techniques. Simulation results on a system of hard spheres indicate that as the strength of the external field is reduced, the transition from solid to fluid is continuous at low and intermediate pressures and discontinuous at high pressures. Fluid-solid coexistence for the hard-sphere model is established by analyzing the phase transition of the modified model in the limit in which the external field vanishes. The coexistence pressure and densities are in excellent agreement with current state-of-the-art techniques.