Self-consistent effective interactions in charged colloidal suspensions

Evidence of electrostatic-enhanced depletion attraction in the structural properties and phase behavior of binary charged colloidal suspensions

In this paper we study the structure and phase behavior of binary mixtures of charged particles at low ionic strength. Due to the large size asymmetry between both species, light scattering measurements give us access only to the partial static structure factor that corresponds to the big particles. We observe that the addition of small charged colloids produces a decrease of the main peak of the measured static structure factor and a shift to larger scattering vector values. This finding is in agreement with theory based on integral equations with the Hypernetted-Chain Closure (HNC) relation. The effective interaction between two big particles due to the presence of small particles is obtained by a HNC inversion scheme and used in numerical simulations that adequately reproduce the experiments. We find that the presence of small particles induces an electrostatic depletion screening among the big colloids, creating around them an exclusion zone for the small charged colloids greater than that caused in the case of neutral small colloids, which in turn augments the depletion effect.

Effective electrostatic interactions in colloid-nanoparticle mixtures

Interparticle interactions and bulk properties of colloidal suspensions can be substantially modified by the addition of nanoparticles. Extreme asymmetries in size and charge between colloidal particles and nanoparticles present severe computational challenges to molecular-scale modeling of such complex systems. We present a statistical mechanical theory of effective electrostatic interactions that can greatly ease large-scale modeling of charged colloid-nanoparticle mixtures. By applying a sequential coarse-graining procedure, we show that a multicomponent mixture of charged colloids, nanoparticles, counterions, and coions can be mapped first onto a binary mixture of colloids and nanoparticles and then onto a one-component model of colloids alone. In a linear-response approximation, the one-component model is governed by a single effective pair potential and a one-body volume energy, whose parameters depend nontrivially on nanoparticle size, charge, and concentration. To test the theory, we perform molecular dynamics simulations of the two-component and one-component models and compute structural properties. For moderate electrostatic couplings, colloid-colloid radial distribution functions and static structure factors agree closely between the two models, validating the sequential coarse-graining approach. Nanoparticles of sufficient charge and concentration enhance screening of electrostatic interactions, weakening correlations between charged colloids and destabilizing suspensions, consistent with experiments.

Ion permeation inside microgel particles induced by specific interactions: from charge inversion to overcharging

In this work we have performed a theoretical study of a system formed by ionic microgels in the presence of monovalent salt with the help of Ornstein-Zernike integral equations within the hypernetted-chain (HNC) approximation. We focus in particular on analysing the role that the short-range specific interactions between the polymer fibres of the microgel and the incoming ions have on the equilibrium ion distribution inside and outside the microgel. For this purpose, a theoretical model based on the equilibrium partitioning effect is developed to determine the interaction between the microgel particle and a single ion. The results indicate that when counterions are specifically attracted to the polymer fibres of the microgel, an enhanced counterion accumulation occurs that induces the charge inversion of the microgel and a strong increase of the microgel net charge (or overcharging). In the case of coions, the specific attraction is also able to provoke the coion adsorption even though they are electrostatically repelled, and so increasing the microgel charge (true overcharging). Moreover, we show that ion adsorption onto the microgel particle is very different in swollen and shrunken states due to the competition between specific attraction and steric repulsion. In particular, ion adsorption occurs preferentially in the internal core of the particle for swollen states, whereas it is mainly concentrated in the external shell for de-swollen configurations. Finally, we observe the existence of a critical salt concentration, where the net charge of the microgels vanishes; above this inversion point the net charge of the microgels increases again, thus leading to reentrant stability of microgel suspensions.

Effective electrostatic interactions in mixtures of charged colloids

We present a theory of effective electrostatic interactions in polydisperse suspensions of charged macroions, generalizing to mixtures a theory previously developed for monodisperse suspensions. Combining linear response theory with a random phase approximation for microion correlations, we coarse grain the microion degrees of freedom to derive general expressions for effective macroion-macroion pair potentials and a one-body volume energy. For model mixtures of charged hard-sphere colloids, we give explicit analytical expressions. The resulting effective pair potentials have the same general form as predicted by linearized Poisson-Boltzmann theory, but consistently incorporate dependence on macroion density and excluded volume via the Debye screening constant. The volume energy, which depends on the average macroion density, contributes to the free energy and so can influence thermodynamic properties of deionized suspensions. To validate the theory, we compute radial distribution functions of binary mixtures of oppositely charged colloidal macroions from molecular dynamics simulations of the coarse-grained model (with implicit microions), taking effective pair potentials as input. Our results agree closely with corresponding results from more computationally intensive Monte Carlo simulations of the primitive model (with explicit microions). Simulations of a mixture with large size and charge asymmetries indicate that charged nanoparticles can enhance electrostatic screening of charged colloids. The theory presented here lays a foundation for future large-scale modeling of complex mixtures of charged colloids, nanoparticles, and polyelectrolytes.

Poisson-Boltzmann theory of charged colloids: limits of the cell model for salty suspensions

Thermodynamic properties of charge-stabilized colloidal suspensions and polyelectrolyte solutions are commonly modelled by implementing the mean-field Poisson-Boltzmann (PB) theory within a cell model. This approach models a bulk system by a single macroion, together with counterions and salt ions, confined to a symmetrically shaped, electroneutral cell. While easing numerical solution of the nonlinear PB equation, the cell model neglects microion-induced interactions and correlations between macroions, precluding modelling of macroion ordering phenomena. An alternative approach, which avoids the artificial constraints of cell geometry, exploits the mapping of a macroion-microion mixture onto a one-component model of pseudo-macroions governed by effective interparticle interactions. In practice, effective-interaction models are usually based on linear-screening approximations, which can accurately describe strong nonlinear screening only by incorporating an effective (renormalized) macroion charge. Combining charge renormalization and linearized PB theories, in both the cell model and an effective-interaction (cell-free) model, we compute osmotic pressures of highly charged colloids and monovalent microions, in Donnan equilibrium with a salt reservoir, over a range of concentrations. By comparing predictions with primitive model simulation data for salt-free suspensions, and with predictions from nonlinear PB theory for salty suspensions, we chart the limits of both the cell model and linear-screening approximations in modelling bulk thermodynamic properties. Up to moderately strong electrostatic couplings, the cell model proves accurate for predicting osmotic pressures of deionized (counterion-dominated) suspensions. With increasing salt concentration, however, the relative contribution of macroion interactions to the osmotic pressure grows, leading predictions from the cell and effective-interaction models to deviate. No evidence is found for a liquid-vapour phase instability driven by monovalent microions. These results may guide applications of PB theory to colloidal suspensions and other soft materials.

Electrical potentials of two identical planar, cylindrical, and spherical colloidal particles in a salt-free medium

A perturbation solution based on the curvature of the surface of a particle is derived for the electrical potential of two planar, cylindrical, and spherical colloidal particles immersed in a salt-free medium. Compared with the available results in the literature, the present approach improves both the accuracy and the applicable range in the level of the surface potential, the valence of counterions, and the thickness of the double-layer surrounding a particle. The zero-order perturbation solution is exact for the case of two planar particles. Under typical conditions, the averaged percentage deviation of the first-order perturbation solution is on the order of ca. 5% for both two cylindrical and two spherical particles.

Electroneutrality and phase behavior of colloidal suspensions

Several statistical mechanical theories predict that colloidal suspensions of highly charged macroions and monovalent microions can exhibit unusual thermodynamic phase behavior when strongly deionized. Density-functional, extended Debye-Hückel, and response theories, within mean-field and linearization approximations, predict a spinodal phase instability of charged colloids below a critical salt concentration. Poisson-Boltzmann cell model studies of suspensions in Donnan equilibrium with a salt reservoir demonstrate that effective interactions and osmotic pressures predicted by such theories can be sensitive to the choice of reference system, e.g., whether the microion density profiles are expanded about the average potential of the suspension or about the reservoir potential. By unifying Poisson-Boltzmann and response theories within a common perturbative framework, it is shown here that the choice of reference system is dictated by the constraint of global electroneutrality. On this basis, bulk suspensions are best modeled by density-dependent effective interactions derived from a closed reference system in which the counterions are confined to the same volume as the macroions. Lower-dimensional systems (e.g., monolayers, clusters), depending on the strength of macroion-counterion correlations, may be governed instead by density-independent effective interactions tied to an open reference system with counterions dispersed throughout the reservoir, possibly explaining the observed structural crossover in colloidal monolayers and anomalous metastability of colloidal crystallites.

Secondary Minimum Coagulation in Charged Colloidal Suspensions from Statistical Mechanics Methods

A statistical mechanics approach is applied to predict the critical parameters of coagulation in the secondary minimum for charged colloidal suspensions. This method is based on the solution of the reference hypernetted chain (RHNC) integral equation, and it is intended to estimate only the locus of the critical point instead of the full computation of the "gas-liquid" coexistence. We have used an extrapolation procedure due to the lack of solution of the integral equation in the vicinity of the critical point. Knowing that the osmotic isothermal compressibility of the colloidal system should ideally diverge in the critical point, we work out the critical salt concentration for which the inverse of the compressibility should be zero. This extrapolation procedure is more rapid than that previously proposed by Morales and co-workers [Morales, V.; Anta, J. A.; Lago, S. Langmuir 2003, 19, 475], and it is shown to give equivalent results. We also present experimental results about secondary minimum coagulation for polystyrene latexes and use our method to reproduce the experimental trends. The comparison between theory and experiment is quite good for all colloidal diameters studied.

Approximate Analytical Expressions for the Electrical Potential between Two Planar, Cylindrical, and Spherical Surfaces

Approximate analytical expressions for the electrical potential of planar, cylindrical, and spherical surfaces are derived for the case in which the dispersion medium contains counterions only. On the basis of the results for single surfaces, those for two identical surfaces can be derived. The curvature effect of a surface on the electrical potential distribution can be neglected when the order of its radius exceeds approximately 100 times the thickness of the corresponding double layer. If this effect needs to be considered, it can be taken into account by multiplying a correction function by the electrical potential of a planar surface. The electrical potential at the center between two derived surfaces is readily applicable to the evaluation of the electrostatic force per unit area between two surfaces, or the osmotic pressure. For the same set of parameters, the magnitudes of the osmotic pressure for various types of surfaces rank as follows: planar surface > cylindrical surfaces > spherical surfaces.

Entropic forces in dilute colloidal systems

Depletion forces can be accounted for by a contraction of the description in the framework of the integral equations theory of simple liquids. This approach includes, in a natural way, the effects of the concentration on the depletion forces, as well as energetic contributions. In this paper we systematically study this approach in a large variety of dilute colloidal systems composed of spherical and nonspherical hard particles, in two and in three dimensions, in the bulk and in front of a hard wall with a relief pattern. We show by this way the form in which concentration and geometry determine the entropic interaction between colloidal particles. The accuracy of our results is corroborated by comparison with computer simulations.

Effective interactions in the colloidal suspensions from hypernetted-chain theory

The hypernetted-chain (HNC) Ornstein-Zernike integral equations are used to determine the properties of simple models of colloidal solutions where the colloids and ions are immersed in a solvent considered as a dielectric continuum and have a size ratio equal to 80 and a charge ratio varying between 1 and 4000. At an infinite dilution of colloids, the effective interactions between colloids and ions are determined for ionic concentrations ranging from 0.001 to 0.1 mol/l and compared to those derived from the Poisson-Boltzmann theory. At finite concentrations, we discuss on the basis of the HNC results the possibility of an unambiguous definition of the effective interactions between the colloidal molecules.

Nonlinear screening and effective electrostatic interactions in charge-stabilized colloidal suspensions

A nonlinear response theory is developed and applied to electrostatic interactions between spherical macroions, screened by surrounding microions, in charge-stabilized colloidal suspensions. The theory describes leading-order nonlinear response of the microions (counterions, salt ions) to the electrostatic potential of the macroions and predicts microion-induced effective many-body interactions between macroions. A linear response approximation [A.R. Denton, Phys. Rev. E 62, 3855 (2000)] yields an effective pair potential of screened-Coulomb (Yukawa) form, as well as a one-body volume energy, which contributes to the free energy. Nonlinear response generates effective many-body interactions and essential corrections to both the effective pair potential and the volume energy. By adopting a random-phase approximation (RPA) for the response functions, and thus neglecting microion correlations, practical expressions are derived for the effective pair and triplet potentials and for the volume energy. Nonlinear screening is found to weaken repulsive pair interactions, induce attractive triplet interactions, and modify the volume energy. Numerical results for monovalent microions are in good agreement with available ab initio simulation data and demonstrate that nonlinear effects grow with increasing macroion charge and concentration and with decreasing salt concentration. In the dilute limit of zero macroion concentration, leading-order nonlinear corrections vanish. Finally, it is shown that nonlinear response theory, when combined with the RPA, is formally equivalent to the mean-field Poisson-Boltzmann theory and that the linear response approximation corresponds, within integral-equation theory, to a linearized hypernetted-chain closure.

Solution of the Percus-Yevick equation for square well spherocylinders

The Percus-Yevick equation for square-well spherocylinders has been numerically solved for some selected orientations following a methodology proposed previously for different fluids of elongated molecules. The equation is solved for particles of aspect ratios ranging from L/sigma=0.3 up to L/sigma=5.0, attractive range lambda/sigma=1.5, and packing fractions within eta=0.1-0.3. The resulting pair correlation functions are checked against isothermal-isobaric Monte Carlo simulations and good agreement is found for the short-range structure, at intermolecular distances within one molecular diameter sigma to contact for each of the selected orientations. At larger distances, the integral equation tends to overestimate the pair correlations. The results confirm the prediction of reference-system average Mayer-function perturbation theory for short aspect ratios, reaching the Onsager limit for the greater aspect ratios. Some instabilities of the solution for the longest models and higher densities are tentatively discussed in terms of their possible relation to frustration phenomena found in some polymer and complex systems.