Surface Protolysis and Its Kinetics Impact the Electrical Double Layer.

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 ζ potential and electro-osmotic flow.

The importance of specifically adsorbed ions for electrokinetic phenomena: Bridging the gap between experiments and MD simulations.

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.

On the transferability of ion parameters to the TIP4P/2005 water model using molecular dynamics simulations.

Countless molecular dynamics studies have relied on available ion and water force field parameters to model aqueous electrolyte solutions. The TIP4P/2005 model has proven itself to be among the best rigid water force fields, whereas many of the most successful ion parameters were optimized in combination with SPC/E, TIP3P, or TIP4P/Ew water. Many researchers have combined these ions with TIP4P/2005, hoping to leverage the strengths of both parameter sets. To assess if this widely used approach is justified and to provide a guide in selecting ion parameters, we investigated the transferability of various commonly used monovalent and multivalent ion parameters to the TIP4P/2005 water model. The transferability is evaluated in terms of ion hydration free energy, hydration radius, coordination number, and self-diffusion coefficient at infinite dilution. For selected ion parameters, we also investigated density, ion pairing, chemical potential, and mean ionic activity coefficients at finite concentrations. We found that not all ions are equally transferable to TIP4P/2005 without compromising their performance. In particular, ions optimized for TIP3P water were found to be poorly transferable to TIP4P/2005, whereas ions optimized for TIP4P/Ew water provided nearly perfect transferability. The latter ions also showed good overall agreement with experimental values. The one exception is that no combination of ion parameters and water model considered here was found to accurately reproduce experimental self-diffusion coefficients. Additionally, we found that cations optimized for SPC/E and TIP3P water displayed consistent underpredictions in the hydration free energy, whereas anions consistently overpredicted the hydration free energy.