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. 2. AC electrokinetic and dielectric properties

The network simulation method is used to solve numerically the equation system that determines the dynamic electrophoretic mobility and the dielectric response of dilute suspensions of soft particles. This system was extensively studied theoretically by Ohshima (H. Ohshima, J. Colloid Interface Sci. 233 (2001) 142-152), who obtained analytical expressions for the static and dynamic electrophoretic mobility. However, the validity of his analytical result is restricted to relatively thick membranes with high drag coefficient and to relatively high electrolyte concentrations. As for the dielectric properties, there are only a few works dealing with particles without a core (ion exchange resins) and, to our knowledge, no numerical studies. Our theoretical model is basically similar to Ohshima's, except that we take into account the mechanical force acting on the surface of the core, which he neglects. The inclusion of this term is crucial when the general problem including arbitrary values of the parameters is analyzed. However, it has little bearing when the membrane is thick and the drag coefficient is high, so that our results for the electrophoretic mobility generally confirm Ohshima's equation when all the required conditions are met.     

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.     

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.     

Specific cation adsorption on protein-covered particles and its influence on colloidal stability

Protein coated particles present an anomalous colloidal stability at high ionic strength when the classical theory (DLVO) predicts aggregation. This observed deviation from DLVO behaviour appears for electrolyte concentrations above some critical bulk value. As we have suggested in previous publications the existence of an additional short-range repulsive 'hydration force' due to specific hydrated cation adsorption could explain this anomalous stability. The overlap of the hydration layers when two particles approach should provoke this repulsive force. New evidence of this mechanism has been observed when electrophoretic mobilities of protein-carrying latex particles were measured at various concentrations of sodium and calcium chloride. In the latter case a sign reversal of zeta-potential was found, probably due to the specific adsorption of Ca(2+) ions on protein molecules. The adsorption increases with the medium pH. These results have been analyzed following the treatment proposed by Ohshima and co-workers for large charged colloidal particles coated with a layer of protein. This study shows an increase in the positive fixed-charge density on the protein caused by the adsorption of cations.     

Electrophoretic mobility of a highly charged colloidal particle in a solution of general electrolytes.

For a highly charged particle in an electrolyte solution, counterions are condensed very near the particle surface. The electrochemical potential of counterions accumulated near the particle surface is thus not affected by the applied electric field, so that the condensed counterions do not contribute to the particle electrophoretic mobility. In the present paper we derive an expression for the electrophoretic mobility mu(infinity) of a highly charged spherical particle of radius a and zeta potential zeta in the limit of very high zeta in a solution of general electrolytes with large ka (where k is the Debye-Hückel parameter) on the basis of our previous theory for the case of symmetrical electrolytes (H. Ohshima, J. Colloid Interface Sci. 263 (2003) 337). It is shown that zeta can formally be expressed as the sum of two components: the co-ion component, zetaco-ion, and the counterion component, zetacounterion (where zeta = zetaco-ion + zetacounterion) and that the limiting electrophoretic mobility mu(infinity) is given by mu(infinity) = epsilonr epsilon0 zetaco-ion(infinity)/eta + 0(1/ka), where zetaco-ion(infinity) is the high zeta-limiting form of zetaco-ion, epsilonr and eta are, respectively, the relative permittivity and viscosity of the solution, and epsilon0 is the permittivity of a vacuum. That is, the particle behaves as if its zeta potential were zetaco-ion(infinity), independent of zeta. For the case of a positively charged particle in an aqueous electrolyte solution at 25 degrees C, the value of zetaco-ion(infinity) is 35.6 mV for 1-1 electrolytes, 46.0 mV for 2-1 electrolytes, and 12.2 mV for 1-2 electrolytes. It is also found that the magnitude of mu(infinity) increases as the valence of co-ions increases, whereas the magnitude of mu(infinity) decreases as the valence of counterions increases.     

Electrokinetics of diffuse soft interfaces. 1. Limit of low Donnan potentials

The current theoretical approaches to electrokinetics of gels or polyelectrolyte layers are based on the assumption that the position of the very interface between the aqueous medium and the gel phase is well defined. Within this assumption, spatial profiles for the volume fraction of polymer segments (phi), the density of fixed charges in the porous layer (rho fix), and the coefficient modeling the friction to hydrodynamic flow (k) follow a step-function. In reality, the "fuzzy" nature of the charged soft layer is intrinsically incompatible with the concept of a sharp interface and therefore necessarily calls for more detailed spatial representations for phi, rho fix, and k. In this paper, the notion of diffuse interface is introduced. For the sake of illustration, linear spatial distributions for phi and rho fix are considered in the interfacial zone between the bulk of the porous charged layer and the bulk electrolyte solution. The corresponding distribution for k is inferred from the Brinkman equation, which for low phi reduces to Stokes' equation. Linear electrostatics, hydrodynamics, and electroosmosis issues are analytically solved within the context of streaming current and streaming potential of charged surface layers in a thin-layer cell. The hydrodynamic analysis clearly demonstrates the physical incorrectness of the concept of a discrete slip plane for diffuse interfaces. For moderate to low electrolyte concentrations and nanoscale spatial transition of phi from zero (bulk electrolyte) to phi o (bulk gel), the electrokinetic properties of the soft layer as predicted by the theory considerably deviate from those calculated on the basis of the discontinuous approximation by Ohshima.     

Modelling of pressure-driven membrane separation of electrolytes.: `High temperature' approximation

A model of pressure-driven membrane process of electrolyte separation is presented. The electric field potential assumed as being known, exact solution for permeate composition is readily obtained. All species are assumed to have the same convection velocity. Local electroneutrality condition is not used. The electric potential has been taken into account under high temperature approximation, thus reducing the problem to algebraic equation in exp(Ψ), where Ψ is dimensionless flow potential, and making it possible to calculate concentrations of ions in permeate. Negative retention is shown to be possible for one-component electrolyte solution. For electrolyte mixtures, concentration of ion with high charge is shown to “govern” the membrane selectivity in respect to low-charge ions. Results obtained are in qualitative accordance with the earlier experimental data on membrane separation of reaction mixtures in homogeneous catalysis.     

Surface Charge Density/Surface Potential Relationship for a Cylindrical Particle in an Electrolyte Solution

An accurate analytic expression of the surface charge density/surface potential relationship for an infinitely long cylindrical colloidal particle in a solution of general electrolytes is derived from an approximate solution to the nonlinear cylindrical Poisson– Boltzmann equation. The mathematical procedure is based on a method developed previously by Ohshima, Healy, and White for the case of a sphere (J. Colloid Interface Sci.90, 17 (1982)). Comparison is made with exact numerical results. Accurate expressions for the potential distribution around a cylinder and the effective surface potential of a cylinder are also derived. Finally, expressions for the double-layer interaction energy and force between two cylinders at large separations are derived on the basis of the method of Brenner and Parsegian (Biophys. J.14, 327 (1974)).     

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

Interaction force between two charged plates immersed in a solution of charged particles. Coupling between double layer and depletion forces

When two parallel plates are immersed in a solution of small charged particles, the center of the particles is excluded from a region of thickness D/2 near the plate, where D is their diameter. The approach which Langmuir developed for the double layer repulsion in the presence of an electrolyte with ions of negligible size is extended to the case in which one of the "ions" is a charged particle of finite, relatively small size. A general expression for the force generated between the two charged plates immersed in an electrolyte solution containing relatively small charged particles is derived. In this expression, only the electrical potential at the middle distance between the plates is required to calculate the force. A Poisson-Boltzmann equation which accounts for the volume exclusion of the charged particles in the vicinity of the surface is solved to obtain the electrical potential at the middle between the two plates. Starting from this expression, some results obtained previously for the depletion force acting between two plates or two spheres are rederived. For charged plates immersed in a solution of an electrolyte and charged small particles, the effects of the particle charge, particle charge sign, particle size, and volume fraction of the particles on the force acting between the two plates are examined.     

Magnetohydrodynamic effects on a charged colloidal sphere with arbitrary double-layer thickness

An analytical study is presented for the magnetohydrodynamic (MHD) effects on a translating and rotating colloidal sphere in an arbitrary electrolyte solution prescribed with a general flow field and a uniform magnetic field at a steady state. The electric double layer surrounding the charged particle may have an arbitrary thickness relative to the particle radius. Through the use of a simple perturbation method, the Stokes equations modified with an electric force term, including the Lorentz force contribution, are dealt by using a generalized reciprocal theorem. Using the equilibrium double-layer potential distribution from solving the linearized Poisson-Boltzmann equation, we obtain closed-form formulas for the translational and angular velocities of the spherical particle induced by the MHD effects to the leading order. It is found that the MHD effects on the particle movement associated with the translation and rotation of the particle and the ambient fluid are monotonically increasing functions of κa, where κ is the Debye screening parameter and a is the particle radius. Any pure rotational Stokes flow of the electrolyte solution in the presence of the magnetic field exerts no MHD effect on the particle directly in the case of a very thick double layer (κa→0). The MHD effect caused by the pure straining flow of the electrolyte solution can drive the particle to rotate, but it makes no contribution to the translation of the particle.     

Suspended particles surrounded by an inhomogeneously charged permeable membrane. Solution of the Poisson-Boltzmann equation by means of the network method

The Poisson-Boltzmann equation is numerically solved for a suspended spherical particle surrounded by a permeable membrane that contains an inhomogeneous distribution of fixed charges. The calculations are carried out using the network simulation method, which makes it possible to solve the problem in the most general case, extending previous results (J.P. Hsu, Y.C. Kuo, J. Membrane Sci. 108 (1995) 107). Approximate analytical expressions for the counterion concentration and the electric potential in the membrane are also presented, together with criteria that determine their ranges of validity. The limiting case of a distribution of fixed charges in the membrane that reduces to a surface charge is also analyzed. It is shown that the solution for this case, considering a vanishingly small radius of the core, reduces to a superposition of solutions corresponding to a charged impermeable particle suspended in an electrolyte solution and to a cavity filled with a charged electrolyte solution.     

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.     

Electrophoretic motion of a spherical particle with a symmetric nonuniform surface charge distribution in a nanotube

The electrophoretic motion of a spherical nanoparticle, subject to an axial electric field in a nanotube filled with an electrolyte solution, has been investigated using a continuum theory, which consists of the Nernst-Planck equations for the ionic concentrations, the Poisson equation for the electric potential in the solution, and the Stokes equation for the hydrodynamic field. In particular, the effects of nonuniform surface charge distributions around the nanoparticle on its axial electrophoretic motion are examined with changes in the bulk electrolyte concentration and the surface charge of the tube's wall. A particle with a nonuniform charge distribution is shown to induce a corresponding complex ionic concentration field, which in turn influences the electric field and the fluid motion surrounding the particle and thus its electrophoretic velocity. As a result, contrary to the relatively simple dynamics of a particle with a uniform surface charge, dominated by the irradiating electrostatic force, that with a nonuniform surface charge distribution shows various intriguing behaviors due to the additional interplay of the nonuniform electro-osmotic effects.     

Contribution of convection,diffusion and migration to electrolyte transport through nanofiltration membranes

Transport mechanisms through nanofiltration membranes are investigated in terms of contribution of convection, diffusion and migration to electrolyte transport. A Donnan steric pore model, based on the application of the extended Nernst-Planck equation and the assumption of a Donnan equilibrium at both membrane-solution interfaces, is used. The study is focused on the transport of symmetrical electrolytes (with symmetric or asymmetric diffusion coefficients). The influence of effective membrane charge density, permeate volume flux, pore radius and effective membrane thickness to porosity ratio on the contribution of the different transport mechanisms is investigated. Convection appears to be the dominant mechanism involved in electrolyte transport at low membrane charge and/or high permeate volume flux and effective membrane thickness to porosity ratio. Transport is mainly governed by diffusion when the membrane is strongly charged, particularly at low permeate volume flux and effective membrane thickness to porosity ratio. Electromigration is likely to be the dominant mechanism involved in electrolyte transport only if the diffusion coefficient of coions is greater than that of counterions.     

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.     

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.     

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 with charged rigid core coated with pH-regulated polyelectrolyte layer

We present a theoretical study on the electrophoresis of a soft particle with a dielectric charged rigid core grafted with a charge-regulated polyelectrolyte layer. The polyelectrolyte layer possesses either an acidic or a basic functional group and the charge dissociation depends on the local pH and ionic concentration of the electrolyte. The dielectric rigid core is considered to possess a uniform volumetric charge density. The electric potential distribution is determined by computing the Poisson-Boltzmann equation outside the core coupled with a Poisson equation inside the impermeable core along with suitable matching conditions at the core-shell interface. The computed electric field is used to determine the mobility of the particle through an existing analytic expression based on the Debye-Huckel approximation. Our results are found to be in good agreement with the existing solutions for the limiting cases. The influence of the core charge density, ionic concentration, and pH of the electrolyte on the particle mobility is studied for different choice of hydrodynamic penetration length of the polyelectrolyte and dissociation constant of the functional group. The critical value of the pH required to achieve zero mobility is estimated. We find that in a monovalent electrolyte solution, the soft particle with a net negative (positive) charge can have positive (negative) mobility.