Revisiting the birth of NaCl crystals using molecular dynamics simulation

Abstract

Molecular dynamics simulations were used to investigate the initial stage of phase separation mechanisms for an oversaturated electrolytic solution. We developed a low computational cost methodology to determine the simulation frames where the first ionic clusters are formed. By discretizing the simulation box, we obtain a density profile in the moments preceding and succeeding the nuclei's formation. The growth of the clusters identified with our methodology was analyzed until the end of the simulation. Calculation of the Steinhardt parameter showed symmetry of the solid, giving indications that the classical nucleation theory explains the mechanism of the solid formation. The methodology developed was useful for identifying phase separation mechanisms in the nucleation process. At lower concentrations, there was no formation of stable clusters. At intermediate concentrations, the analyses indicate a transition of phases in one stage, from a oversaturate electrolytic solution to a crystalline solid. At high concentration, a transition of phases in two stages, initially, is the formation of a dense liquid, and only after that, crystalline solid formed inside the dense liquid. The change in phase separation mechanism due to increasing oversaturation underscores the importance of precise determination of the driving force for phase separation and concentration limits for each mechanism.

Keywords: Classic nucleation theory; Molecular dynamics; Phase separation mechanisms; Solidification; Spinodal decomposition.