TY - JOUR
T1 - Molecular Origins of the Zeta Potential
AU - Předota, Milan
AU - Machesky, Michael L.
AU - Wesolowski, David J.
N1 - M.P. was supported by the Czech Science Foundation Project 13-08651S. Access to the CERIT-SC computing and storage facilities provided under the programme Center CERIT Scientific Cloud, part of the Operational Program Research and Development for Innovations, reg. no. CZ. 1.05/3.2.00/ 08.0144, is greatly appreciated. M.L.M. and D.J.W. were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences and Biosciences Division.
PY - 2016/10/11
Y1 - 2016/10/11
N2 - The zeta potential (ZP) is an oft-reported measure of the macroscopic charge state of solid surfaces and colloidal particles in contact with solvents. However, the origin of this readily measurable parameter has remained divorced from the molecular-level processes governing the underlying electrokinetic phenomena, which limits its usefulness. Here, we connect the macroscopic measure to the microscopic realm through nonequilibrium molecular dynamics simulations of electroosmotic flow between parallel slabs of the hydroxylated (110) rutile (TiO2) surface. These simulations provided streaming mobilities, which were converted to ZP via the commonly used Helmholtz-Smoluchowski equation. A range of rutile surface charge densities (0.1 to -0.4 C/m2), corresponding to pH values between about 2.8 and 9.4, in RbCl, NaCl, and SrCl2 aqueous solutions, were modeled and compared to experimental ZPs for TiO2 particle suspensions. Simulated ZPs qualitatively agree with experiment and show that "anomalous" ZP values and inequalities between the point of zero charge derived from electrokinetic versus pH titration measurements both arise from differing co- and counterion sorption affinities. We show that at the molecular level the ZP arises from the delicate interplay of spatially varying dynamics, structure, and electrostatics in a narrow interfacial region within about 15 Å of the surface, even in dilute salt solutions. This contrasts fundamentally with continuum descriptions of such interfaces, which predict the ZP response region to be inversely related to ionic strength. In reality the properties of this interfacial region are dominated by relatively immobile and structured water. Consequently, viscosity values are substantially greater than in the bulk, and electrostatic potential profiles are oscillatory in nature.
AB - The zeta potential (ZP) is an oft-reported measure of the macroscopic charge state of solid surfaces and colloidal particles in contact with solvents. However, the origin of this readily measurable parameter has remained divorced from the molecular-level processes governing the underlying electrokinetic phenomena, which limits its usefulness. Here, we connect the macroscopic measure to the microscopic realm through nonequilibrium molecular dynamics simulations of electroosmotic flow between parallel slabs of the hydroxylated (110) rutile (TiO2) surface. These simulations provided streaming mobilities, which were converted to ZP via the commonly used Helmholtz-Smoluchowski equation. A range of rutile surface charge densities (0.1 to -0.4 C/m2), corresponding to pH values between about 2.8 and 9.4, in RbCl, NaCl, and SrCl2 aqueous solutions, were modeled and compared to experimental ZPs for TiO2 particle suspensions. Simulated ZPs qualitatively agree with experiment and show that "anomalous" ZP values and inequalities between the point of zero charge derived from electrokinetic versus pH titration measurements both arise from differing co- and counterion sorption affinities. We show that at the molecular level the ZP arises from the delicate interplay of spatially varying dynamics, structure, and electrostatics in a narrow interfacial region within about 15 Å of the surface, even in dilute salt solutions. This contrasts fundamentally with continuum descriptions of such interfaces, which predict the ZP response region to be inversely related to ionic strength. In reality the properties of this interfacial region are dominated by relatively immobile and structured water. Consequently, viscosity values are substantially greater than in the bulk, and electrostatic potential profiles are oscillatory in nature.
UR - http://www.scopus.com/inward/record.url?scp=84991094784&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84991094784&partnerID=8YFLogxK
U2 - 10.1021/acs.langmuir.6b02493
DO - 10.1021/acs.langmuir.6b02493
M3 - Article
AN - SCOPUS:84991094784
SN - 0743-7463
VL - 32
SP - 10189
EP - 10198
JO - Langmuir
JF - Langmuir
IS - 40
ER -