Surface conductivity in the electrical double layer (EDL) is known to be affected by proton hopping and diffusion at solid-liquid interfaces. Yet, the role of surface protolysis and its kinetics on the thermodynamic and transport properties of the EDL are usually ignored as physical models consider static surfaces. Here, using a novel molecular dynamics method mimicking surface protolysis, we unveil the impact of such chemical events on the system’s response. Protolysis is found to strongly affect the EDL and electrokinetic aspects with major changes in Formula Presented potential and electro-osmotic flow.

Molecular Dynamics (MD) simulations are uniquely suitable for providing molecular-level insights into the Electric Double Layer (EDL) that forms when a charged surface is in contact with an aqueous solution. However, simulations are only as accurate in predicting EDL properties as permitted by the atomic interaction models. Experimental ζ-potential values and surface charges could provide a potentially suitable reference to validate and tune the interaction models, if not for the fact that they themselves are a product of imperfect models used to interpret the raw measurement data. Here, we present an approach to tune an interaction model by comparing Electro-Osmotic Flow (EOF) MD simulations against experimental Streaming Current (SC) measurements while minimizing potential modeling errors arising from both approaches. The point that is least susceptible to interpretation and modeling errors is argued to be at the concentration for which zero flow velocity is observed in EOF simulations and a net zero electric current is measured in SC experiments. At this concentration, the ζ-potential is also zero. We were able to match the experimental concentration at which ζ = 0 in MD simulations for a CaCl2 solution at pH 7.5 in contact with fused silica by tuning the ion-surface Lennard-Jones cross interactions. These interactions were found to greatly affect the ion distribution within the EDL and particularly the formation of inner-sphere surface-complexes, which, in turn, affects the electrokinetic flow. With the ion distribution determined explicitly, a series of properties can be calculated unambiguously, such as the capacitance needed for surface complexation models.

Preferential ion adsorption in mixed electrolytes plays a crucial role in many practical applications, such as ion sensing and separation and in colloid science. Using all-atom molecular dynamics simulations of aqueous NaCl, CaCl₂, and NaCl-CaCl₂ solutions confined by charged amorphous silica, we show that Na⁺ ions can adsorb preferentially over Ca²⁺ ions, depending on the surface structure. We propose that this occurs when the local surface structure sterically hinders the first hydration shell of the Ca²⁺ ion. Introducing a protrusion metric as a function of protrusion of deprotonated silanols, ion-specificity is successfully predicted on isolated, vicinal, and geminal silanols alike, provided that no other deprotonated silanols are found nearby. Furthermore, we introduce a new strategy to analyze the results as a function of distance from the surface. This approach effectively removes surface roughness effects allowing for direct comparison with classical electric double layer theory and distinction of specifically adsorbed ions and electrostatically adsorbed ions.