The renormalized Jellium model of colloidal suspensions with multivalent counterions
Abstract
An extension of the renormalized Jellium model which allows to study colloidal suspensions containing trivalent counterions is proposed. The theory is based on a modified PoissonBoltzmann equation which incorporates the effects of counterion correlations near the colloidal surfaces using a new boundary condition. The renormalized charges, the counterion density profiles, and osmotic pressures can be easily calculated using the modified renormalized Jellium model. The results are compared with the ones obtained using the traditional WignerSeitz (WS) cell approximation also with a new boundary condition. We find that while the thermodynamic functions obtained within the renormalized Jellium model are in a good agreement with their WS counterpart, the effective charges predicted by the two theories can be significantly different.
I Introduction
Strongly acidic or basic groups are often used to stabilize colloidal suspensions against flocculation and precipitation. The resulting electrostatic repulsion prevents the likecharged particles from coming sufficiently near one another where a shortrange van der Waals attraction can result in an irreversible association. It is well known that the PoissonBoltzmann (PB) theory describes accurately the properties of such colloidal suspensions containing a symmetric 1:1 electrolyte (1); (2). For these systems the correlations between small, strongly hydrated ions are weak, making the meanfield PB equation quasiexact. This is no longer true when the counterions are multivalent. In this case, electrostatic correlations between the ions can no longer be neglected and the meanfield PB equation looses its validity. Although the PB equation predicts very similar behaviors for suspensions containing either monovalent or multivalent counterions (3), both experiments and simulations show that this is not correct. Many interesting phenomena inherent to colloidal suspensions with multivalent counterions, such as colloidal charge reversal (4); (5); (6); (7), likecharged colloidal attraction (8), or the reversal of the electrophoretic mobility (9) can not be captured by the simple mean field approach. The description of these phenomena requires the use of more sophisticated and numerically more demanding approaches, such as the Monte Carlo simulations (10), the integral equations (11), or the density functional theories (12).
To account for the electrostatic correlations in suspensions containing multivalent counterions, while preserving the simplicity of the PB theory, dos Santos et al. proposed a WignerSeitz (WS) model with a modified boundary condition at the colloidal surface (13); (14). The model is based on the Shklovskii’s idea (15) that the validity of the meanfield theory can be extended if the correlations between the condensed counterions are properly taken into account. In the context of the PB equation, this can be done using a modified boundary condition which must be satisfied by the meanfield potential at the colloidal surface. Using this new boundary condition, the PB equation can be solved inside the WignerSeitz cell to yield the ionic density profiles from which the effective colloidal charge can be calculated using the Alexander prescription (13); (16).
The WS cell model is based on the assumption that the liquid state structure of a colloidal suspension is not very different from a crystal one. This allows one to avoid the complicated manybody problem by considering only the electrostatic interaction between one colloidal particle and its counterions and coions. As a further approximation, the geometry of the WS cell is taken to match the colloidal one (e. g. spherical). The electroneutrality condition is imposed by requiring that the electric field vanishes at the cell boundary. Thus, within the cell model there no direct electrostatic interactions between the different colloidal particles. Nevertheless, the effective colloidal charge obtained by matching the linear and nonlinear solutions of the PB equation at the cell boundary is often used within the DerjaguinLandauVerweyOverbeek (DLVO) potential (17) to account for particleparticle interaction.
The cell model should be particularly appropriate for suspensions with large volume fractions of colloidal particles, while its validity for dilute systems is questionable. In this paper we propose a more liquidstatelike point of view, which is particularly appropriate for suspensions with low volume fraction of colloidal particles.
Inside a suspension each colloidal particle feels the field produced by other particles and their condensed counterions. If we neglect the correlations between the colloidal particles they, together with their condensed counterions, can be thought to provide a uniform background in which free microions move. The effective colloidal charge and the background charge must be determined selfconsistently, resulting in a Renormalized Jellium Model (RJM) (18). This theory is particularly appropriate for calculating the effective colloidal charges because within the RJM, unlike in the cell model, the interaction potential between the colloidal particles has precisely the DLVO form (18). The RJM has been successfully used to calculate the effective charges and the structural and thermodynamic properties of colloidal suspensions containing 1:1 electrolyte (17); (19); (20); (21). In this paper the RJM will be extended to account for the electrostatic correlations in suspensions with trivalent counterions.
The paper is organized as follows. In section II, the general aspects of the model will be outlined. In section III we will describe the RJM with a modified boundary condition. In section IV the results of the theory will be presented. Finally, the conclusions and discussion will be given in section V.
Ii The Model
We consider a system of colloidal particles of radius and (negative) charge (uniformly distributed on the surfaces), and counterions of radius and charge , ( is the charge of proton) inside an aqueous solution of volume . We adopt a primitive model (PM) description in which the solvent is treated as a uniform continuum of dielectric constant . The overall charge neutrality requires , where and are the particle numbers of counterions and colloids, respectively. The typical length scale that characterizes the system is the Bjerrum length, defined as , which is Å, in water at room temperature. In order to maintain the simplicity of the model, image charge effects which can become nontrivial for multivalent counterions (22), are neglected within the primitive model approach adopted here — colloidal particle has the same dielectric constant as the solvent.
Because of strong electrostatic interaction between the counterions and the colloidal particles, many of the counterions become condensed onto colloidal surface (1); (2). A ’complex’ composed of one colloidal particles with a layer of its condensed counterions can then be regarded as a single entity carrying an effective charge (1); (2). The charge neutrality condition then becomes , where , and is the number density of free unassociated counterions.
Traditionally colloidal suspensions have been modeled using a crystallike approximation of a single macroion with its counterions inside a WignerSeitz cell, the radius of which is determined by the volume fraction of colloidal particles. Although this picture is appropriate at large concentrations — when strong correlations between the charged particles lead to a crystallike ordering — it might not be appropriate for dilute suspension. To avoid the WS cell approximation we will instead use a liquidstate RJM.
Suppose we fix one colloidal particle at the origin. Far from this particle the counterion density profile will have a Boltzmannlike form, , where is the mean electrostatic potential. Note that refers only to free, uncondensed, counterions. The density of other colloidal particles and of their condensed counterions provides a uniform neutralizing background . The mean electrostatic potential satisfies the JelliumPoissonBoltzmann (JPB) equation:
(1) 
where is the reduced potential, and defines the inverse effective Debye screening length. The selfconsistency condition requires that the effective charge, calculated from the far field solution of this equation, and the background charge must have the same value, . We should note that the screening of the electrostatic potential in the far field is produced only by the free (uncondensed) ions. A similar behavior is implicit within in the cell model if one tries to define the effective colloidal charge, as is done within the Alexander prescription (16). However, this renormalization is less transparent within the cell model than within the Jellium formalism.
Although the above meanfield equation works very well for monovalent ions , it becomes a rather poor approximation when suspension contains multivalent counterions (). In such cases, the strong electrostatic correlations between the condensed counterions lead to significant deviations from the PB theory. The deviations can be so strong that they qualitative modify the behavior of suspensions containing multivalent counterions. In such suspensions, one finds that the counterion condensation can become so strong as to reverse the sign of the effective colloidal charge. Furthermore, addition of a multivalent electrolyte can result in attraction between likecharged colloidal particles, thus destabilizing suspension against flocculation.
To include the effects of the counterion correlations one can proceed in a number of different ways. One approach is to use a weighteddensity functional theory (23) to account for the corrections to the meanfield electrostatic potential. Another approach is to use the integral equations theory (24). All of these methods, however, have their own drawbacks and are significantly more computationally demanding than the simple WS cell PB theory or the RJM.
We note that the counterion correlations are the strongest among the condensed counterions, since these ions are in the closest proximity of each other. Following Shklovskii we will, therefore, attempt to include the counterion correlations within the RJM using a modified boundary condition at the colloidal surface. The condensed counterions will be treated as a strongly correlated fluid — a concentrated quasitwo dimensional plasma. On the other hand, in the bulk the concentration of counterions is quite small, so that the correlations can be neglected and the mean field approach is still sufficient. By matching the two regimes, the JPB equation is recovered, but with a new boundary condition at the colloidal surface. From now on, we will restrict our attention to trivalent counterions, .
Iii The Theory
Strong electrostatic interactions between the colloidal particles and their counterions lead to counterion condensation. The condensed ions are in thermodynamic equilibrium with the free ions of suspension. Close to the colloidal surface the counterion chemical potential is,
(2) 
where is the de Broglie thermal wavelength and is the coursegrained density of condensed counterions. The correlational chemical potential is given by that of a twodimensional one component plasma (13); (25):
(3) 
were is the plasma parameter. Far away, in the bulk solution, the chemical potential is well approximated using the mean field electrostatic potential,
(4) 
The thermodynamic equilibrium between the condensed counterions and the free ions requires equality of Eqs. (2) and (4), from which follows
(5) 
This equation correctly describes the density profile after a short distance from the colloidal surface. This cutoff distance delimits the region where the microion correlations are important and the mean field approximation breaks down. However, since the range of counterion correlations is quite small and the PB density profile vary smoothly, Eq. (5) can be extrapolated all the way to the colloidal surface. This simplification then results in a new boundary condition for the JPB equation at the colloidal surface,
(6) 
The value of can be obtained using the strong coupling theory (26) and coarse graining procedure (13). We find
(7) 
The strong dependence of the coarse grained density on the colloidal radius is a direct consequence of the contact theorem which states that the difference between the contact and the bulk density of counterions is proportional to the square of the electric field (27). Since the bulk counterion density is much lower than the counterion concentration at the colloidal surface, the theorem requires that the density near the colloidal surface scale as (13) . Together, Eqs. (3), (6) and (7) provide a new relation between the density profile at the colloidal surface and the bare colloidal charge . This should be contrasted with the usual boundary condition for the PB equation, which, however, does not capture the strong counterion condensation resulting from electrostatic correlations at short distance from the colloidal surface.
The calculation of the electrostatic potential from Eq. (1) is now quite straight forward. Far from the colloidal particle the electrostatic potential has exactly the DLVO form,
(8) 
Now suppose that we know , , , and the volume fraction . Eq. (8) then provides the electrostatic potential and the electric field far from the colloidal surface. Using these as the initial conditions, we can numerically integrate the JPB equation (1) to obtain the density . Numerical integration is performed using the usual RungeKutta algorithm. Eqs. (6) and (7) can then be used to calculate the bare colloidal charge. In practice, of course, we know the bare charge and would like to calculate . This can be easily done by coupling the JPB solver with a root finding subroutine, such as the Newton–Raphson method. For each there is a corresponding . The root finding subroutine allows us to efficiently search the values of to find the one that corresponds to the given .
Iv Results


Fig. 1 compares the effective charge as a function of the bare charge for several colloidal concentrations, using the modified boundary condition within the RJM (solid curves) and within the WS cell model (dashed curves) (13). In all the calculations, the counterion and the colloidal radii used are Å and Å, respectively.
Although both models predict similar qualitative behaviors, there is a significant quantitative difference between the two effective charges. The effective charges calculated using the RJM lie below the ones calculated in the WS cell approach. A similar result was found for the monovalent ions: apparently, the colloidcolloid correlations, implicit in the cell model, reduce the counterion condensation (28). The difference, however, becomes smaller as the colloidal concentration decreases, as is also observed in the case of 1:1 electrolyte.
Unlike the monovalent situation, the effective charge in the case of trivalent counterions is not a monotonically increasing function of the bare charge, as can be clearly seen in Fig.1. Instead, after reaching a maximum, the effective charge decreases with further increase of the bare charge. This general trend is in perfect agreement with the findings of the Monte Carlo simulations (13).
In Fig. 2, the counterion density profiles calculated using both Jellium and WS cell models for a fixed bare charge , are displayed for different colloidal concentrations. We conclude that the difference in the density profiles is even less pronounced when compared with the corresponding discrepancy in the effective charges calculated using the two models. It is important to remember that the density profiles obtained using the above theory are only valid after some distance from the colloidal surface. At short distances, the strongcoupling regime dominates over the meanfield. The present theory coarsegrains the whole nearfield region into a modified boundary condition for the JPB equation.
Besides the effective charges and the ionic distributions, the thermodynamic properties can also be easily calculated in the framework of the renormalized Jellium and the cell models. The osmotic pressure within the RJM is a function of the bulk counterion concentration (17); (29):
(9) 
Employing Eq. (9), we have calculated the osmotic pressure as a function of the volume fraction for two fixed bare colloidal charges, and , using both the WS cell model and the RJM. The results are shown in Fig. 3. Again, we see only a very small difference at moderate volume fractions. For small volume fractions, the osmotic pressures are identical within the two models.
V Summary and Conclusions
We have extended the range of applicability of the renormalized Jellium model to describe suspensions with multivalent counterions in the absence of added salt. The model uses Shklovskii’s idea to include the counterion correlations as a modified boundary condition for the JPB equation. A similar strategy has already proven to be successful for the modified WS PB model (13); (14). Comparing the predictions of the renormalized Jellium model and the WS model, we find a quantitative difference in the values of the effective charge of colloidal particles. Since the far field potential within the renormalized Jellium formalism has precisely the DLVO form, we expect that the effective charges calculated using this formalism should be more reliable for structural calculations. On the other hand the measurable thermodynamic quantities such as the osmotic pressure come out to be practically identical in the two models.
Besides providing an alternative to the cell model, the RJM can be further extended to take into account the colloidcolloid correlations. The homogeneous background charge distribution can be replaced by a nonhomogeneous one related to the colloid distribution function. This approach has been successfully implemented for suspensions with monovalent electrolyte (20); (21). Unfortunately, for the case of multivalent counterions, inclusion of colloidal correlations is not so straightforward. To have an accurate structure function one needs to know the interaction potential not only in the far field, but also in the near field. The DLVO effective interaction potential, however, is not valid at short distances where it becomes strongly modified by the counterion correlations (30); (31). The work in this direction is now in progress.
Vi Acknowledgments
T.E.C. would like to acknowledge useful conversations with Alexandre P. dos Santos. This work was partially supported by the CNPq, INCTFCx, and by the USAFOSR under the grant FA95500910283.
References
 Y. Levin, Rep. Prog. Phys. 65, 1577 (2002).
 L. Belloni, J. Phys.: Condens. Matter 12, R549 (2000).
 A. Diehl and Y. Levin, J. Phys.: Condens. Matter 17, S3309 (2005).
 G. N. Patey, J. Chem. Phys. 72, 5763 (1980).
 A. Diehl and Y. Levin, J. Chem. Phys. 125, 054902 (2006).
 S. Pianegonda, M. C. Barbosa, and Y. Levin, Europhys. Lett. 71, 831 (2005).
 A. Y. Grosberg, T. T. Nguyen and B. I. Shklovskii, Rev. Mod. Phys. 74, 329 (2002).
 M. Deserno, C. Holm, and S. May, Macromolecules 33, 199 (2000); I. Rouzina and V. Bloomfield, J. Chem. Phys. 100, 9977 (1996); W. M. Gelbart, R. F. Bruinsma, P. A. Pincus, and V. A. Parsegian, Phys. Today 53, 38 (2000); F. J. Solis and M. O. de la Cruz, ibid. 54, 71 (2001).
 A. FernándezNieves, A. FernándezBarbero, F. J. de las Nieves, and B. Vincent, J. Chem. Phys. 123, 054905 (2005); A. MartinMolina, M. QuesadaPerez, F. GalisteoGonzalez and R. HidalgoAlvarez, J. Phys.: Condens. Matter 15, S3475 (2003); M. QuesadaPérez, J. CallejasFernández and R. HidalgoÁlvarez, Adv. Colloid Interface Sci. 95, 295 (2002).
 R. Messina, J. Phys.: Condens. Matter 21, 113102 (2009); J. Z. Wu, D. Bratko, H. W. Blanch, and J. M. Prausnitz, J. Chem. Phys. 111, 7084 (1999); V. Lobaskin and K. Qamhieh, J. Phys. Chem. B 107, 8022 (2003).
 M. D. CarbajalTinoco and P. GonzálezMozuelos, J. Chem. Phys. 117, 2344 (2002); J. A. Anta and S. Lago, J. Chem. Phys. 116, 10514 (2002); J. A. Anta, F. Bresme and S. Lago, J. Phys.: Condens. Matter 15 S3491 (2003); D. Léger and D. Levesque, J. Chem. Phys. 123, 124910 (2005).
 M. J. Stevens and M. O. Robbins, Europhys. Lett. 12, 81 (1990).
 A. P. dos Santos, A. Diehl, and Y. Levin, J. Chem. Phys. 130, 124110 (2009).
 A. P. dos Santos, A. Diehl, and Y. Levin, J. Chem. Phys. 132, 104105 (2010).
 B. I. Shklovskii, Phys. Rev. E 60, 5802 (1999).
 S. Alexander, P. M. Chaikin, P. Grant, G. J. Morales and P. Pincus, J. Chem. Phys. 80, 5776 (1984).
 B. V. Derjaguin and L. Landau, Acta Physicochim. URSS 14, 633 (1941); E. J. Verwey and J. T. G. Overbeek, Theory of the Stability of Lyophobic Colloids (Elsevier, Amsterdam, 1948).
 E. Trizac and Y. Levin, Phys. Rev. E 69, 031403 (2004).
 S. Pianegonda, E. Trizac, and Y. Levin, J. Chem. Phys. 126, 014702 (2007).
 R. CastañedaPriego, L. F. RojasOchoa, V. Lobaskin and J. C. MixtecoSánchez, Phys. Rev. E 74, 051408 (2006).
 T. E. Colla, Y. Levin and E. Trizac, J. Chem. Phys. 131, 074115 (2009).
 R. Messina, J. Chem. Phys. 117, 11062 (2002); Z. Y. Wang and Y. Q. Ma, J. Phys. Chem. B 114, 13386 (2010)
 A. Diehl, M. Tamashiro, M. C. Barbosa and Y. Levin, Physica A 274, 433 (1999).
 L. Belloni, J. Phys.: Condens. Matter 14, 9323 (2002); G. I. GuerreroGarcía, E. GonzálezTovar, M. LozadaCassou and F. de J. GuevaraRodríguez, J. Chem. Phys 123, 034703 (2005); A. MartínMolina, M. QuesadaPérez and R. HidalgoAlvarez, J. Phys. Chem. B 110, 1326 (2006).
 H. Totsuji, J. Phys. Soc. Jpn. 39, 253 (1975).
 A. G. Moreira and R. R. Netz, Europhys. Lett. 52, 705 (2000); A. Naji and R. R. Netz, Europhys. Lett. 57, 911 (2002); A. Naji and R. R. Netz, Eur. Phys. J. E 13, 43 (2004); A. Naji, S. Jungblut, A. G. Moreira and R. R. Netz, Physica A 352, 131 (2005); M. Kanduč, A. Naji and R. Podgornik, J. Chem. Phys. 132, 224703 (2010); A. Naji, M. Kanduc, R. R. Netz and R. Podgornik, arXiv.org:condmat/1008.0357v1 (2010).
 D. Henderson and L. Blum, J. Chem. Phys. 69, 5441 (1978); H. Wennerström, B. Jönsson and P. Linse, J. Chem. Phys. 76, 4665 (1982).
 J. Dobnikar, R. CastañedaPriego, H. H. von Grünberg and E. Trizac, New J. Phys. 8, 277 (2006).
 E. Trizac, J. Dobnikar, L. Belloni, H. H. von Grünberg and R. CastañedaPriego, Phys. Rev. E 75, 011401 (2007).
 P. GonzálezMozuelos and M. D. CarbajalTinoco, J. Chem. Phys. 109, 11074 (1998).
 V. Lobaskin, A. Lyubartsev and P. Linse, Phys. Rev. E 63, 020401 (2001).