Molecular dynamics study of the vaporization of an ionic drop

The melting of a microcrystal in vacuum and subsequent vaporization of a drop of NaCl were studied through molecular dynamics simulations with the Born-Mayer-Huggins-Tosi-Fumi rigid-ion effective potential. The vaporization was studied for a single isochor at increasing temperatures until the drop completely vaporized, and gaseous NaCl formed. Examination of the vapor composition shows that the vapor of the ionic drop and gaseous NaCl are composed of neutral species, the most abundant of which, ranging from simple NaCl monomers (ion pairs) to nonlinear polymers, (Na(n)Cl(n))(n=2-4). The enthalpies of sublimation, vaporization, and dissociation of the different vapor species are found to be in reasonable agreement with available experimental data. The decrease of the enthalpy of vaporization of the vapor species, with the radius of the drop decrease, accounts for a larger fraction of trimers and tetramers than that inferred from experiments. Further, the rhombic dimer is significantly more abundant than its linear isomer although the latter increases with the temperature. The present results suggest that both trimers and linear dimers may be important to explain the vapor pressure of molten NaCl at temperatures above 1500 K.