Dynamic electrophoretic mobility of concentrated dispersions of spherical colloidal particles. On the consistent use of the cell model

This paper outlines a complete and self-consistent cell model theory of the electrokinetics of dense spherical colloidal suspensions for general electrolyte composition, frequency of applied field, zeta potential, and particle size. The standard electrokinetic equations, first introduced for any given particle configuration, are made tractable to computation by averaging over particle configurations. The focus of this paper is on the systematic development of suitable boundary conditions at the outer cell boundary obtained from global constraints on the suspension. The approach is discussed in relation to previously published boundary conditions that have often been introduced in an ad hoc manner. Results of a robust numerical calculation of high-frequency colloidal transport properties, such as dynamic mobility, using the present model are presented and compared with some existing dense suspension models.     

Electroviscous effect of moderately concentrated colloidal suspensions: Stern-layer influence

A previous model for the viscosity of moderately concentrated suspensions has been extended. The influence of a dynamic Stern layer (DSL), which produces an additional surface conductance at the electrolyte-particle interface, is included. The theoretical treatment is based on Happel's cell model with Simha's boundary conditions for the interparticle hydrodynamic interactions and on a dynamic Stern-layer model for ionic conduction on the particle surface according to Mangelsdorf and White (ref 39). The results are valid for arbitrary zeta potentials and double-layer thickness. Extensive theoretical predictions are shown and interesting new behaviors are found. The comparison with the results in the absence of additional surface conductance shows a great influence of this mechanism in the energy dissipation during the laminar flow of these suspensions. We conclude that the inclusion of a dynamic Stern layer will be required to match the predictions with the experimental results.     

Electrokinetics of concentrated suspensions of spherical colloidal particles with surface conductance, arbitrary zeta potential, and double-layer thickness in static electric fields

In this paper the electrophoretic mobility and the electrical conductivity of concentrated suspensions of spherical colloidal particles have been numerically studied under arbitrary conditions including zeta potential, particle volume fraction, double-layer thickness (overlapping of double layers is allowed), surface conductance by a dynamic Stern layer model (DSL), and ionic properties of the solution. We present an extensive set of numerical data of both the electrophoretic mobility and the electrical conductivity versus zeta potential and particle volume fraction, for different electrolyte concentrations. The treatment is based on the use of a cell model to account for hydrodynamic and electrical interactions between particles. Other theoretical approaches have also been considered for comparison. Furthermore, the study includes the possibility of adsorption and lateral motion of ions in the inner region of the double layers (DSL model), according to the theory developed by C. S. Mangelsdorf and L. R. White (J. Chem. Soc. Faraday Trans.86, 2859 (1990)). The results show that the correct limiting cases of low zeta potentials and thin double layers for dilute suspensions are fulfilled by our conductivity formula. Moreover, the presence of a DSL causes very important changes, even dramatic, on the values of both the electrophoretic mobility and the electrical conductivity for a great range of volume fractions and zeta potentials, specially when double layers of adjacent cells overlap, in comparison with the standard case (no Stern layer present). It can be concluded that in general the presence of a dynamic Stern layer causes the electrophoretic mobility to decrease and the electrical conductivity to increase in comparison with the standard case for every volume fraction, zeta potential, and double-layer thickness.     

The Electrophoretic Mobility and Electric Conductivity of a Concentrated Suspension of Colloidal Spheres with Arbitrary Double-Layer Thickness

The electrophoresis in a monodisperse suspension of dielectric spheres with an arbitrary thickness of the electric double layers is analytically studied. The effects of particle interactions are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations, which govern the ionic concentration distributions, the electric potential profile, and the fluid flow field in the electrolyte solution surrounding the charged sphere in a unit cell, are linearized assuming that the system is only slightly distorted from equilibrium. Using a perturbation method, these linearized equations are solved with the surface charge density (or zeta potential) of the particle as the small perturbation parameter. Analytical expressions for the electrophoretic mobility of the colloidal sphere in closed form correct to O(zeta) are obtained. Based on the solution of the electrokinetic equations in a cell, a closed-form formula for the electric conductivity of the suspension up to O(zeta(2)) is derived from the average electric current density. Comparisons of the results of the cell model with different conditions at the outer boundary of the cell are made for both the electrophoretic mobility and the electric conductivity. Copyright 2001 Academic Press.     

Colloid Vibration Potential in a Concentrated Suspension of Spherical Colloidal Particles

A relation between the dynamic electrophoretic mobility of spherical colloidal particles in a concentrated suspension and the colloid vibration potential (CVP) generated in the suspension by a sound wave is obtained from the analogy with the corresponding Onsager relation between electrophoretic mobility and sedimentation potential in concentrated suspensions previously derived on the basis of Kuwabara's cell model. The obtained expression for CVP is applicable to the case where the particle zeta potential is low, the particle relative permittivity is very small, and the overlapping of the electrical double layers of adjacent particles is negligible. It is found that CVP shows much stronger dependence on the particle volume fraction φ than predicted from the φ dependence of the dynamic electrophoretic mobility. It is also suggested that the same relation holds between the electrokinetic sonic amplitude of a concentrated suspension of spherical colloidal particles and the dynamic electrophoretic mobility. Copyright 1999 Academic Press.     

Analysis of the response of suspended colloidal soft particles to a constant electric field

A network model, originally designed for an electrokinetic study of soft particle suspensions, has been used for an in-depth analysis of the physical behavior of these systems under the action of an externally applied DC electric field. The versatility of the network simulation method used makes it possible to obtain information readily not only about the electrophoretic mobility, but also about any physical variable of interest at all points around the suspended particle: electric potential, ion concentrations, fluid velocity. The field-induced polarization of the double layer is described in terms of the dependence of these and other derived variables (volume charge density, electric field components, ion flux components) on the distance to the membrane-solution interface. In contrast to colloidal suspensions of hard particles, which basically depend on just two parameters (the reciprocal Debye length multiplied by the particle radius, kappaa, and the zeta potential, zeta), soft particle suspensions require a wider parameter set. First, there are two characteristic diffusion lengths in the system (one inside the membrane and the other in the solution) and two geometrical lengths (the core radius a and the membrane thickness (b-a)). Furthermore, there is the fixed charge density inside the membrane (and possibly a surface charge density over the core) that cannot be represented by a zeta potential. Finally, the parameter that characterizes the interaction between the fluid and the permeable membrane, gamma, strongly influences the behavior of the system. Dependences on all these parameters (except the geometrical ones) are included in this study.     

Diffusiophoresis in a suspension of charge-regulating colloidal spheres

An analytical study of diffusiophoresis in a homogeneous suspension of identical spherical charge-regulating particles with an arbitrary thickness of the electric double layers in a solution of a symmetrically charged electrolyte with a uniform prescribed concentration gradient is presented. The charge regulation due to association/dissociation reactions of ionogenic functional groups on the particle surface is approximated by a linearized regulation model, which specifies a linear relationship between the surface charge density and the surface potential. The effects of particle-particle electrohydrodynamic interactions are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations that govern the electric potential profile, the ionic concentration distributions, and the fluid flow field in the electrolyte solution surrounding the particle in a unit cell are linearized assuming that the system is only slightly distorted from equilibrium. Using a regular perturbation method, these linearized equations are solved with the equilibrium surface charge density (or zeta potential) of the particle as the small perturbation parameter. Closed-form formulas for the diffusiophoretic velocity of the charge-regulating sphere correct to the second order of its surface charge density or zeta potential are derived. Our results indicate that the charge regulation effect on the diffusiophoretic mobility is quite sensitive to the boundary condition for the electric potential specified at the outer surface of the unit cell. For the limiting cases of a very dilute suspension and a very thin or very thick electric double layer, the particle velocity is independent of the charge regulation parameter.     

Numerical calculation of the electrophoretic mobility of concentrated suspensions of soft particles

The electrophoretic mobility of spherical soft particles in concentrated colloidal suspensions is numerically calculated. The particle is modeled as a hard core coated with an ion-penetrable membrane bearing a uniform distribution of fixed charges, while the high particle concentration is taken into account by means of a cell model. The network simulation method used makes it possible to solve the problem without any restrictions on the values of the parameters such as particle concentration, membrane thickness, fixed charge density in the membrane, viscous drag in the membrane, number and valence of ionic species, electrolyte concentration, etc. The theoretical model used is similar to the one presented by Ohshima [H. Ohshima, J. Colloid Interface Sci. 225 (2000) 233], except for the use of the Shilov-Zharkikh, rather than the Levine-Neale, boundary condition for the electric potential, and the inclusion in the force balance equation of an additional term corresponding to the force exerted by the liquid on the core of the moving particle [J.J. López-García, C. Grosse, J. Horno, J. Colloid Interface Sci. 265 (2003) 327]. The obtained results only coincide with existing analytical expressions for low particle concentrations, low particle charge, and when the electrolyte concentration is high, the membrane is thick, and its resistance to the fluid flow is high. This suggests that most interpretations of the electrophoretic mobility of soft particles in concentrated suspensions require numerical calculations.     

Diffusiophoresis of concentrated suspensions of spherical particles with identical ionic diffusion velocities

Diffusiophoresis of concentrated suspensions of spherical particles subject to a small electrolyte gradient is analyzed theoretically at arbitrary levels of zeta potential and double-layer thickness. The Kuwabara unit cell model is adopted to describe the system under consideration. The effect of double-layer polarization is taken into account. It is found that the diffusiophoretic mobility exhibits a local maximum as well as a local minimum with varying zeta potential or double-layer thickness, similar to the corresponding dilute dispersion. The direction of the particle movement may even change back and forth. The previous low-zeta-potential approach is found to significantly overestimate the diffusiophoretic mobility as the zeta potential goes high. The deviation may be several fold sometimes. The effect of the volume fraction ratio of colloids is also examined. The higher the ratio, the lower the mobility.     

Diffusiophoresis in a suspension of spherical particles with arbitrary double-layer thickness

The diffusiophoresis in a homogeneous suspension of identical dielectric spheres with an arbitrary thickness of the electric double layers in a solution of a symmetrically charged electrolyte with a constant imposed concentration gradient is analytically studied. The effects of particle interactions (or particle volume fraction) are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations that govern the ionic concentration distributions, the electrostatic potential profile, and the fluid flow field in the electrolyte solution surrounding the charged sphere in a unit cell are linearized assuming that the system is only slightly distorted from equilibrium. Using a perturbation method, these linearized equations are solved with the surface charge density (or zeta potential) of the particle as the small perturbation parameter. Analytical expressions for the diffusiophoretic velocity of the dielectric sphere in closed form correct to the second order of its surface charge density or zeta potential are obtained from a balance between its electrostatic and hydrodynamic forces. Comparisons of the results of the cell model with different conditions at the outer boundary of the cell are made.     

Cell model calculation for electrokinetic phenomena in concentrated suspensions: an Onsager relation between sedimentation potential and electrophoretic mobility

Cell model calculations for the electrophoretic mobility, electrical conductivity and sedimentation potential in concentrated suspensions of colloidal particles with low zeta potentials are reviewed with particular emphasis on an Onsager relation between sedimentation potential and electrophoretic mobility. A general Onsager relation is derived on the basis of the thermodynamics of irreversible processes. This relation, which involves the ratio of the electrical conductivity K* of the suspension to the conductivity Kinfinity in the absence of the particles, reproduces the Onsager relation derived from cell model calculations at low zeta potentials, where K*/Kinfinity becomes (1 - phi)/(1 + phi/2), phi being the particle volume fraction.     

Numerical and analytical studies of the electrical conductivity of a concentrated colloidal suspension

In the past few years, different models and analytical approximations have been developed facing the problem of the electrical conductivity of a concentrated colloidal suspension, according to the cell-model concept. Most of them make use of the Kuwabara cell model to account for hydrodynamic particle-particle interactions, but they differ in the choice of electrostatic boundary conditions at the outer surface of the cell. Most analytical and numerical studies have been developed using two different sets of boundary conditions of the Neumann or Dirichlet type for the electrical potential, ionic concentrations or electrochemical potentials at that outer surface. In this contribution, we study and compare numerical conductivity predictions with results obtained using different analytical formulas valid for arbitrary zeta potentials and thin double layers for each of the two common sets of boundary conditions referred to above. The conductivity will be analyzed as a function of particle volume fraction, phi, zeta potential, zeta, and electrokinetic radius, kappaa (kappa(-1) is the double layer thickness, and a is the radius of the particle). A comparison with some experimental conductivity results in the literature is also given. We demonstrate in this work that the two analytical conductivity formulas, which are mainly based on Neumann- and Dirichlet-type boundary conditions for the electrochemical potential, predict values of the conductivity very close to their corresponding numerical results for the same boundary conditions, whatever the suspension or solution parameters, under the assumption of thin double layers where these approximations are valid. Furthermore, both analytical conductivity equations fulfill the Maxwell limit for uncharged nonconductive spheres, which coincides with the limit kappaa --> infinity. However, some experimental data will show that the Neumann, either numerical or analytical, approach is unable to make predictions in agreement with experiments, unlike the Dirichlet approach which correctly predicts the experimental conductivity results. In consequence, a deeper study has been performed with numerical and analytical predictions based on Dirichlet-type boundary conditions.     

Salt concentration and particle density dependence of electrophoretic mobilities of spherical colloids in aqueous suspension

Using laser Doppler velocimetry in the superheterodyne mode, we conducted a systematic study of the electrophoretic mobility of dispersions of small silica spheres (a=18 nm) suspended in water at different salinities and particle concentrations. The concentration of NaCl was varied from 40 microM up to 16 mM, while the particle concentrations were varied between 4.2x10(18) and 2.1x10(20) m-3. We find a decrease of mobility with increasing salt concentrations and an increase with increased particle number densities. The latter observation is not backed by the standard cell model of electrophoresis with Shilov-Zharkikh boundary conditions. Rather, if the experimental data are interpreted within that model, an unexpected change of the zeta potential at constant added salt concentration results. Interestingly, all experimental data collapse onto a single master curve, if plotted versus the ratio C* of particle counterions to added salt ions. We obtain a logarithmic increase of mobility for C*<1 and a plateau for C*>1. This may indicate a change of the Stern layer structure not yet included in the theoretical model.     

CE characterization of semiconductor nanocrystals encapsulated with amorphous silicium dioxide

The electrophoretic mobility of silica-encapsulated semiconductor nanocrystals (quantum dots) dependent on the pH and the ionic strength of the separation electrolyte has been determined by CE. Having shown the viability of the approach, the electrophoretic mobility mu of the nanoparticles investigated is calculated for varied zeta potential zeta, particle radius r, and ionic strength I employing an approximate analytical expression presented by Ohshima (J. Colloid Interface Sci. 2001, 239, 587-590). The comparison of calculated with measured data shows that the experimental observations exactly follow what would be expected from theory. Within the parameter range investigated at fixed zeta and I there is an increase in mu with r which is a nonlinear function. This dependence of mu on size parameters can be used for the size-dependent separation of particles. Modeling of mu as function of I and zeta makes it possible to calculate the size distribution of nanoparticles from electrophoretic data (using the peak shape of the particle zone in the electropherogram) without the need for calibration provided that zeta is known with adequate accuracy. Comparison of size distributions calculated via the presented method with size histograms determined from transmission electron microscopy (TEM) micrographs reveals that there is an excellent matching of the size distribution curves obtained with the two independent methods. A comparison of calculated with measured distributions of the electrophoretic mobility showed that the observed broad bands in CE studies of colloidal nanoparticles are mainly due to electrophoretic heterogeneity resulting from the particle size distribution.     

The Dynamic Mobility of Colloidal Semiconducting Antimony-Doped Tin(IV) Oxide Particles

The dynamic mobility spectra of suspensions of semiconducting tin(IV) oxide particles doped with antimony have been measured with the technique of electroacoustics. The magnitude of the complex mobility decreases essentially monotonically with increasing frequency, just as for a nonconducting (dielectric) particle under the same conditions. Unlike the case for a dielectric particle, however, the magnitudes at low frequency increase with increasing conductivity. The phase angle behavior is also different from that of a normal dielectric particle. The change in the phase angle behavior is most obvious at low suspension conductivity and high frequency where the phase angles showed a much smaller phase lag than at high conductivities. Reasonable agreement was found between the experimental mobility and the theoretical dynamic mobility spectra obtained with O'Brien's theory for the enhanced permittivity of semiconductors. Copyright 2001 Academic Press.     

Electrophoresis and electric conduction in a salt-free suspension of charged particles

The electrophoresis and electric conduction of a suspension of charged spherical particles in a salt-free solution are analyzed by using a unit cell model. The linearized Poisson-Boltzmann equation (valid for the cases of relatively low surface charge density or high volume fraction of the particles) and Laplace equation are solved for the equilibrium electric potential profile and its perturbation caused by the imposed electric field, respectively, in the fluid containing the counterions only around the particle, and the ionic continuity equation and modified Stokes equations are solved for the electrochemical potential energy and fluid flow fields, respectively. Explicit analytical formulas for the electrophoretic mobility of the particles and effective electric conductivity of the suspension are obtained, and the particle interaction effects on these transport properties are significant and interesting. The scaled zeta potential, electrophoretic mobility, and effective electric conductivity increase monotonically with an increase in the scaled surface charge density of the particles and in general decrease with an increase in the particle volume fraction, keeping each other parameter unchanged. Under the Debye-Hückel approximation, the dependence of the electrophoretic mobility normalized with the surface charge density on the ratio of the particle radius to the Debye screening length and particle volume fraction in a salt-free suspension is same as that in a salt-containing suspension, but the variation of the effective electric conductivity with the particle volume fraction in a salt-free suspension is found to be quite different from that in a suspension containing added electrolyte.     

Diffusiophoresis of a moderately charged spherical colloidal particle

An analytic expression is obtained for the diffusiophoretic mobility of a charged spherical colloidal particle in a symmetrical electrolyte solution. The obtained expression, which is expressed in terms of exponential integrals, is correct to the third order of the particle zeta potential so that it is applicable for colloidal particles with low and moderate zeta potentials at arbitrary values of the electrical double-layer thickness. This is an improvement of the mobility formula derived by Keh and Wei, which is correct to the second order of the particle zeta potential. This correction, which is related to the electrophoresis component of diffusiophoresis, becomes more significant as the difference between the ionic drag coefficients of electrolyte cations and anions becomes larger and vanishes in the limit of thin or thick double layer. A simpler approximate mobility expression is further obtained that does not involve exponential integrals.     

Dynamic electrophoretic mobility of spherical colloidal particles in salt-free concentrated suspensions

In this contribution, the dynamic electrophoretic mobility of spherical colloidal particles in a salt-free concentrated suspension subjected to an oscillating electric field is studied theoretically using a cell model approach. Previous calculations focusing the analysis on cases of very low or very high particle surface charge are analyzed and extended to arbitrary conditions regarding particle surface charge, particle radius, volume fraction, counterion properties, and frequency of the applied electric field (sub-GHz range). Because no limit is imposed on the volume fractions of solids considered, the overlap of double layers of adjacent particles is accounted for. Our results display not only the so-called counterion condensation effect for high particle charge, previously described in the literature, but also its relative influence on the dynamic electrophoretic mobility throughout the whole frequency spectrum. Furthermore, we observe a competition between different relaxation processes related to the complex electric dipole moment induced on the particles by the field, as well as the influence of particle inertia at the high-frequency range. In addition, the influences of volume fraction, particle charge, particle radius, and ionic drag coefficient on the dynamic electrophoretic mobility as a function of frequency are extensively analyzed.     

Electrophoretic mobility of a colloidal particle with constant surface charge density

When the electrophoretic mobility of a particle in an electrolyte solution is measured, the obtained electrophoretic mobility values are usually converted to the particle zeta potential with the help of a proper relationship between the electrophoretic mobility and the zeta potential. For a particle with constant surface charge density, however, the surface charge density should be a more characteristic quantity than the zeta potential because for such particles the zeta potential is not a constant quantity but depends on the electrolyte concentration. In this article, a systematic method that does not require numerical computer calculation is proposed to determine the surface charge density of a spherical colloidal particle on the basis of the particle electrophoretic mobility data. This method is based on two analytical equations, that is, the relationship between the electrophoretic mobility and zeta potential of the particle and the relationship between the zeta potential and surface charge density of the particle. The measured mobility values are analyzed with these two equations. As an example, the present method is applied to electrophoretic mobility data on gold nanoparticles (Agnihotri, S. M.; Ohshima, H.; Terada, H.; Tomoda, K.; Makino, K. Langmuir 2009, 25, 4804).