Skip to main content
Article thumbnail
Location of Repository

Calculation of the melting point of NaCl by molecular simulation.

By Jamshed Anwar, D. Frenkel and M.G. Noro

Abstract

NoWe 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ß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/dT~3 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

Topics: Sodium compounds, Melting point, Chemical potential, Digital simulation, Molecular dynamics method, Free energy, Solid-liquid transformations
Year: 2009
DOI identifier: 10.1063/1.1522375
OAI identifier: oai:bradscholars.brad.ac.uk:10454/3975
Provided by: Bradford Scholars
Download PDF:
Sorry, we are unable to provide the full text but you may find it at the following location(s):
  • http://hdl.handle.net/10454/39... (external link)
  • http://dx.doi.org/10.1063/1.15... (external link)
  • Suggested articles


    To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.