Next Article in Journal
On the Adequacy of Some Low-Order Moments Method to Simulate Certain Particle Removal Processes
Previous Article in Journal
AFM Slip Length Measurements for Water at Selected Phyllosilicate Surfaces
Previous Article in Special Issue
Effect of Produced Sand Particles and Fines on Scale Inhibitor: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electrophoretic Mobility and Electric Conductivity of Salt-Free Suspensions of Charged Soft Particles

Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
*
Author to whom correspondence should be addressed.
Colloids Interfaces 2021, 5(4), 45; https://doi.org/10.3390/colloids5040045
Submission received: 27 September 2021 / Revised: 11 October 2021 / Accepted: 14 October 2021 / Published: 17 October 2021
(This article belongs to the Special Issue Locomotion of Colloidal Particles)

Abstract

:
A unit cell model is employed to analyze the electrophoresis and electric conduction in a concentrated suspension of spherical charged soft particles (each is a hard core coated with a porous polyelectrolyte layer) in a salt-free medium. The linearized Poisson–Boltzmann equation applicable to a unit cell is solved for the equilibrium electrostatic potential distribution in the liquid solution containing the counterions only surrounding a soft particle. The counterionic continuity equation and modified Stokes/Brinkman equations are solved for the ionic electrochemical potential energy and fluid velocity distributions, respectively. Closed-form formulas for the electrophoretic mobility of the soft particles and effective electric conductivity of the suspension are derived, and the effect of particle interactions on these transport characteristics is interesting and significant. Same as the case in a suspension containing added electrolytes under the Debye–Hückel approximation, the scaled electrophoretic mobility in a salt-free suspension is an increasing function of the fixed charge density of the soft particles and decreases with increases in the core-to-particle radius ratio, ratio of the particle radius to the permeation length in the porous layer, and particle volume fraction, keeping the other parameters unchanged. The normalized effective electric conductivity of the salt-free suspension also increases with an increase in the fixed charge density and with a decrease in the core-to-particle radius ratio, but is not a monotonic function of the particle volume fraction.

1. Introduction

When charged particles are suspended in an ionic solution and subjected to an external electric field, both the particles and their neighboring ions move owing to electrophoresis and electric migration (together with diffusion), respectively. As a result, the disperse fluid is dragged to flow, and the electric current in the suspension is different from that in a bulk solution. Numerous formulas for the electrophoretic mobility and electric conductivity in suspensions of charged particles have been obtained in the past [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16].
Salt-free solutions are liquid media that contain no added electrolytes, but only counterions dissociated from ionogenic groups at adjoining solid surfaces [17,18,19]. In a salt-free solution suspending colloidal particles, most of the counterions are condensed in the electric double layers when the fixed charge density of the particles is high, whereas this effect of counterion condensation and electrostatic screening disappears when the fixed charge density is low. Thus, the electrophoresis and electric conduction in a salt-free suspension of particles in general differ from those in salt-containing suspensions.
The electrophoresis and electric conduction of dilute salt-free suspensions of charged hard and soft spheres were analyzed by Ohshima [20,21,22] using a unit cell model and the electrophoretic mobility and electric conductivity in the suspensions are found to be proportional to the fixed charges of the particles (and coincide with those in salt-containing suspensions) when these charges are lower than some critical values, but approach constants due to the effect of counterion condensation around the fixed charges when they are higher than the critical values. On the other hand, Carrique et al. [23] numerically solved for the electrophoretic mobility and electric conductivity of concentrated salt-free suspensions of hard spherical particles via the unit cell model and confirmed the effect of counterion condensation at high surface charge density of the particles. Chiang et al. [24] and He et al. [25] also numerically predicted the electrophoretic mobility of concentrated salt-free suspensions of charged hard and porous spheres, respectively, using the cell model. Recently, the electrophoresis and electric conduction of salt-free suspensions of hard spherical particles have been analyzed by Luo and Keh [26] using the cell model to investigate the effects of particle concentration. However, the particle concentration effects on the electrophoretic mobility and electric conductivity in a salt-free suspension of charged soft particles [27,28] of arbitrary volume fraction have not been analytically studied.
In this article, the electrophoresis and electric conduction in concentrated suspensions of spherical charged soft particles (polyelectrolyte-coated particles) in salt-free solutions are analyzed. The linearized Poisson–Boltzmann equation, counterionic continuity equation, and modified Stokes/Brinkman equations are solved for the electrostatic potential, ionic electrochemical potential energy, and fluid flow profiles, respectively. Explicit formulas for the electrophoretic velocity and effective electric conductivity are obtained in Equations (45) and (47), respectively, in terms of the dimensionless fixed charge density and concentration of the particles and other relevant parameters.

2. Analysis

Consider a suspension of charged soft spheres of radius a in a salt-free medium containing only counterions with the valence Z . Each charged soft sphere refers to an uncharged dielectric hard core of radius r 0 coated with a uniform solution-permeable charged porous layer (e.g., polyelectrolyte layer) of thickness a r 0 , and it reduces to an entirely charged porous sphere of radius a in the limit r 0 = 0 . When the suspension is subject to an imposed electric field E e z , the particles translate with a velocity U e z by electrophoresis, where e z is the unit vector in the z direction. We use a unit cell model, as shown in Figure 1, in which each particle is encompassed by a concentric solution shell of outside radius b and the particle-to-cell volume ratio equals the particle volume fraction φ of the whole suspension, i.e., φ = ( a / b ) 3 . The origin of the spherical coordinates ( r ,   θ ,   φ ) is set at the center of the particle/cell, and the problem is axially symmetric about the z axis.
To obtain the electrophoretic mobility of the soft particles and effective electric conductivity of the salt-free suspension, we first need to determine the distributions of the electrostatic potential, fluid velocity, and counterionic electrochemical potential energy in a unit cell.

2.1. Electrostatic Potential Profile

The electrostatic potential distribution ψ ( r , θ ) in the ionic solution between the hard core of the soft particle and the cell border ( r 0 r b ) can be expressed as the equilibrium potential profile ψ eq ( r ) arising from the fixed charges in the porous layer ( r 0 r a ) of the particle together with counterions in the ambient fluid superimposed by the perturbed potential profile ψ a ( r , θ ) resulting from the applied electric field E e z [6],
ψ = ψ eq + ψ a .
With adopting the Debye–Hückel approximation at equilibrium, the electrostatic potential ψ eq satisfies the linearized Poisson–Boltzmann equation [13],
1 r 2 d d r ( r 2 d ψ eq d r ) = k T Z e κ 2 exp ( Z e ψ eq k T ) h ( r ) Q ε κ 2 ( k T Z e + ψ eq ) h ( r ) Q ε ,
where ε is the dielectric permittivity of the fluid, Q is the fixed-charge density in the porous layer of the soft particle, κ = ( n av / ε k T ) 1 / 2 Z e is the Debye screening parameter, h ( r ) is a step function equal to unity if r 0 r a and zero otherwise, k is Boltzmann’s constant, T is the absolute temperature, e is the elementary electric charge, and n av is the average concentration (number density) of the counterions in the salt-free solution, which is related to the fixed-charge density in the porous layer of the soft particle by Q = ( b 3 r 0 3 ) Z e n av / ( a 3 r 0 3 ) with the electric neutrality in a unit cell containing a particle and its ambient counterions only. In Equation (2), the value of ψ eq has been set to be zero at the position where the local ionic concentration is identical to n av . The boundary conditions for the equilibrium potential, which allow overlapping of the neighboring electric double layers, are
r = r 0 ,   b : d ψ eq d r = 0 ,
r = a : ψ eq   and   d ψ eq d r   are   continuous .
The solution to Equations (2)–(4) can be obtained as
ψ eq ( r ) = Q ε κ 2 [ g 1 ( r 0 ) g 2 + ( r ) 2 g 0 ( r 0 ) κ r a 3 r 0 3 b 3 r 0 3 ]   for   a r b ,
ψ eq ( r ) = Q e κ ( 2 b a + 2 r 0 r ) 2 g 0 ( r 0 ) ε κ 3 r { 1 g 0 ( r 0 ) ( κ b 1 ) [ κ a ( κ r 0 1 ) + κ r 0 ] + κ b
+ ( κ r 0 1 ) ( κ b + 1 ) e κ ( a 2 b ) [ ( κ a 1 ) e κ a 2 κ r e κ r ]
( κ r 0 + 1 ) ( κ b 1 ) e κ ( r 2 r 0 ) [ ( κ a + 1 ) e κ r 2 κ r e κ a ]
+ ( κ r 0 + 1 ) ( κ a 1 ) ( κ b + 1 ) e 2 κ ( a b r 0 + r ) } Q ( a 3 r 0 3 ) ε κ 2 ( b 3 r 0 3 )   for   r 0 r a ,
where
g 0 ( r ) = ( κ b 1 ) ( κ r + 1 ) e 2 κ b ( κ b + 1 ) ( κ r 1 ) e 2 κ r ,
g 1 ( r ) = ( κ a 1 ) ( κ r + 1 ) ( κ a + 1 ) ( κ r 1 ) e 2 κ ( r a ) ,
g 2 ± ( r ) = ( κ b + 1 ) e κ ( a + r ) ± ( κ b 1 ) e κ ( a + 2 b r ) .
The relative surface potential of the particle can be expressed as ψ s = ψ eq ( a ) ψ eq ( b ) , the potential difference across the fluid shell in the unit cell, which is positive. Equation (5) results in
Z e ψ s k T = Q ¯ g 1 ( r 0 ) 2 ( κ a ) 3 g 0 ( r 0 ) [ g 2 + ( a ) 2 κ a e κ ( a + b ) ] ,
where the dimensionless fixed charge density of the porous layer of the soft particle defined by Q ¯ = a 2 Z e Q / ε k T (which equals ( κ a ) 2 ( a 3 r 0 3 φ ) / ( a 3 r 0 3 ) φ with the electric neutrality in a unit cell containing a particle and its ambient counterions) is also positive ( Z has the same sign as Q ) but low (less than about 10 and valid for the Debye–Hückel approximation). Note that only two of the three parameters Q ¯ , φ , and κ a in a salt-free suspension may be specified independently for a given value of the core-to-particle radius ratio r 0 / a , unlike suspensions containing added electrolytes in which all of them are independent parameters tunable in experiments. Equation (10) predicts that, for a nonzero value of φ , Z e ψ s / k T increases with an increase in Q ¯ from zero at Q ¯ = 0 .
The electric potential ψ a induced by the applied electric field E e z satisfies the Laplace equation
2 ψ a = 0
and boundary conditions
r = r 0 : ψ a r = 0 ,
r = b : ψ a = E r cos θ   ( for   the   Dirichlet   approach ) ,
ψ a r = E cos θ   ( for   the   Neumann   approach ) .
Equation (12) assumes that the permittivity of the hard core of the soft particle is small relative to that of the fluid. The solution to Equations (11)–(14) is
ψ a = E 2 ν ( 2 + r 0 3 r 3 ) r cos θ ,
where ν = 1 + r 0 3 / 2 b 3 if Equation (13) is employed and ν = 1 r 0 3 / b 3 if Equation (14) is used. Note that both boundary conditions (13) and (14) lead to the same result of ν = 1 for a suspension of porous spheres ( r 0 = 0 ). Additionally, the electric field strength at infinity E in Equation (13) might be replaced more appropriately by the volume-mean electric field strength [11].

2.2. Ionic Electrochemical Potential Energy Profile

The perturbed electrochemical potential energy distribution μ ( r , θ ) of the counterions, which will be used to calculate the effective electric conductivity of the suspension of particles, satisfies the continuity equation [6]
2 μ = Z 2 e 2 E ν k T ( 1 r 0 3 r 3 ) d ψ eq d r cos θ
and boundary conditions
r = r 0 : μ r = 0 ,
r = a :   μ   and   μ r   are   continuous ,
r = b : μ = Z e E r cos θ   ( if   Equation   ( 13 )   is   used ) ,
μ r = Z e E cos θ   ( if   Equation   ( 14 )   is   used ) .
The solution to Equations (16)–(20) is
μ = Z e E ν [ 1 + r 0 3 2 r 3 + F μ ( r ) ] r cos θ ,
where
F μ ( r ) = 1 6 r 3 { [ 2 r 0 3 I 3 ( a , b ) + b 3 I 0 ( a , b ) ] 2 r 3 + r 0 3 ν b 3 + [ 2 I 3 ( r 0 , a ) + I 0 ( r 0 , a ) ]
× 2 ( 1 ν ) r 3 + r 0 3 ν + 2 r 0 3 I 3 ( a , r ) 2 r 3 I 0 ( a , r ) }   for   a r b ,
F μ ( r ) = 1 6 r 3 { [ 2 r 0 3 I 3 ( r 0 , b ) + b 3 I 0 ( r 0 , b ) ] 2 r 3 + r 0 3 ν b 3
+ 2 r 0 3 I 3 ( r 0 , r ) 2 r 3 I 0 ( r 0 , r ) }   for   r 0 r a ,
I n ( r 1 , r 2 ) = Z e k T r 1 r 2 ( 1 r 0 3 r 3 ) ( r r 0 ) n d ψ eq d r d r ,
and ψ eq was obtained in Equations (5) and (6).

2.3. Fluid Flow Field

The fluid velocity profile v ( r , θ ) and dynamic pressure distribution p ( r , θ ) satisfy the modified Stokes/Brinkman equations
2 v h ( r ) λ 2 v = 1 η [ p ε κ 2 ( k T Z e + ψ eq ) ψ a ] ,
v = 0 ,
where λ is the reciprocal of the permeation length (shielding coefficient) in the porous layer of the soft particle and η is the fluid viscosity. Letting the coordinate frame travel with the particle, one may write the boundary conditions of the fluid velocity as
r = r 0 : v r = v θ = 0 ,
r = a : v r ,   v θ ,   τ r r ,   and   τ r θ   are   continuous ,
r = b : v r = U cos θ ,
τ r θ η = r r ( v θ r ) + 1 r v r θ = 0   ( for   the   Happel   model ) ,
[ × v ] φ = 1 r r ( r v θ ) 1 r v r θ = 0   ( for   the   Kuwabara   model ) ,
where v r and v θ are the nontrivial fluid velocity components, τ r r and τ r θ are the normal and shear stresses, respectively, and U is the electrophoretic velocity of the particle to be determined.
The solution to Equations (25)–(31) and Equation (15) is
v r = Q E ν κ 2 η F r ( r ) cos θ ,
v θ = ( r 2 v r ) 2 r r tan θ ,
p = Q E ν κ 2 a F p ( r ) cos θ ,
where
F r ( r ) = C 1 J 2 ( r ) + [ C 2 + J 3 ( r ) ] a r + [ C 3 1 5 J 5 ( r ) ] ( a r ) 3
+ [ C 4 + 1 5 J 0 ( r ) ] ( r a ) 2   for   a r b ,
F r ( r ) = C 5 + [ C 6 + C 7 α ( λ r ) + C 8 β ( λ r ) ] ( a r ) 3 2 ( λ a ) 2 [ J 0 ( r ) J 3 ( r ) ( a r ) 3
3 J α ( r ) β ( λ r ) ( λ r ) 3 + 3 J β ( r ) α ( λ r ) ( λ r ) 3 ] for   r 0 r a ,
F p ( r ) = [ C 2 + J 3 ( r ) ] ( a r ) 2 + { 10 C 4 + 2 J 0 ( r )
ε a 2 κ 4 Q ( 1 + r 0 3 2 r 3 ) [ k T Z e + ψ eq ( r ) ] } r a for   a r b ,
F p ( r ) = J 3 ( r ) ( a r ) 2 + { 2 J 0 ( r ) ( λ a ) 2 [ C 5 C 6 2 ( a r ) 3 ]
ε a 2 κ 4 Q ( 1 + r 0 3 2 r 3 ) [ k T Z e + ψ eq ( r ) ] } r a for   r 0 r a ,
J n ( r ) = Z e ( a 3 r 0 3 ) 6 k T ( b 3 r 0 3 ) ( κ a ) 2 a r ( r a ) n ( 2 + r 0 3 r 3 ) d ψ eq d r d r ,
J α ( r ) = Z e ( a 3 r 0 3 ) 6 k T ( b 3 r 0 3 ) ( κ a ) 2 a r α ( λ r ) ( 2 + r 0 3 r 3 ) d ψ eq d r d r ,
J β ( r ) = Z e ( a 3 r 0 3 ) 6 k T ( b 3 r 0 3 ) ( κ a ) 2 a r β ( λ r ) ( 2 + r 0 3 r 3 ) d ψ eq d r d r ,
α ( x ) = x cosh x sinh x ,
β ( x ) = x sinh x cosh x ,
The expressions for the dimensionless constants C n in both the Happel and Kuwabara cell models are lengthy and available elsewhere [29], and the equilibrium potential ψ eq is given by Equations (5) and (6).

2.4. Electrophoretic Velocity

The drag force exerted by the fluid on the outer border of a unit cell is
F h = 2 π b 2 0 π { p e r + η [ v + ( v ) T ] e r } r = b sin θ d θ ,
where e r is the unit vector in the r direction. As an entire cell is electrically neutral, this drag force vanishes. Applying this constraint after the substitution of Equations (32)–(38) into Equation (44), we obtain the particle velocity as
U = Q a 2 E ν η { J 2 ( b ) ( κ a ) 2 + 1 ( κ a ) 2 A 0 [ A 1 ε κ 4 6 a Q { ( 2 b 3 + r 0 3 ) [ k T Z e + ψ eq ( b ) ] } + A 2 J 0 ( b )
+ ( A 1 + A 3 ) J 3 ( b ) + A 4 J 5 ( b ) + A 5 J 0 ( r 0 ) + A 6 J 3 ( r 0 ) + A 7 J α ( r 0 ) + A 8 J β ( r 0 ) ] } ,
and the expressions of the dimensionless constants A n are lengthy and available elsewhere [29]. Since the perturbation to the electric potential distribution ψ ( r , θ ) by the fluid flow is not considered in the analysis, the relaxation (polarization) effect of the mobile ions around the soft sphere is not contained in this result. For a suspension of charged porous spheres ( r 0 = 0 and ν = 1 ), Equation (45) reduces to Equation (A1) or (A2) in Appendix A.

2.5. Electric Conductivity

The volume average of the electric current density in a suspension of charged spherical particles leads to its effective electric conductivity, expressed by [9]
Λ = Λ + 3 Z e D n av 2 b k T E 0 π (   r μ r μ ) r = b sin θ cos θ d θ ,
where Λ = Z 2 e 2 n av D / k T is a characteristic electric conductivity of the counterionic solution, which is an increasing function of the particle volume fraction φ for a specified value of the dimensionless fixed charge density Q ¯ , and D is the diffusion coefficient of the counterions. Substituting Equations (21)–(23) for the perturbed ionic electrochemical potential energy μ as well as Equations (5) and (6) for the equilibrium potential ψ eq into Equation (46), we obtain the effective electric conductivity as
Λ = Λ [ 1 3 r 0 3 2 ν b 3 + b ν ( d F μ d r ) r = b ] = Λ { 1 3 r 0 3 2 ν b 3
+ Q ¯   32 ν 2 g 0 ( r 0 ) κ 5 a 6 b 6 [ b 3 g 1 ( b ) K 1 8 a 3 g 1 ( r 0 ) K 2 ] } ,
where
K 1 = κ 6 r 0 6 a 4 e κ a { e 2 κ r 0 ( κ r 0 1 ) [ E ( κ a ) E ( κ r 0 ) ] + ( κ r 0 + 1 ) [ E ( κ a ) E ( κ r 0 ) ] }
+ 2 e κ ( a + r 0 ) { 2 κ 3 r 0 3 a 4 ( κ 2 r 0 2 + 4 ) + κ [ 48 a 4 ( a r 0 ) 2 κ 2 r 0 ( 8 a 6 4 r 0 3 a 3 + r 0 5 a r 0 6 )
κ 4 r 0 6 a 2 ( a + r 0 ) ] cosh ( κ a κ r 0 ) + [ 48 a 4 + 2 κ 2 ( 8 a 6 24 a 5 r 0 4 a 3 r 0 3 r 0 6 )
+ κ 4 r 0 6 a ( a + 2 r 0 ) + κ 6 r 0 7 a 3 ] sinh ( κ a κ r 0 ) } ,
K 2 = ( 2 r 0 3 + b 3 + 2 ν b 3 ) [ 2 κ 3 a ( b 3 r 0 3 ) e κ ( b + a )
( 3 a + κ 2 a 3 κ 2 r 0 3 ) g 2 + ( a ) + 3 κ a 2 g 2 ( a ) ] ,  
E ( x ) = 1 t 1 e x t d t ,
and the functions g 0 ( r ) , g 1 ( r ) , and g 2 ± ( r ) are defined by Equations (7)–(9). As Equation (21) for μ is not affected by the fluid flow, Equation (47) for the effective conductivity does not depend on the reciprocal permeation length λ and the hydrodynamic condition at the outer border of the cell (i.e., applicable to both Happel and Kuwabara models). For a suspension of charged porous spheres ( r 0 = 0 and ν = 1 ), Equation (47) reduces to Equation (A3) in Appendix A.

3. Results and Discussion

3.1. Equilibrium Electrostatic Potential

The equilibrium electrostatic potential distribution ψ eq ( r ) inside a unit cell for a suspension of charged soft spheres in a salt-free solution is given by Equations (5) and (6). In Figure 2a,b, the dimensionless electric potential Z e ψ eq / k T is plotted versus the scaled coordinate r / a around a porous sphere ( r 0 = 0 ) and a soft sphere with r 0 / a = 1 / 2 , respectively, for various values of the scaled fixed charge density Q ¯ of the particles in the case of the particle volume fraction φ = ( a / b ) 3 = 0.125 . For specified values of Q ¯ and φ , as expected, the value of Z e ψ eq / k T is positive near the center of the porous particle and the hard-core surface of the soft particle, decreases monotonically with an increase in r / a (especially in the vicinity of the particle surface), and becomes negative after a critical value of r / a , which corresponds to n av for the local counterionic concentration and is a weak decreasing function of φ and Q ¯ . The values of ψ eq ( r 0 ) , ψ eq ( b ) , ψ eq ( r 0 ) ψ eq ( r ) , and ψ eq ( r ) ψ eq ( b ) (where r 0 < r < b ) are greater for a soft sphere with a smaller value of the core-to-particle radius ratio r 0 / a (i.e., greatest for a porous sphere) and increase with an increase in Q ¯ from zero at Q ¯ = 0 , keeping the other parameters unchanged.
The values of the dimensionless relative surface potential Z e ψ s / k T for a charged porous sphere ( r 0 = 0 ) and a charged soft sphere with r 0 / a = 4 / 5 in their corresponding salt-free suspensions calculated from Equation (10) are plotted versus the scaled fixed charge density Q ¯ in Figure 3a,b, respectively, for various values of φ . As expected, ψ s (or ψ eq ( r 0 ) or ψ eq ( b ) ) for a constant value of Q ¯ decreases with an increase in φ (i.e., a decrease in the relative cell size) and the effect of interactions among the particles in the suspensions on ψ s or ψ eq is significant. Note that, if Q ¯ is relatively high, our results obtained with the Debye–Hückel approximation may not apply for the salt-free suspensions with a low φ , since in this case the electric potential becomes high because of the presence of very small amount of ions. When Q ¯ is small (say, Q ¯ 10 , which corresponds to Z e ψ s / k T < 0.95 for φ 0.1 ), Z e ψ s / k T is nearly proportional to Q ¯ . When Q ¯ is relatively high, however, the increase of Z e ψ s / k T with Q ¯ is somewhat suppressed due to the counterion condensation occurring near the fixed charges of the particles (which is obviously underestimated with the linearized Poisson–Boltzmann equation for obtaining Equations (5) and (6)). For fixed values of Q ¯ and φ , the value of Z e ψ s / k T for a soft sphere decreases with an increase in the core-to-particle radius ratio r 0 / a from r 0 / a = 0 (a porous sphere).

3.2. Electrophoretic Mobility

The scaled electrophoretic mobility η Z e U / ε k T E in a bounded salt-free suspension of charged soft spheres, calculated from Equation (45) for the Happel and Kuwabara models, is plotted versus the scaled fixed charge density Q ¯ , shielding parameter λ a , particle volume fraction φ , and core-to-particle radius ratio r 0 / a in Figure 4, Figure 5, Figure 6 and Figure 7, respectively. For the case of porous spheres ( r 0 = 0 ), the coefficient ν in Equation (45) equals unity, independent of the boundary condition (13) or (14) for the external electric potential applied at the border of the unit cell. Consistent with the previous results [22,25] without the relaxation effect, η Z e U / ε k T E increases substantially with an increase in Q ¯ , a decrease in λ a , a decrease in φ , and a decrease in r 0 / a , keeping the other parameters unchanged. Additionally, η Z e U / ε k T E is nearly proportional to Q ¯ as Q ¯ 10 and its increase with Q ¯ is somewhat suppressed due to the counterion condensation effect (which is underestimated again with the Debye–Hückel approximation) when the value of Q ¯ is relatively high. For fixed values of Q ¯ , λ a , r 0 / a , and φ , the Kuwabara model leads to a less value of η Z e U / ε k T E than the Happel model does, and the Neumann condition for the electrostatic potential at the outer boundary of the unit cell results in a greater value of η Z e U / ε k T E than the Dirichlet condition does, but both differences are unsubstantial.
For comparison with the case of a suspension containing added electrolytes, η U / a 2 Q E (or η κ 2 U / Q E or η λ 2 U / Q E ), the electrophoretic mobility of the salt-free suspension of soft spheres normalized with the fixed charge density Q in the porous surface layer can be plotted versus the parameters κ a , λ a , r 0 / a , and φ . These figures, in which the value of η U / a 2 Q E decreases monotonically with increases in κ a , λ a , r 0 / a , and φ , are the same as those for a salt-containing suspension under the Debye–Hückel approximation [12,13,30].

3.3. Effective Electric Conductivity

The normalized effective electric conductivity Λ / Λ of a salt-free suspension of charged soft spheres calculated using Equation (47) is plotted versus the particle volume fraction φ and core-to-particle radius ratio r 0 / a in Figure 8 and Figure 9, respectively. The effect of particle interactions on the effective conductivity is significant. For a specified value of φ , Λ / Λ is a monotonic increasing function of the scaled fixed charge density Q ¯ and almost proportional to Q ¯ as Q ¯ 10 . For a given value of Q ¯ , Λ / Λ first increases with an increase in φ , attains a maximum, and then decreases with a further increase in φ . This behavior is different from that for the effective conductivity of a salt-containing suspension of charged soft spheres, where Λ / Λ (with Λ being the electric conductivity of the ionic solution in the absence of particles) as a function of φ depends on the valence, diffusivity, and bulk concentration of each ionic species [12,13]. For fixed values of Q ¯ and φ , the value of Λ / Λ decreases with an increase in r 0 / a from r 0 / a = 0 (porous spheres). Note that Λ / Λ can be less than unity when Q ¯ is small, φ is large, and r 0 / a is large.

4. Summary

The electrophoresis and electric conduction of a salt-free suspension of charged soft spheres with arbitrary values of the particle volume fraction φ , core-to-particle radius ratio r 0 / a , and shielding parameter λ a are analyzed by using a unit cell model in this work. Each soft sphere reduces to an entirely porous sphere if its hard core vanishes. The linearized Poisson–Boltzmann equation is solved for the equilibrium electric potential profile in the fluid containing the counterions only around the particle in a unit cell, and the counterionic continuity equation and modified Stokes/Brinkman equations are solved for the ionic electrochemical potential energy and fluid velocity distributions, respectively. Closed-form expressions for the electrophoretic velocity of the suspended particles with a low value of the scaled fixed charge density Q ¯ and the effective electric conductivity of the suspension are obtained in Equations (45) and (47), respectively. The effect of particle interactions on these transport characteristics is significant. Same as the case in a salt-containing suspension under the Debye–Hückel approximation, the scaled electrophoretic mobility in a salt-free suspension is a monotonic increasing function of Q ¯ and decreases with increases in r 0 / a , λ a , and φ , keeping the other parameters unchanged. The normalized effective electric conductivity of the salt-free suspension also increases with an increase in Q ¯ and with a decrease in r 0 / a , but is not a monotonic function of φ .

Author Contributions

Conceptualization, H.J.K.; methodology, H.J.K. and W.C.L.; investigation, H.J.K. and W.C.L.; writing—original draft preparation, H.J.K. and W.C.L.; writing—review and editing, H.J.K.; supervision, H.J.K.; funding acquisition, H.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology, Taiwan (Republic of China) grant number MOST 109-2221-E-002-062.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

For a salt-free suspension of charged porous spheres ( r 0 = 0 and ν = 1 ), Equation (45) for the electrophoretic velocity reduces to
U = Q E η κ 2 λ 2 a 2 [ 3 ( λ a ) 2 α ( λ a ) + 2 B 2 φ 5 / 3 ] { ε κ 4 a 2 3 Q φ [ k T Z e + ψ eq ( b ) ] [ 3 B 4 φ 5 / 3
3 ( λ a ) 4 α ( λ a ) φ 1 / 3 + ( λ a ) 2 ( B 1 2 B 2 φ 2 ) ] + 3 ( λ a ) 4 α ( λ a ) J 2 ( b )
+ 3 [ 2 B 2 J 3 ( 0 ) B 4 J 3 ( b ) ] φ 5 / 3 + ( λ a ) 2 [ 6 J α ( 0 ) + 3 α ( λ a ) J 3 ( 0 ) B 1 J 3 ( b )
B 3 J 0 ( b ) ] + ( λ a ) 2 [ 2 B 2 J 2 ( b ) + B 3 J 5 ( b ) + 6 α ( λ a ) J 3 ( 0 ) 6 J α ( 0 ) ] φ 5 / 3 }
for the Happel model and
U = Q E 15 η κ 2 λ 2 a 2 α ( λ a ) { ε κ 4 a 2 3 Q φ [ k T Z e + ψ eq ( b ) ] [ 5 B 1 18 ( λ a ) 2 α ( λ a ) φ 1 / 3
+ 2 ( 5 B 3 B 2 φ ) φ ] + 3 ( λ a ) 2 α ( λ a ) [ 5 J 2 ( b ) + J 5 ( b ) φ ] + 15 α ( λ a ) J 3 ( 0 )
5 B 1 J 3 ( b ) + 2 B 2 J 0 ( b ) φ 5 B 3 [ J 0 ( b ) + J 3 ( b ) φ ]
+ 30 [ J α ( 0 ) J α ( 0 ) φ + α ( λ a ) J 3 ( 0 ) φ ] }
for the Kuwabara model, and Equation (47) for the effective electric conductivity becomes
Λ = Λ + Λ Q ¯   e κ a 4 κ 5 a 2 b 3 g 2 + ( 0 ) { 2 [ e 2 κ a ( κ 2 a 2 3 κ a + 3 )
( κ 2 a 2 + 3 κ a + 3 ) ] g 0 ( a ) + 3 [ κ a + 1 + e 2 κ a ( κ a 1 ) ]
× [ ( κ 2 a 2 + 3 ) g 2 + ( a ) 3 κ a g 2 ( a ) 2 κ 3 b 3 e κ ( a + b ) ] } ,
where
B 1 = 3 α ( λ a ) + 2 ( λ a ) 3   cos h ( λ a ) ,
B 2 = 15 α ( λ a ) + ( λ a ) 2 [ λ a   cos h ( λ a ) 6   sin h ( λ a ) ] ,
B 3 = 6 α ( λ a ) + ( λ a ) 2 [ λ a   cos h ( λ a ) 3   sin h ( λ a ) ] ,
B 4 = 30 α ( λ a ) + 2 ( λ a ) 2 [ 7 λ a   cos h ( λ a ) 12   sin h ( λ a ) ]
+ ( λ a ) 4 [ λ a   cos h ( λ a ) 5   sin h ( λ a ) ] ,
and the functions g 0 ( r ) and g 2 ± ( r ) are defined by Equations (7) and (9).

References

  1. Levine, S.; Neale, G.H. The Prediction of Electrokinetic Phenomena within Multiparticle Systems I. Electrophoresis and Electroosmosis. J. Colloid Interface Sci. 1974, 47, 520–529. [Google Scholar] [CrossRef]
  2. O’Brien, R.W. The Electric Conductivity of a Dilute Suspension of Charged Particles. J. Colloid Interface Sci. 1981, 81, 234–248. [Google Scholar] [CrossRef]
  3. Zharkikh, N.I.; Shilov, V.N. Theory of Collective Electrophoresis of Spherical Particles in the Henry Approximation. Colloid J. USSR (Eng. Transl.) 1982, 43, 865–870. [Google Scholar]
  4. Ohshima, H.; Healy, T.W.; White, L.R. Approximate Analytic Expressions for the Electrophoretic Mobility of Spherical Colloidal Particles and the Conductivity of Their Dilute Suspensions. J. Chem. Soc. Faraday Trans. 2 1983, 79, 1613–1628. [Google Scholar] [CrossRef]
  5. Chen, S.B.; Keh, H.J. Axisymmetric Electrophoresis of Multiple Colloidal Spheres. J. Fluid Mech. 1992, 238, 251–276. [Google Scholar] [CrossRef]
  6. Liu, Y.C.; Keh, H.J. Electric Conductivity of a Dilute Suspension of Charged Composite Spheres. Langmuir 1998, 14, 1560–1574. [Google Scholar] [CrossRef]
  7. Ohshima, H. Electrical Conductivity of a Concentrated Suspension of Spherical Colloidal Particles. J. Colloid Interface Sci. 1999, 212, 443–448. [Google Scholar] [CrossRef]
  8. Borkovec, M.; Behrens, S.H.; Semmler, M. Observation of the Mobility Maximum Predicted by the Standard Electrokinetic Model for Highly Charged Amidine Latex Particles. Langmuir 2000, 16, 5209–5212. [Google Scholar] [CrossRef]
  9. Ding, J.M.; Keh, H.J. The Electrophoretic Mobility and Electric Conductivity of a Concentrated Suspension of Colloidal Spheres with Arbitrary Double-Layer Thickness. J. Colloid Interface Sci. 2001, 236, 180–193. [Google Scholar] [CrossRef] [Green Version]
  10. Carrique, F.; Cuquejo, J.; Arroyo, F.J.; Jimenez, M.L.; Delgado, Á.V. Influence of Cell-Model Boundary Conditions on the Conductivity and Electrophoretic Mobility of Concentrated Suspensions. Adv. Colloid Interface Sci. 2005, 118, 43–50. [Google Scholar] [CrossRef]
  11. Zholkovskij, E.K.; Masliyah, J.H.; Shilov, V.N.; Bhattacharjee, S. Electrokinetic Phenomena in Concentrated Disperse Systems: General Problem Formulation and Spherical Cell Approach. Adv. Colloid Interface Sci. 2007, 134–135, 279–321. [Google Scholar] [CrossRef] [PubMed]
  12. Keh, H.J.; Liu, C.P. Electric Conductivity and Electrophoretic Mobility in Suspensions of Charged Porous Spheres. J. Phys. Chem. C 2010, 114, 22044–22054. [Google Scholar] [CrossRef]
  13. Liu, H.C.; Keh, H.J. Electrophoresis and Electric Conduction in a Suspension of Charged Soft Spheres. Colloid Polym. Sci. 2016, 294, 1129–1141. [Google Scholar] [CrossRef]
  14. Galli, M.; Sáringer, S.; Szilágyi, I.; Trefalt, G. A Simple Method to Determine Critical Coagulation Concentration from Electrophoretic Mobility. Colloids Interfaces 2020, 4, 20. [Google Scholar] [CrossRef]
  15. Lin, W.C.; Keh, H.J. Diffusiophoresis in Suspensions of Charged Soft Particles. Colloids Interfaces 2020, 4, 30. [Google Scholar] [CrossRef]
  16. Lai, Y.C.; Keh, H.J. Transient Electrophoresis of a Charged Porous Particle. Electrophoresis 2020, 41, 259–265. [Google Scholar] [CrossRef]
  17. Delgado, Á.V.; Carrique, F.; Roa, R.; Ruiz-Reina, E. Recent Developments in Electrokinetics of Salt-Free Concentrated Suspensions. Curr. Opin. Colloid Interface Sci. 2016, 24, 32–43. [Google Scholar] [CrossRef]
  18. Ohshima, H. Ion Size Effect on Counterion Condensation around a Spherical Colloidal Particle in a Salt-Free Medium Containing Only Counterions. Colloid Polym. Sci. 2018, 296, 1293–1300. [Google Scholar] [CrossRef]
  19. Su, Y.W.; Keh, H.J. Electrokinetic Flow of Salt-Free Solutions in a Fibrous Porous Medium. J. Phys. Chem. B 2019, 123, 9724–9730. [Google Scholar] [CrossRef]
  20. Ohshima, H. Electrophoretic Mobility of a Spherical Colloidal Particle in a Salt-Free Medium. J. Colloid Interface Sci. 2002, 248, 499–503. [Google Scholar] [CrossRef]
  21. Ohshima, H. Electrokinetic Phenomena in a Dilute Suspension of Spherical Colloidal Particles in a Salt-Free Medium. Colloids Interfaces A 2003, 222, 207–211. [Google Scholar] [CrossRef]
  22. Ohshima, H. Electrophoretic Mobility of a Soft Particle in a Salt-Free Medium. J. Colloid Interface Sci. 2004, 269, 255–258. [Google Scholar] [CrossRef]
  23. Carrique, F.; Ruiz-Reina, E.; Arroyo, F.J.; Delgado, Á.V. Cell Model of the Direct Current Electrokinetics in Salt-Free Concentrated Suspensions: The Role of Boundary Conditions. J. Phys. Chem. B 2006, 110, 18313–18323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chiang, C.-P.; Lee, E.; He, Y.-Y.; Hsu, J.-P. Electrophoresis of a Spherical Dispersion of Polyelectrolytes in a Salt-Free Solution. J. Phys. Chem. B 2006, 110, 1490–1498. [Google Scholar] [CrossRef] [PubMed]
  25. He, Y.-Y.; Wu, E.; Lee, E. Electrophoresis in Suspensions of Charged Porous Spheres in Salt-Free Media. Chem. Eng. Sci. 2010, 65, 5507–5516. [Google Scholar] [CrossRef]
  26. Luo, R.H.; Keh, H.J. Electrophoresis and electric conduction in a salt-free suspension of charged particles. Electrophoresis 2021, 42. [Google Scholar] [CrossRef]
  27. Kong, C.Y.; Muthukumar, M. Monte Carlo study of adsorption of a polyelectrolyte onto charged surfaces. J. Chem. Phys. 1998, 109, 1522. [Google Scholar] [CrossRef]
  28. Hierrezuelo, J.; Szilagyi, I.; Vaccaro, A.; Borkovec, M. Probing Nanometer-Thick Polyelectrolyte Layers Adsorbed on Oppositely Charged Particles by Dynamic Light Scattering. Macromolecules 2010, 43, 9108–9116. [Google Scholar] [CrossRef]
  29. Lin, W.-C. Electrophoresis and Electric Conduction in Salt-Free Suspensions of Charged Soft Particles. Master’s Thesis, National Taiwan University, Taipei, Taiwan, 2021. [Google Scholar]
  30. Keh, H.J.; Chen, W.C. Sedimentation Velocity and Potential in Concentrated Suspensions of Charged Porous Spheres. J. Colloid Interface Sci. 2006, 296, 710–720. [Google Scholar] [CrossRef]
Figure 1. Geometrical sketch for the electrophoresis of a soft sphere in a unit cell for a suspension.
Figure 1. Geometrical sketch for the electrophoresis of a soft sphere in a unit cell for a suspension.
Colloids 05 00045 g001
Figure 2. Dimensionless equilibrium electric potential Z e ψ eq / k T for a soft sphere in a unit cell of φ = 0.125 ( a / b = 0.5 ) versus the scaled coordinate r / a with Q ¯ as a parameter: (a) r 0 / a = 0 (porous sphere); (b) r 0 / a = 1 / 2 .
Figure 2. Dimensionless equilibrium electric potential Z e ψ eq / k T for a soft sphere in a unit cell of φ = 0.125 ( a / b = 0.5 ) versus the scaled coordinate r / a with Q ¯ as a parameter: (a) r 0 / a = 0 (porous sphere); (b) r 0 / a = 1 / 2 .
Colloids 05 00045 g002
Figure 3. Dimensionless relative surface potential Z e ψ s / k T for a suspension of soft spheres versus Q ¯ with φ as a parameter: (a) r 0 / a = 0 (porous spheres); (b) r 0 / a = 4 / 5 .
Figure 3. Dimensionless relative surface potential Z e ψ s / k T for a suspension of soft spheres versus Q ¯ with φ as a parameter: (a) r 0 / a = 0 (porous spheres); (b) r 0 / a = 4 / 5 .
Colloids 05 00045 g003
Figure 4. Scaled electrophoretic mobility η Z e U / ε k T E versus Q ¯ for a suspension of soft spheres with λ a = 10 : (a) r 0 / a = 0 (porous spheres) with the solid and dashed curves representing calculations for the Happel and Kuwabara models, respectively; (b) φ = 0.125 for the Happel model with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Figure 4. Scaled electrophoretic mobility η Z e U / ε k T E versus Q ¯ for a suspension of soft spheres with λ a = 10 : (a) r 0 / a = 0 (porous spheres) with the solid and dashed curves representing calculations for the Happel and Kuwabara models, respectively; (b) φ = 0.125 for the Happel model with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Colloids 05 00045 g004
Figure 5. Scaled electrophoretic mobility η Z e U / ε k T E versus λ a for a suspension of soft spheres with φ = 0.125 : (a) r 0 / a = 0 (porous spheres) with the solid and dashed curves representing calculations for the Happily and Kuwabara models, respectively; (b) r 0 / a = 1 / 2 for the Happel model with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Figure 5. Scaled electrophoretic mobility η Z e U / ε k T E versus λ a for a suspension of soft spheres with φ = 0.125 : (a) r 0 / a = 0 (porous spheres) with the solid and dashed curves representing calculations for the Happily and Kuwabara models, respectively; (b) r 0 / a = 1 / 2 for the Happel model with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Colloids 05 00045 g005
Figure 6. Scaled electrophoretic mobility η Z e U / ε k T E versus φ for a suspension of soft spheres with Q ¯ = 4 : (a) r 0 / a = 0 (porous spheres) with the solid and dashed curves representing calculations for the Happel and Kuwabara models, respectively; (b) r 0 / a = 1 / 2 for the Happel model with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Figure 6. Scaled electrophoretic mobility η Z e U / ε k T E versus φ for a suspension of soft spheres with Q ¯ = 4 : (a) r 0 / a = 0 (porous spheres) with the solid and dashed curves representing calculations for the Happel and Kuwabara models, respectively; (b) r 0 / a = 1 / 2 for the Happel model with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Colloids 05 00045 g006
Figure 7. Scaled electrophoretic mobility η Z e U / ε k T E versus r 0 / a for a suspension of soft spheres with Q ¯ = 4 , λ a = 10 , and φ as a parameter for the Happel model. The solid and dashed curves denote results obtained from conditions (13) and (14), respectively, at the cell border.
Figure 7. Scaled electrophoretic mobility η Z e U / ε k T E versus r 0 / a for a suspension of soft spheres with Q ¯ = 4 , λ a = 10 , and φ as a parameter for the Happel model. The solid and dashed curves denote results obtained from conditions (13) and (14), respectively, at the cell border.
Colloids 05 00045 g007
Figure 8. Normalized effective electric conductivity Λ / Λ versus φ with Q ¯ as a parameter for a suspension of soft spheres: (a) r 0 / a = 0 (porous spheres); (b) r 0 / a = 1 / 2 with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Figure 8. Normalized effective electric conductivity Λ / Λ versus φ with Q ¯ as a parameter for a suspension of soft spheres: (a) r 0 / a = 0 (porous spheres); (b) r 0 / a = 1 / 2 with the solid and dashed curves denoting results obtained from conditions (13) and (14), respectively, at the cell border.
Colloids 05 00045 g008
Figure 9. Normalized effective electric conductivity Λ / Λ versus r 0 / a for a suspension of soft spheres with Q ¯ = 4 and φ as a parameter. The solid and dashed curves denote results obtained from conditions (13) and (14), respectively, at the cell border.
Figure 9. Normalized effective electric conductivity Λ / Λ versus r 0 / a for a suspension of soft spheres with Q ¯ = 4 and φ as a parameter. The solid and dashed curves denote results obtained from conditions (13) and (14), respectively, at the cell border.
Colloids 05 00045 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, W.C.; Keh, H.J. Electrophoretic Mobility and Electric Conductivity of Salt-Free Suspensions of Charged Soft Particles. Colloids Interfaces 2021, 5, 45. https://doi.org/10.3390/colloids5040045

AMA Style

Lin WC, Keh HJ. Electrophoretic Mobility and Electric Conductivity of Salt-Free Suspensions of Charged Soft Particles. Colloids and Interfaces. 2021; 5(4):45. https://doi.org/10.3390/colloids5040045

Chicago/Turabian Style

Lin, Wei C., and Huan J. Keh. 2021. "Electrophoretic Mobility and Electric Conductivity of Salt-Free Suspensions of Charged Soft Particles" Colloids and Interfaces 5, no. 4: 45. https://doi.org/10.3390/colloids5040045

APA Style

Lin, W. C., & Keh, H. J. (2021). Electrophoretic Mobility and Electric Conductivity of Salt-Free Suspensions of Charged Soft Particles. Colloids and Interfaces, 5(4), 45. https://doi.org/10.3390/colloids5040045

Article Metrics

Back to TopTop