Solubility of NaCl in water and its melting point by molecular dynamics in the slab geometry and a new BK3-compatible force field

Saturated concentration of rock salt in water is determined by a simulation of brine in contact with a crystal in the slab geometry. The NaCl crystals are rotated to expose facets with higher Miller indices than [001] to brine. The rock salt melting point is obtained by both the standard and adiabatic simulations in the slab geometry with attention paid to finite size effects as well as to a possible influence of facets with higher Miller indices and applied stress. Two force fields are used, the Lennard-Jones-based model by Young and Cheatham with SPC/E water and the Kiss and Baranyai polarizable model with BK3 water. The latter model is refitted to thermomechanical properties of crystal NaCl leading to better values of solubility and the melting point.