Phase diagrams of alkali halides using two interaction models: a molecular dynamics and free energy study

Phase diagrams for potassium and sodium chlorides are determined by molecular dynamics and free energy calculations. Two rigid-ion interaction models, namely, the Born-Mayer-Huggins (BMH) and Michielsen-Woerlee-Graaf (MWG) effective pair potentials, have been used. The critical and triple point properties are discussed and compared with available experimental and simulation data. The MWG model reproduces the experimental liquid-gas equilibria better than the BMH model, being the accordance very good in the lowest temperature region of the coexistent liquids, particularly for NaCl. However, both models underestimate the critical temperatures of KCl and NaCl. Relatively to the solid-gas equilibria, the models do not reproduce well the experimental data. As for the solid-liquid coexistences either the BMH or the MWG models appear unrealistic.