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.     

Numerical study of colloidal suspensions of soft spherical particles using the network method. 1. DC electrophoretic mobility

The electrophoretic mobility of a spherical particle coated with a uniformly charged permeable membrane and suspended in a general electrolyte solution is calculated numerically. The network simulation method used makes it possible to solve the problem without any restrictions on the values of the parameters such as the membrane thickness, fixed charge density in the membrane, viscous drag in the membrane, number and valence of the ionic species, and electrolyte concentration. The theoretical model used is similar to the one presented by Ohshima (H. Ohshima, J. Colloid Interface Sci. 228 (2000) 190), except for 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. This inclusion is theoretically proven in the limiting case of a nonconducting suspending medium, in which the equation system can be analytically solved. The results obtained coincide with existing analytical expressions when the electrolyte concentration is high, the membrane is thick, and its resistance to the fluid flow is high.     

Electrophoresis of soft particles at high electrolyte concentrations: an interpretation by the Henry theory

The existence of electrophoretic mobility at high electrolyte concentrations defines a remarkable peculiarity in the electrosurface characteristics of soft particles. According to Ohshima [H. Ohshima, Colloids Surf. 103 (1995) 249], this effect is caused by the electroosmotic flow within the soft particle shell. An explanation supporting Ohshima's conclusion can be derived from classic electrokinetic theories. Based on the Henry theory [D.C. Henry, Proc. R. Soc. London Ser. A 133 (1931) 106], we demonstrate that the electrophoretic mobility of soft particles does not disappear at decinormal concentration.     

Influence of membrane layer properties on the electrophoretic behavior of a soft particle

The influence of the physical properties of the membrane layer of a soft particle, which comprises a rigid core and a porous membrane layer, on its electrophoretic behavior, is investigated. Because that influence was almost always neglected in the previous studies, the corresponding results can be unrealistic. The applicability of the model proposed is verified by the available theoretical and experimental results. The electrophoretic mobility of the particle under various conditions is simulated through varying the dielectric constant, the thickness, and the drag coefficient of the membrane layer, and the bulk ionic concentration. We show that under typical conditions, the deviation in the electrophoretic mobility arising from assuming that the dielectric constant of the membrane layer is the same as that of the bulk liquid phase can be in the order of 50%. In addition, the thicker the membrane layer and/or the higher the bulk ionic concentration, the larger the deviation. If the surface of the core of the particle is charged, as in the case of inorganic particles covered by an artificial membrane layer, the deviation at constant core surface potential is larger than that under other types of charged conditions. However, if the surface of the core is uncharged, as in the case of biocolloids, then that deviation becomes negligible. These findings are of fundamental significance to theoreticians in their analysis on the electrokinetic behaviors of soft particles, and to experimentalists in the interpretation of their data.     

Electrokinetic phenomena at grafted polyelectrolyte layers

During the last decades the electrokinetic theory of Smoluchowski (Z. Phys. Chem. 92 (1918) 129) was extended to be applicable for soft surfaces (grafted polyelectrolyte layers (PL), biological and artificial membranes, etc.) by either using the Debye approximation or numerical solutions. In the theory of Ohshima (Colloids Surf. A 103 (1995) 249) the nonlinearized Poisson-Boltzmann (PB) equation for thick and uniform PL is solved analytically and a general hydrodynamic equation is derived in an integral form. These advantages in the theory of Ohshima provided a base for the further development of a generalized electrokinetic theory for soft surfaces. In his theory the final equation for the electroosmotic (electrophoretic) velocity is specified for the case of the complete dissociation of ionic sites within PL. Accordingly, the equation may be used only if the difference between pK and pH is very large. However, it turned out that an analytical solution of the nonlinearized PB equation for thick PL is possible for any degree of dissociation. This was achieved using the approximation of excluded coions if the absolute value of the reduced Donnan potential is larger than 2 and due to the simplification in the case of weak dissociation, when the absolute value of the reduced Donnan potential is less than 2. Combining this generalized double layer (DL) theory for PL and the theory of Ohshima enables to obtain an analytical equation for electroosmosis for the general case of any degree of dissociation. This equation creates for the first time a theoretical base for the interpretation of electrokinetic fingerprinting (EF) for the characterization of soft surfaces.     

Electrophoresis of soft particles: analytic approximations

An approximate analytic expression is derived for the electrophoretic mobility of a weakly charged spherical soft particle (i.e., a hard particle covered with a weakly charged polyelectrolyte layer) on the basis of the general mobility expression for soft particles (Ohshima, H., J. Colloid Interface Sci. 2000, 228, 190-193). The obtained mobility expression, which reproduces various approximate results so far derived and gives some new mobility formulas, covers all types of weakly charged soft particles with arbitrary values of the thickness of polymer layer, the radius of the particle core, the electrophoretic softness, and the Debye length, including spherical polyelectrolytes with no particle core as well as spherical hard particles with no polyelectrolyte layer.     

Effects of dielectric gradients-mediated ions partitioning on the electrophoresis of composite soft particles: An analytical theory

In this work, we report original analytical expressions defining the electrophoretic mobility of composite soft particles comprising an inner core and a surrounding polymer shell with differentiated permeabilities to ions from aqueous background electrolyte and to fluid flow developed under applied DC field conditions. The existence of dielectric permittivity gradients operational at the core/shell and shell/solution interfaces is accounted for within the Debye–Hückel approximation and flat plate configuration valid in the thin double layer regime. The proposed electrophoretic mobility expressions, applicable to weakly to moderately charged particles with size well exceeding the Debye layer thickness, involve the relevant parameters describing the particle core/shell structure and the electrohydrodynamic features of the core and shell particle components. It is shown that the analytical expressions reported so far in literature for the mobility of hard (impermeable) or porous particles correspond to asymptotic limits of the more generic results detailed here. The impacts of dielectric-mediated effects of ions partitioning between bulk solution and particle body on the electrophoretic response are further discussed. The obtained expressions pave the way for a refined quantitative, analytical interpretation of electrophoretic mobility data collected on soft (nano)particles (e.g., functionalized dendrimers and multilayered polyelectrolytic particles) or biological cells (e.g., viruses) for which the classical hard core-soft shell representation is not appropriate.     

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.     

Effect of a Dynamic Stern Layer on the Sedimentation Velocity and Potential in a Dilute Suspension of Colloidal Particles

In this paper the theory of the sedimentation velocity and potential (gradient) in a dilute suspension of charged spherical colloidal particles developed by Ohshima et al. (H. Ohshima, T. W. Healy, L. R. White, and R. W. O'Brien, J. Chem. Soc., Faraday Trans. 2, 80, 1299 (1984)) has been modified to include the presence of a dynamic Stern layer on the particle surfaces. The starting point has been the theory that Mangelsdorf and White (C. S. Mangelsdorf, and L. R. White, J. Chem. Soc., Faraday Trans. 86, 2859 (1990)) developed to calculate the electrophoretic mobility of a colloidal particle allowing for the lateral motion of ions in the inner region of the double layer (dynamic Stern layer). The effects of varying the different Stern layer parameters on the sedimentation velocity and potential are discussed and compared to the case when a Stern layer is absent. The influence of electrolyte concentration and zeta potential of the particles is also analyzed. The results show that regardless of the chosen set of Stern layer and solution parameters, the presence of a dynamic Stern layer causes the sedimentation velocity to increase and the sedimentation potential to decrease, in comparison with the standard case (no Stern layer present). These changes are almost negligible when sedimentation velocity is concerned, but they are very important when it comes to the sedimentation potential. A justification for this fact can be given in terms of an Onsager reciprocal relation, connecting the magnitudes of the sedimentation potential and the electrophoretic mobility. As previously reported, the presence of a dynamic Stern layer exerts a great influence on the electrophoretic mobility of a colloidal particle, and by means of the Onsager relation, the same is confirmed to occur when the sedimentation potential is concerned. Copyright 2000 Academic Press.     

Electrokinetics of diffuse soft interfaces. 2. Analysis based on the nonlinear Poisson-Boltzmann equation

In a previous study (Langmuir 2004, 20, 10324), the electrokinetic properties of diffuse soft layers were theoretically investigated within the framework of the Debye-Hückel approximation valid in the limit of sufficiently low values for the Donnan potential. In the current paper, the electrokinetics is tackled on the basis of the rigorous nonlinearized Poisson-Boltzmann equation, the numerical evaluation of the electroosmotic velocity profile, and the analytically derived hydrodynamic velocity profile. The results are illustrated and discussed for a diffuse soft interface characterized by a linear gradient for the friction coefficient and the density of hydrodynamically immobile ionogenic groups in the transition region separating the bulk soft layer and the bulk electrolyte solution. In particular, it is shown how the strong asymmetry for the potential distribution, as met for high values of the bulk fixed charge density and/or low electrolyte concentrations, is reflected in the electrokinetic features of the diffuse soft layer. The analysis clearly highlights the shortcomings of the discontinuous approximation by Ohshima and others for the modeling of the friction and electrostatic properties of soft layers exhibiting high Donnan potentials. This is in line with reported electrokinetic measurements of various soft particles and permeable gels at low electrolyte concentrations which fail to match predictions based on Ohshima's theory.     

Electrophoresis of soft particles with hydrophobic inner core grafted with pH-regulated and highly charged polyelectrolyte layer

Electrophoresis of core–shell composite soft particles possessing hydrophobic inner core grafted with highly charged polyelectrolyte layer (PEL) has been studied analytically. The PEL bears pH-dependent charge properties due to the presence of zwitterionic functional groups. The dielectric permittivity of the PEL and bulk aqueous medium were taken to be different, which resulted in the ion-partitioning effect. Objective of this study was to provide a simple expression for the mobility of such core–shell soft particles under Donnan limit where the thickness of the PEL well exceeds the electric double layer thickness. Going beyond the widely used Debye–Hückel linearization, the nonlinear Poisson–Boltzmann equation coupled with Stokes–Darcy–Brinkman equations was solved to determine the electrophoretic mobility. The derived expression further recovers all the existing results for the electrophoretic mobility under various simplified cases. The graphical presentation of the results illustrated the impact of pertinent parameters on the electrophoretic mobility of such a soft particle.     

Transient electrophoresis of dielectric spheres

The dynamic electrophoretic response of a spherical dielectric particle suspended in an electrolyte solution to a step change in the applied electrics field is analytically studied. The electrical double layer surrounding the particle may have either a small but finite thickness or a very large thickness relative to the particle radius. For the case of electrophoresis of a particle with a thin double layer, the local electroosmotic velocity at the outer edge of the double layer evolving with time after the external field is imposed is used as an apparent slip boundary condition at the particle surface so that the unsteady equation of motion for the fluid flow outside the double layer is solved. Closed-form formulas for the transient electrophoretic mobility of the particle are derived as functions of relevant parameters. The results demonstrate that, when the double layer surrounding the particle is relatively thin, the normalized electrophoretic mobility at a given dimensionless time decreases monotonically with a decrease in the parameter kappaa, where kappa(-1) is the Debye screening length and a is the particle radius. When the double layer of the particle is relatively thick, the particle mobility can have magnitudes comparable to those for a particle with a thin double layer in the initial stage, but will become much smaller afterward. In general, the effect of the relaxation time for transient electrophoresis is negligible, regardless of the value of kappaa.     

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).     

Parameter identifiability in application of soft particle electrokinetic theory to determine polymer and polyelectrolyte coating thicknesses on colloids

Soft particle electrokinetic models have been used to determine adsorbed nonionic polymer and polyelectrolyte layer properties on nanoparticles or colloids by fitting electrophoretic mobility data. Ohshima first established the formalism for these models and provided analytical approximations ( Ohshima, H. Adv. Colloid Interface Sci.1995, 62, 189 ). More recently, exact numerical solutions have been developed, which account for polarization and relaxation effects and require fewer assumptions on the particle and soft layer properties. This paper characterizes statistical uncertainty in the polyelectrolyte layer charge density, layer thickness, and permeability (Brinkman screening length) obtained from fitting data to either the analytical or numerical electrokinetic models. Various combinations of particle core and polymer layer properties are investigated to determine the range of systems for which this analysis can provide a solution with reasonably small uncertainty bounds, particularly for layer thickness. Identifiability of layer thickness in the analytical model ranges from poor confidence for cases with thick, highly charged coatings, to good confidence for cases with thin, low-charged coatings. Identifiability is similar for the numerical model, except that sensitivity is improved at very high charge and permeability, where polarization and relaxation effects are significant. For some poorly identifiable cases, parameter reduction can reduce collinearity to improve identifiability. Analysis of experimental data yielded results consistent with expectations from the simulated theoretical cases. Identifiability of layer charge density and permeability is also evaluated. Guidelines are suggested for evaluation of statistical confidence in polymer and polyelectrolyte layer parameters determined by application of the soft particle electrokinetic theory.     

A simple model of the high-frequency dynamic mobility in concentrated suspensions

Because electroacoustic techniques are gaining interest in many fields of colloid science, a number of theories dealing with the phenomenon of electrophoresis in high-frequency (on the order of the MHz) electric fields have been developed. In the present work we propose a straightforward derivation of a simple formula for the dynamic mobility of colloidal particles in mildly concentrated systems. Starting with a simple expression for the electrophoretic mobility in dilute suspensions, given as a function of the zeta potential and of the dipole coefficient, we introduce successive corrections related to: (i) the back flow of fluid induced by the electrophoretic motion of the particles; (ii) the electrostatic interactions among particles; (iii) the difference between the macroscopic and the external electric fields; (iv) the difference between the zero-momentum and the laboratory reference frames. Considering furthermore that the frequency dependence of the dipole coefficient is due to the Maxwell-Wagner-O'Konski double-layer relaxation, we obtain a mobility expression that compares well with other (semi)analytical models and (in proper conditions) with numerical cell-model calculations. However, its main merit is that it allows to understand, to a large extent, the physical origin of the frequency and volume fraction dependences of the dynamic mobility.     

Effect of core charge density on the electrophoresis of a soft particle coated with polyelectrolyte layer

The electrophoresis of a charged soft particle with charged rigid core is considered under a weak imposed field condition. The rigid core of the soft particle is considered to have a finite dielectric permittivity and a fixed volume charge density. The electric potential distribution is determined by solving the Poisson-Boltzman equation out side the rigid core and a Poisson equation within the core along with continuity conditions on the core-shell interface. We have extended the analytic expression of Ohshima (Electrophoresis 27:526–533, 2006) for the electrophoretic mobility of a soft particle with a charged shell to include the effect of the volume charge density of the rigid core. Mobility based on the present expression matches exactly with the existing analytical solutions for a soft particle with an uncharged core. We have also made a comparison of our solution for mobility with an uncharged rigid core with the existing experimental results. The impact of the core charge density on the soft particle mobility is analyzed.     

Dynamic electrophoresis of a spherical dispersion of soft particles subject to a stress-jump condition

The dynamic mobility of a spherical dispersion of soft particles, where a particle comprises a rigid core and a membrane layer, is evaluated for the case when the shear stress across the membrane layer-liquid interface is discontinuous, the so-called stress-jump condition. We show that, due to the effect of double-layer deformation, the magnitude of the dynamic mobility of a particle has a local maximum and the corresponding phase angle has a negative (phase lead) local minimum at a low to medium level of the frequency of the applied electric field. This effect becomes insignificant if the frequency of the applied electric field is sufficiently high. The stress-jump condition may lead to a significant influence on the drag, and consequently the mobility of a particle. The degree of influence depends upon the sign of the stress-jump coefficient and the charged conditions of the membrane layer of the particle.     

Numerical calculation of the electrophoretic mobility of a spherical particle in a salt-free medium

By extending an approximate theory of the electrophoretic mobility of dilute spherical colloidal particles in a salt-free medium containing only counterions (H. Ohshima, J. Colloid Interface Sci. 248 (2002) 499--503), a systematic numerical method is given for the calculation of the electrophoretic mobility, which is based on an iteration method. We assume that each sphere is surrounded by a spherical free volume, within which counterions are distributed so that electro-neutrality is satisfied. The electrophoretic mobility is found to be determined mainly by the pressure due to the counterions at the outer surface of the free volume. It is shown how the mobility values deviate from those expected from Hückel's formula for high particle charges or zeta potentials because of the counterion condensation effect.     

Dielectric response of a concentrated colloidal suspension in a salt-free medium

In this paper the complex dielectric constant of a concentrated colloidal suspension in a salt-free medium is theoretically evaluated using a cell model approximation. To our knowledge this is the first cell model in the literature addressing the dielectric response of a salt-free concentrated suspension. For this reason, we extensively study the influence of all the parameters relevant for such a dielectric response: the particle surface charge, radius, and volume fraction, the counterion properties, and the frequency of the applied electric field (subgigahertz range). Our results display the so-called counterion condensation effect for high particle charge, previously described in the literature for the electrophoretic mobility, and also the relaxation processes occurring in a wide frequency range and their consequences on the complex electric dipole moment induced on the particles by the oscillating electric field. As we already pointed out in a recent paper regarding the dynamic electrophoretic mobility of a colloidal particle in a salt-free concentrated suspension, the competition between these relaxation processes is decisive for the dielectric response throughout the frequency range of interest. Finally, we examine the dielectric response of highly charged particles in more depth, because some singular electrokinetic behaviors of salt-free suspensions have been reported for such cases that have not been predicted for salt-containing suspensions.