TY - JOUR

T1 - Calculation of the melting point of NaCl by molecular simulation

AU - Anwar, J

AU - Frenkel, D

AU - Noro, MG

PY - 2003/1/8

Y1 - 2003/1/8

N2 - We report a numerical calculation of the melting point of NaCl. The solid-liquid transition was located by determining the point where the chemical potentials of the solid and liquid phases intersect. To compute these chemical potentials, we made use of free energy calculations. For the solid phase the free energy was determined by thermodynamic integration from the Einstein crystal. For the liquid phase two distinct approaches were employed: one based on particle insertion and growth using the Kirkwood coupling parameter, and the other involving thermodynamic integration of the NaCl liquid to a Lennard-Jones fluid. The latter approach was found to be significantly more accurate. The coexistence point at 1074 K was characterized by a pressure of -30+/-40 MPa and a chemical potential of -97.9+/-0.2k(beta)T. This result is remarkably good as the error bounds on the pressure enclose the expected coexistence pressure of about 0.1 MPa (ambient). Using the Clausius-Clapyron relation, we estimate that dP/dTapproximate to3 MPa/K. This yields a melting point of 1064+/-14 K at ambient pressure, which encompasses the quoted range for the experimental melting point (1072.45-1074.4 K). The good agreement with the experimental melting-point data provides additional evidence that the Tosi-Fumi model for NaCl is quite accurate. Our study illustrates that the melting point of an ionic system can be calculated accurately by employing a judicious combination of free energy techniques. The techniques used in this work can be directly extended to more complex, charged systems. (C) 2003 American Institute of Physics.

AB - We report a numerical calculation of the melting point of NaCl. The solid-liquid transition was located by determining the point where the chemical potentials of the solid and liquid phases intersect. To compute these chemical potentials, we made use of free energy calculations. For the solid phase the free energy was determined by thermodynamic integration from the Einstein crystal. For the liquid phase two distinct approaches were employed: one based on particle insertion and growth using the Kirkwood coupling parameter, and the other involving thermodynamic integration of the NaCl liquid to a Lennard-Jones fluid. The latter approach was found to be significantly more accurate. The coexistence point at 1074 K was characterized by a pressure of -30+/-40 MPa and a chemical potential of -97.9+/-0.2k(beta)T. This result is remarkably good as the error bounds on the pressure enclose the expected coexistence pressure of about 0.1 MPa (ambient). Using the Clausius-Clapyron relation, we estimate that dP/dTapproximate to3 MPa/K. This yields a melting point of 1064+/-14 K at ambient pressure, which encompasses the quoted range for the experimental melting point (1072.45-1074.4 K). The good agreement with the experimental melting-point data provides additional evidence that the Tosi-Fumi model for NaCl is quite accurate. Our study illustrates that the melting point of an ionic system can be calculated accurately by employing a judicious combination of free energy techniques. The techniques used in this work can be directly extended to more complex, charged systems. (C) 2003 American Institute of Physics.

KW - MODEL

KW - HIGH-TEMPERATURE

KW - DIAGRAM

KW - SYSTEM

KW - FREE-ENERGY

KW - CRYSTALS

KW - NITROGEN

KW - PHASE COEXISTENCE

KW - TRANSITIONS

KW - INTEGRATION

U2 - 10.1063/1.1522375

DO - 10.1063/1.1522375

M3 - Journal article

VL - 118

SP - 728

EP - 735

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 2

ER -