Multiple Pathways in NaCl Homogeneous Crystal Nucleation

Abstract: NaCl crystal nucleation from metastable solutions has long been considered to occur according to a single-step mechanism. Recent experimental observations suggest that NaCl crystals emerge from disordered intermediates which is seemingly at odds with this established view. Here, we performed biased and unbiased molecular dynamics simulations to analyse and characterise the pathways to crystalline phases from solutions far into the metastable region. We find that large liquid-like NaCl clusters emerge as the solution concentration is increased and two-step nucleation pathways are more dominant than one-step pathways to phase separation. Analyses of cluster size populations and the ion pair association constant show that these clusters are transient, unlike the thermodynamically stable prenulceation cluster solute species that were suggested in other mineralising systems. A Markov State Model was developed to analyse the mechanisms and timescales for nucleation from unbiased molecular dynamics trajectories in a reaction coordinate space characterising the dense regions in clusters and crystalline order. This allowed calculation of the committor probabilities for the system to relax to the solution or crystal states, and to estimate the rate of nucleation, which shows excellent agreement with literature values. From a fundamental nucleation perspective, our work highlights the need to extend the attribute 'critical' to an ensemble of clusters which can display a broad range of structures and include sizeable disordered domains depending upon the reaction conditions.Moreover, our recent simulation studies demonstrated that carbon surfaces catalyse the formation of liquid-like NaCl networks which, combined with the observations here, suggests that alternative pathways beyond the single-step mechanism can be exploited to control the crystallisation of NaCl.