The surface chemistry of divalent metal carbonate minerals; a critical assessment of surface charge and potential data using the charge distribution multi-site ion complexation model

Abstract

The Charge Distribution MUltiSite Ion Complexation or CD–MUSIC modeling approach is used to describe the chemical structure of carbonate mineralaqueous solution interfaces. The new model extends existing surface complexation models of carbonate minerals, by including atomic scale information on the surface lattice and the adsorbed water layer. In principle, the model can account for variable proportions of face, edge and kink sites exposed at the mineral surface, and for the formation of inner- and outer-sphere surface complexes. The model is used to simulate the development of surface charges and surface potentials on divalent carbonate minerals as a function of the aqueous solution composition. A comparison of experimental data and model output indicates that the large variability in the observed pH trends of the surface potential for calcite may in part reflect variable degrees of thermodynamic disequilibrium between mineral, solution and, when present, gas phase during the experiments. Sample preparation and non-stoichiometric surfaces may introduce further artifacts that complicate the interpretation of electrokinetic and surface titration measurements carried out with carbonate mineral suspensions. The experimental artifacts, together with the high sensitivity of the model toward parameters describing hydrogen bridging and bond lengths at the mineralwater interface, currently limit the predictive application of the proposed CD–MUSIC model. The results of this study emphasize the need for internally consistent experimental data sets obtained with well-characterized mineral surfaces and in situ aqueous solution compositions (that is, determined during the charge or potential measurements), as well as for further molecular dynamic simulations of the carbonate mineral-water interface to better constrain the bond lengths and the number plus valence contribution of hydrogen bridges associated with different structural surface sites.