Jun-Sik
Sin
*,
Hak-Chol
Pak
,
Kwang-Il
Kim
,
Kuk-Chol
Ri
,
Dok-Yong
Ju
,
Nam-Hyok
Kim
and
Chung-Sik
Sin
College of Physics, Kim Il Sung University, Taesong District, Pyongyang, DPR Korea. E-mail: ryongnam24@yahoo.com
First published on 12th November 2015
The response of a suspension under a variety of static or alternating external fields strongly depends on the equilibrium electric double layer that surrounds the colloidal particles in the suspension. The theoretical models for salt-free suspensions can be improved by incorporating non-uniform size effects and orientational ordering of water dipoles neglected in previous mean-field approaches, which are based on the Poisson–Boltzmann approach. Our model including non-uniform size effects and orientational ordering of water dipoles seems to have quite a promising effect because the model can predict the phenomena like a heavy decrease in relative permittivity of the suspension and counterion stratification near highly charged colloidal particles. In this work we numerically obtain the electric potential, the counterion concentration and the relative permittivity around a charged particle in a concentrated salt-free suspension corrected by non-uniform size effects and orientational ordering of water dipoles. The results show the worth of such corrections for medium to high particle charges at every particle volume fraction. We conclude that non-uniform size effects and orientational ordering of water dipoles are necessary for the development of new theoretical models to study non-equilibrium properties in concentrated colloidal suspensions.
The unique behaviours of salt-free suspensions can be understood by using physical quantities like electrophoretic mobility, electrical conductivity and dielectric response, which are electrohydrodynamic properties of suspension in an electrolyte solution.8–12 The physical quantities can be obtained by solving electrohydrodynamic equations with initial and boundary conditions. All of static or dynamic quantities are closely related to the properties of the equilibrium electric double layer surrounding the particles. If the salt-free suspension is dilute, for typical cases, particle–particle electrohydrodynamic interactions can be negligible. However, in the high regime of particle volume fraction, the response to external fields is affected by the interactions which give rise to the mathematical difficulty related with many-body interactions as well as the numerical problems for solving the equations without any approximation.
To overcome this problem, the authors of ref. 13–16 introduced a cell model approach, and the approach successfully was applied to dealing with electrokinetic and rheological properties in concentrated suspension. It was confirmed in ref. 17 that up to moderately strong electrostatic couplings, the cell model accurately predicts osmotic pressures of deionised suspensions in agreement with Monte Carlo simulations and renormalized-effective interaction approaches. In spite of the help of the cell model, because the theoretical models are based on the Poisson–Boltzmann (PB) approach, the models maintain the shortcomings of the Poisson–Boltzmann approach.
In ref. 18–22, the authors investigated dilute colloidal suspensions by classical PB equation, which neglects ionic correlation, ionic size effects and solvent polarization. They confirmed the fact that considering the distance of closest approach to the charged particle is important for explaining the overcharging mechanism near the charged particle without explicit consideration of ionic correlation.
Roa and coworkers23–26 considered the finite sizes of ions by extending Borukhov's approach27 to concentrated salt-free suspension with the help of the cell model. The consideration of finite ion sizes yielded reasonable counterion concentration profiles and electric potential profiles and allowed one to understand reasonably electrophoretic mobility, electric conductivity and dielectric response.
To well understand the properties of the electric double layer, different computational approaches such as Monte Carlo and numerical solutions of integral equations were introduced, but they are unable to make feasible predictions out of equilibrium and involve complex calculations.28,29
In the case of electrolytes, electric double layer theory has a lengthy history compared to salt-free suspension, dating back to Helmholtz.32 The original PB approach proposed by Gouy and Chapmann does not consider the finite volumes of ions, and neglects ionic correlations in the electrolyte.30,31 In fact, for typical situations including low to medium particle charges and monovalent ions in the electrolyte solution, this approach seems to be quite reasonable for representation of the equilibrium problem. The limitations of the PB approach yield unphysical ionic concentration profiles near highly charged interfaces and also inability to predict the overcharging phenomena.
To eliminate such a shortcoming of the PB approach, various modifications of the PB approach were pioneered and successfully applied to the practical problems. In order to account for the ionic size effect, Stern33 modified the PB approach considering the finite size effect of ions by combining the Helmholtz model32 with the Gouy–Chapmann mode l.31 Bikerman34 empirically extended Boltzmann distribution by correcting ion concentration for the volume excluded by ions. In the last two decades, researchers considered finite volumes of ions and water molecules within the lattice statistical mechanics approach.35–40 However, most of them are based on the assumption that different species of ions in electrolytes have equal size. Recent studies41–48 indicate that difference in sizes of ions enables one to understand phenomena like the asymmetric differential electric capacitance and the stratification of counterions.
It is well known that near a charged surface in an electrolyte the relative permittivity of the electrolyte solution varies according to the distance from the surface.49,50 Although the Booth model51,52 is widely used in many practices, it has the drawback that the model does not take into account the sizes of both ions and water molecules in the electrolyte solution. In ref. 53–55 the authors developed a water polarization model including size effects of ions and water molecules. The permittivity model of electrolytes well represents the fact that the permittivity of an electrolyte solution may be strongly decreased by orientational ordering of water dipoles and depletion of water molecules.
Recently, we ref. 56 have incorporated not only the non-uniform size effects of ions and water molecules in the electrolyte but also the orientational ordering of water dipoles into the Poisson–Boltzmann approach. It is evident that in the case of salt-free suspension as well as electrolytes, such effects will influence the equilibrium and non-equilibrium properties of the electric double layer.
Summarizing, in this study, we will focus on the influence of orientational ordering of water dipoles and non-uniform size effects on the electrostatic properties of the electric double layer of a charged particle in a concentrated salt-free suspension. In other words, we extend the Poisson–Boltzmann approach developed for aqueous electrolytes by us56 to the case of salt-free suspension. As in ref. 20, 21 and 23, our model also incorporates an excluded region in contact with the particle of a hydrated radius size, which yields more realistic representation of the solid–liquid interface and also predicts results in good agreement with experimental electrokinetic data. Solving numerically new equations, we will analyze the equilibrium surface potential, relative permittivity of the suspension and counterion concentration profiles inside a cell upon changing particle volume fraction, particle surface charge density and size of the counterions. In order to show the importance of the orientational ordering of water dipoles, the results will be compared with the ones of ref. 23 that takes into account only finite ion size.
![]() | (1) |
![]() | (2) |
Consider a unit volume of the cell around the charged particle. The entropy density is the logarithm of the number of translational and orientational arrangements of non-interacting nc counterions and ρ(ωi)ΔΩi(i = 1⋯N) water molecules, where ΔΩi = 2πsin
ωiΔω is an element of a solid angle and Δω = π/N. The counterion and the water molecule occupy volumes of Vc and Vw, respectively. We assume that the volumes are independent of the ionic concentrations. Within a lattice statistics approach each particle in the suspension occupies more than one cell of a lattice as in ref. 48 and 56. Unlike in ref. 48, orientational ordering of water dipoles as well as translational arrangements of counterions is explicitly considered. We first place nc counterions of the volume Vc. Accounting for orientational ordering of water dipoles, we put in ρ(ωi) (i = 0, 1,…) water molecules of the volume Vw in the lattice. The number of arrangements W is written as
![]() | (3) |
![]() | (4) |
ns = ncvc + nwvw. | (5) |
![]() | (6) |
The Euler–Lagrange equations for the Lagrangian are obtained and solved with respect to the functions nc and ρ(ω). The variation of the Lagrangian with respect to nc yields an equation from which nc and ρ(ω) in the suspension can be obtained.
![]() | (7) |
![]() | (8) |
![]() | (9) |
![]() | (10) |
In the case when the counterions and water molecules have the same sizes, that is, when Vw = Vc, we can recover the result of ref. 54 in eqn (9) and (10). When we neglect orientational ordering of water dipoles, our approach is identical to that of ref. 48.
The Euler–Lagrange equation for ψ(r) yields the Poisson–Boltzmann equation
![]() | (11) |
![]() | (12) |
![]() | (13) |
Differentiation of eqn (9), (10) and (5) with respect to the distance from the charged surface provides linear algebraic equations in terms of dnc/dr, dnw/dr, dg/dx, where g ≡ (λ − λ0)/kBT:
![]() | (14) |
![]() | (15) |
![]() | (16) |
Solving the system of eqn (14)–(16) for dnc/dr, dnw/dr, dg/dr results in the following coupled differential equations:
![]() | (17) |
![]() | (18) |
The electrostatic potential, the counterion concentration and the number density of water molecules are obtained by solving eqn (11), (17) and (18), respectively.
![]() | (19) |
To completely specify the problem the electric potential should be forced to be continuous at the surface r = a + Rc, and also its first derivative, which is related to the continuity of the normal component of the electric displacement at that surface. The boundary conditions needed for solving the problem are
![]() | (20) |
![]() | ||
Fig. 1 Cell model including the distance of closest approach of the counterions to the particle surface. |
Fig. 2 shows the electric potential distribution, Fig. 2(a), and the counterion concentration profiles, Fig. 2(b), along the cell. Solid lines represent the predictions of the equation considering only finite ion sizes without a distance of closest approach to the particle surface (MPB). Dashed lines stand for the results of the equation accounting for finite ion sizes and a distance of closest approach to the particle surface (MPBL). Circles correspond to the outcome of the equation taking into account finite ion sizes and orientational ordering of water dipoles (DPB). Plus signs exhibit the equation accounting for finite ion sizes, orientational ordering of water dipoles and a distance of closest approach to the particle surface (DPBL). The particle surface charge density has been chosen to be equal to 0.4 C m−2 and the particle volume fraction is ϕ = 0.5 (very concentrated suspension), which implies a normalized cell size of b/a = 1.26. In fact, the authors of ref. 23 also used the same size of parameters, but did not consider orientational ordering of water dipoles.
Fig. 2(a) shows that for our cases, the electric potentials are higher than ones of ref. 23. This is attributed to the fact that orientational ordering of water dipoles requires an additional electric energy. In the case when we take into account the distance of closest approach of the counterions to the particle surface (DPBL model), there is an higher increase of the surface potential as in ref. 23. Although in the case of the DPBL model the surface electric potential reaches 2 V that appears unfeasible, in the case of the DPB model the potential is about 1 V. It means that for the case of DPBL, in the region of a + Rc ≤ r ≤ b the electric potential difference is about 1 V, that is below the decomposition potential of water. The rest of potential difference is applied to the region of closest approach to the particle surface. As a result, σ = 0.4 C m−2 is the maximum possible surface charge density for the stable presence of suspension. This is consistent with the fact that in ref. 23 and 40 calculations were performed within the region of σ ≤ 0.4 C m−2, which is called as the physiological range. This shows that our approach well represents practical situations. We can observe in Fig. 2(b) that accounting for orientational ordering of water dipoles causes non-monotonous counterion concentration profiles, i.e. the position of the maximum counterion concentration is not located at the position of the surface of the particle. This is due to the statistical mechanical competition effect between orientational ordering of water dipoles and counterion condensation near the particle surface.56,57 For the counterion concentration profiles, our approach mathematically predicts two peaks. The higher one of them is equal to one of the case considering only finite ion size,23nmaxc. The higher peak occurs at a larger particle surface charge density than corresponding one of the lower peak. The formation of the higher peak requires quite a large particle surface charge density which corresponds to the unphysical value above the decomposition potential of water.
Fig. 3 shows the permittivity of the suspension for considering (circles, DPBL) or not (solid lines, DPB) the excluded region in contact with the particle, Fig. 3(a), and for a particle volume fraction ϕ = 0.5 for different ion sizes, Fig. 3(b), along the cell when orientational ordering of water dipoles is taken into account. The parameters are the same as in Fig. 2. Fig. 3(a) shows that the permittivity of the suspension near the particle surface decreases towards the particle surface. This decrease of εr is attributed to the increased orientational ordering of water dipoles and increased depletion of water molecules due to the accumulated counterions near the particle surface. The spatial region with decreased permittivity extends to a larger distance when we take into account the excluded region in contact with the particle. As shown in Fig. 3(b), it is obvious that the larger the counterion size, the wider the spatial region with permittivity lower than one of pure water. Eqn (12) and (13) explain that a small number of water dipoles and increased orientational ordering results in decreased permittivity of the suspension, as discussed in Fig. 3(a). Fig. 3(a) and (b) well represent that close to the charge surface the orientation of water molecules and depletion of water molecules result in spatial variation of permittivity.54,58,59
![]() | ||
Fig. 3 (a) Relative permittivity along the cell, considering (circles, DPBL) or not (solid lines, DPB) the excluded region in contact with the particle in the case when orientational ordering of water dipoles is taken into account. (b) Relative permittivity along the cell for different ion sizes, considering the excluded region in contact with the particle (DPBL). Triangles, circles and solid lines stand for nmaxc = 22 M, 10 M and 4 M. Other parameters are the same as in Fig. 2. |
Fig. 4 shows the effective charge divided by the particle charge along the cell for a particle volume fraction ϕ = 0.5, for different ion sizes and for MPBL and DPBL. The parameters are the same as in Fig. 2. It is shown in Fig. 4 that the effective charges for the method (DPBL) studied by considering orientational ordering of water dipoles are larger than the corresponding ones of ref. 23 (MPBL). In other words, orientational ordering of water dipoles diminishes the screening of the suspension. This is due to the non-monotonous counterion concentration profiles as shown in Fig. 2(a).
![]() | ||
Fig. 4 Effective charge divided by the particle charge along the cell for different ion sizes, considering or not orientational ordering of water dipoles in the case when the excluded region in contact with the particle is taken into account. For the case when orientational ordering of water dipoles is not taken into account (MPBL), squares, dashed lines and solid lines stand for nmaxc = 4 M, 10 M and 22 M, respectively. For the case when orientational ordering of water dipoles is taken into account (MPBL), crosses, circles and triangles correspond to nmaxc = 4 M, 10 M and 22 M, respectively. Other parameters are the same as in Fig. 3(b). |
Fig. 5 displays the equilibrium electric surface potential as a function of particle surface charge density for different ion sizes and ϕ = 0.5, Fig. 5(a), and for different particle volume fractions and nmaxc = 22 M, Fig. 5(b). It is shown in Fig. 5(a) that although in all cases the surface electric potential increases with particle surface charge density, the shapes of the profiles clearly differ from each other. In the cases when orientational ordering of water dipoles is taken into account (DPB and DPBL), the surface potential with particle surface charge density rapidly increases compared to the case of no accounting for water dipoles (MPB and MPBL), respectively. As shown in ref. 23, the PB equation predicts that in the region of low surface charge density, the surface potential sharply increases with the surface charge density. This is in agreement with the predictions of our theory as well as MPB. However, in the region of medium and high charge density, the surface electric potential predicted by MPB increases more rapidly than the ones of PB as proved in ref. 23. In the cases of DPB and DPBL, the behavior of the classical condensation effect is severely limited, because in the high regime of the surface charge the additional counterions are located in farther positions from the particle surface than ones predicted by MPB and MPBL as shown in Fig. 2(b). Consequently, in the case when orientational ordering of water dipoles is explicitly considered, further increase of the surface potential is attributed to the larger extent of counterion condensation compared to the case of no accounting for orientational ordering of water dipoles. It can be intuitively explained that from the medium regime, the results of DPBL and DPB theories are significantly deviated from the corresponding ones of MPBL and MPB, respectively. In fact, a high surface charge density yields a wide spatial region with low permittivity compared to the case of the low surface charge density. Combining the above fact with the boundary condition and the definition of electric field strength
, we can easily verify that the results of DPBL and DPB theories are higher than the ones of MPBL and MPB, respectively. Fig. 5(b) shows that in the case of orientational ordering of water dipoles is taken into account (DPBL), the surface electric potential increases with the particle volume fraction as in the case of MPBL. When the particle concentration rises, the space for the counterions inside the cell gets small and, consequently the screening of the particle charge is enhanced, thus reducing the value of the surface potential. This fact will have significant consequences on the electrokinetic properties of colloidal particles in concentrated salt-free suspensions, as shown in ref. 20 and 21 for dilute suspensions in electrolyte solutions.
![]() | ||
Fig. 5 (a) Surface electric potential for different values of particle surface charge density for different ion sizes, considering or not orientational ordering of water dipoles and/or the excluded region in contact with the particle (MPB, MPBL, DPB, and DPBL). (b) Surface potential against the particle surface charge density for different particle volume fraction values. Other parameters are the same as in Fig. 2. |
Fig. 6(a) and (b) show the number density of water molecules and the relative electric permittivity of the suspension at the particle surface according to the particle surface charge density for different ion sizes, respectively. Fig. 6(a) exhibits that accounting for simultaneously orientational ordering of water dipoles and finite size effects yields a non-monotonous behavior of the number density of water molecules at the particle surface. It should be noted that the smaller counterion size involves the shallower valley of the number density of water molecules at the particle surface. Then we note that all the number densities of water molecules are nearly equal in the high regime of surface charge density. According to eqn (12) and (13), the number density of water molecules is proportional to the relative permittivity at any position. On the other hand, the electric permittivity decreases with the magnitude of the electric field strength at any position. Even though in the case of larger number density of water molecules, the electric field strength gets slightly smaller, the sequence of magnitude of the electric permittivities for all cases studied is in agreement with one of the magnitude of corresponding number densities of water molecules as shown in Fig. 6(b).
![]() | ||
Fig. 6 Surface number density of water molecules (a) and surface relative permittivity (b) as functions of particle surface charge density for different ion sizes, considering the excluded region in contact with the particle and orientational ordering of water dipoles. Circles, dashed lines and solid lines stand for nmaxc = 22 M, 10 M and 4 M. Other parameters are the same as in Fig. 2. |
We have used a cell model approach to deal with particle–particle interaction, and extended Poisson–Boltzmann equations to account for such non-uniform size effects and orientational ordering of water dipoles. The theoretical procedure has followed that used by the authors56 but with the additional inclusion of a region of closest approach for counterions to the particle surface. Unlike in ref. 23, our PB approach can consider different sizes of different species of counterions as well as the difference in sizes of counterions and water molecules. We also have calculated the counterion concentration, equilibrium electric potential and relative permittivity for concentrated suspensions. The results have shown that orientational ordering of water dipoles is quite important for medium to high particle charges at every particle volume fraction, and even more if the distance of closest approach of the counterions to the particle surface is taken into account.
This equilibrium model presented in this study will be used as the base of non-equilibrium models for the response of a suspension to external electric fields. Experimental results concerning the static or dynamic electrophoretic mobilities, electrical conductivity and dielectric response should be compared with the predictions of non-equilibrium models. Although the stratification of counterions associated with the difference in size of the counterions has not been accessed in the present study, consideration of the difference in size of the counterions is very important for realistic salt-free suspensions which contain different species of counterions. Recently, the authors of ref. 44 and 45 studied the competition of multiple counterions of different valences and different sizes in binding to the surface of a spherical colloidal particle by both a mean-field theory and Monte Carlo simulations. Both the methods predicted the stratification of counterions around the charged surface and found that the ionic valence-to-volume ratios, instead of ionic valences alone, are the key factors that determine the binding of counterions to the charged surface. Because the authors of ref. 44 and 45 did not consider decreased permittivity of suspension near the charge particle, it will be interesting to study the phenomena by using our method. Additionally, the Monte Carlo approach accounting for both non-uniform size effects and orientational ordering of water dipoles is necessary for studying in-depth the phenomena. As a final conclusion, we emphasize that we can apply this model to predict the non-equilibrium properties in realistic salt-free concentrated suspensions. The task will be performed by the authors in the near future.
This journal is © the Owner Societies 2016 |