How Ion Pair Formation Drives Adsorption In The Electrical Double Layer: Molecular Dynamics Of Charged Silica-Water Interfaces In The Presence Of Divalent Alkaline Earth Ions

The study by numerical methods of ionic distributions in charged solid-liquid interfaces allows the interpretation of many concepts and phenomena such as the zeta potential or ion adsorption. Molecular dynamics (MD) can reveal detailed information about electric double layer (EDL), especially the Stern layer, by including electrostatic, van der Waals, and molecular forces. Here, we aim at analyzing the chemical species and the correlations between the ions and the surface including three-body effects, one particle belonging to the surface and two to the solution. Correlations are specified on the basis of interionic distance screening. An extensive description is provided from simulations of three alkaline earth metal chlorides (Mg2+, Ca2+, and Ba2+) aqueous solutions at 0.6 mol.L-1 in the centers of negatively charged silica nanochannels. The resulting McMillan-Mayer potentials of mean force (PMF) exhibit a decreasing affinity of deprotonated silanol along with the series Mg2+ > Ca2+ > Ba2+, while the formation of bulk M2+-Cl- pairs is in reverse order. A similar trend is obtained for the association constants and the residence times. Over 40% of surface-bound Ba2+ ions are correlated with surface-bound Cl-, while the other two cations do not show such trend. When the surface-bound and surface-correlated ions are taken apart, the remaining free ion distributions fit the Poisson-Boltzmann equation well, i.e., the Gouy-Chapman model. This work demonstrates the necessity to account for three-body associations on oxide surface at least for divalent ions.