The polydisperse cell model: Nonlinear screening and charge renormalization in colloidal mixtures

A self-consistent Ornstein-Zernike jellium for highly charged colloids (microgels) in suspensions with added salt

This work discusses a jellium scheme, built within the framework of the multicomponent Ornstein-Zernike (OZ) equation, which is capable of describing the collective structure of suspensions of highly charged colloids with added salt, even in the presence of finite-size multivalent microions. This approach uses a suitable approximation to decouple the microion-microion correlations from the macroion-microion profiles, which in combination with the methodology from the dressed ion theory (DIT) gives a full account of the electrostatic effective potential among the colloids. The main advantages of the present contribution reside in its ability to manage the short-range potentials and non-linear correlations among the microions, as well as its realistic characterization of the ionic clouds surrounding each macroion. The structure factors predicted by this jellium scheme are contrasted with previously reported experimental results for microgel suspensions with monovalent salts (2019Phys. Rev. E100032602), thus validating its high accuracy in these situations. The present theoretical analysis is then extended to microgel suspensions with multivalent salts, which reveals the prominent influence of the counterion valence on the makeup of the effective potentials. Although the induced differences may be difficult to identify through the mesoscopic structure, our results suggest that the microgel collapsing transition may be used to enhance these distinct effects, thus giving a feasible experimental probe for these phenomena.

Jellium and cell model for titratable colloids with continuous size distribution

A good understanding and determination of colloidal interactions is paramount to comprehend and model the thermodynamic and structural properties of colloidal suspensions. In concentrated aqueous suspensions of colloids with a titratable surface charge, this determination is, however, complicated by the density dependence of the effective pair potential due to both the many-body interactions and the charge regulation of the colloids. In addition, colloids generally present a size distribution which results in a virtually infinite combination of colloid pairs. In this paper, we develop two methods and describe the corresponding algorithms to solve this problem for arbitrary size distributions. An implementation in Nim is also provided. The methods, inspired by the seminal work of Torres et al., [J. Chem. Phys. 128, 154906 (2008)] are based on a generalization of the cell and renormalized jellium models to polydisperse suspensions of spherical colloids with a charge regulating boundary condition. The latter is described by the one-pK-Stern model. The predictions of the models are confronted to the equations of state of various commercially available silica dispersions. The renormalized Yukawa parameters (effective charges and screening lengths) are also calculated. The importance of size and charge polydispersity as well as the validity of these two models is discussed in light of the results.

Short- and long-time diffusion and dynamic scaling in suspensions of charged colloidal particles

We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in concentrated suspensions of charge-stabilized colloidal spheres. In theory and simulation, the spheres are assumed to interact directly by a hard-core plus screened Coulomb effective pair potential. The intermediate scattering function, f_c(q, t), is calculated by elaborate accelerated Stokesian dynamics (ASD) simulations for Brownian systems where many-particle hydrodynamic interactions (HIs) are fully accounted for, using a novel extrapolation scheme to a macroscopically large system size valid for all correlation times. The study spans the correlation time range from the colloidal short-time to the long-time regime. Additionally, Brownian Dynamics (BD) simulation and mode-coupling theory (MCT) results of f_c(q, t) are generated where HIs are neglected. Using these results, the influence of HIs on collective and self-diffusion and the accuracy of the MCT method are quantified. It is shown that HIs enhance collective and self-diffusion at intermediate and long times. At short times self-diffusion, and for wavenumbers outside the structure factor peak region also collective diffusion, are slowed down by HIs. MCT significantly overestimates the slowing influence of dynamic particle caging. The dynamic scattering functions obtained in the ASD simulations are in overall good agreement with our dynamic light scattering (DLS) results for a concentration series of charged silica spheres in an organic solvent mixture, in the experimental time window and wavenumber range. From the simulation data for the time derivative of the width function associated with f_c(q, t), there is indication of long-time exponential decay of f_c(q, t), for wavenumbers around the location of the static structure factor principal peak. The experimental scattering functions in the probed time range are consistent with a time-wavenumber factorization scaling behavior of f_c(q, t) that was first reported by Segrè and Pusey [Phys. Rev. Lett. 77, 771 (1996)] for suspensions of hard spheres. Our BD simulation and MCT results predict a significant violation of exact factorization scaling which, however, is approximately restored according to the ASD results when HIs are accounted for, consistent with the experimental findings for f_c(q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, f_s(q, t), and its non-Gaussian parameter α₂(t) and for the particle mean squared displacement W(t) and its time derivative. Since self-diffusion properties are not assessed in standard DLS measurements, a method to deduce W(t) approximately from f_c(q, t) is theoretically validated.

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.

Drying paint: from micro-scale dynamics to mechanical instabilities

Charged colloidal dispersions make up the basis of a broad range of industrial and commercial products, from paints to coatings and additives in cosmetics. During drying, an initially liquid dispersion of such particles is slowly concentrated into a solid, displaying a range of mechanical instabilities in response to highly variable internal pressures. Here we summarize the current appreciation of this process by pairing an advection-diffusion model of particle motion with a Poisson-Boltzmann cell model of inter-particle interactions, to predict the concentration gradients in a drying colloidal film. We then test these predictions with osmotic compression experiments on colloidal silica, and small-angle X-ray scattering experiments on silica dispersions drying in Hele-Shaw cells. Finally, we use the details of the microscopic physics at play in these dispersions to explore how two macroscopic mechanical instabilities-shear-banding and fracture-can be controlled.This article is part of the themed issue 'Patterning through instabilities in complex media: theory and applications.'

Effective pair potential between charged nanoparticles at high volume fractions

Simulations of charged colloidal dispersions are technically challenging. One possible workaround consists in reducing the system to the colloids only, whose interactions are described through an effective pair potential, w^f. Still, the determination of w^f is difficult mainly because it depends on the colloidal density, ϕ. Here we propose to calculate w^f from simulations of a pair of colloids placed in a cubic box with periodic boundary conditions. The variation in ϕ is mimicked by an appropriate change in the concentration of counterions neutralized by an homogeneous background charge. The method is tested at the level of the primitive model. A good description of the structure of the colloidal dispersion is obtained in the low and high coupling regimes, even at high ϕ (≈30%). Furthermore, the method can easily be used in popular molecular simulation program packages and extended to non-spherical objects.

Renormalized jellium model for colloidal mixtures

In an attempt to quantify the role of polydispersity in colloidal suspensions, we present an efficient implementation of the renormalized jellium model for a mixture of spherical charged colloids. The different species may have different size, charge, and density. Advantage is taken from the fact that the electric potential pertaining to a given species obeys a Poisson's equation that is species independent; only boundary conditions do change from one species to the next. All species are coupled through the renormalized background (jellium) density, that is determined self-consistently. The corresponding predictions are compared to the results of Monte Carlo simulations of binary mixtures, where Coulombic interactions are accounted for exactly, at the primitive model level (structureless solvent with fixed dielectric permittivity). An excellent agreement is found.

Alternating strings and clusters in suspensions of charged colloids

We report the formation of alternating strings and clusters in a binary suspension of repulsive charged colloids with double layers larger than the particle size. Within a binary cell model we include many-body and charge-regulation effects under the assumption of a constant surface potential, and consider their repercussions on the two-particle interaction potential. We find that the formation of induced dipoles close to a charge-reversed state may explain the formation of these structures. Finally, we will touch upon the formation of dumbbells and small clusters in a one-component system, where the effective electrostatic interaction is always repulsive.

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.

Renormalized jellium mean-field approximation for binary mixtures of charged colloids

In this work the renormalized jellium model of colloidal suspensions, originally proposed by Trizac and Levin [Phys. Rev. E 69, 031403 (2004)], is extended to study mechanisms of charge renormalization in binary mixtures of charged colloids. We here apply our recent reformulation that introduces the requirement of self-consistency directly into the Poisson-Boltzmann equation, i.e., the background charge is explicitly replaced by the effective one, thus facilitating the whole charge renormalization scheme. We briefly discuss the reformulated model for monodisperse charged suspensions composed of either spheres or rods. In particular, we put emphasis on the effects of the surface charge variation, mixture composition, and particle size on the charge regulation of charge-stabilized colloidal suspensions.

Nonadditivity in the effective interactions of binary charged colloidal suspensions

Based on primitive model computer simulations with explicit microions, we calculate the effective interactions in a binary mixture of charged colloids with species A and B for different size and charge ratios. An optimal pairwise interaction is obtained by fitting the many-body effective forces. This interaction is close to a Yukawa (or Derjaguin-Landau-Verwey-Overbeek (DLVO)) pair potential but the AB cross-interaction is different from the geometric mean of the two direct AA and BB interactions. As a function of charge asymmetry, the corresponding nonadditivity parameter is first positive, then significantly negative and is then positive again. We finally show that an inclusion of nonadditivity within an optimal effective Yukawa model gives better predictions for the fluid pair structure than DLVO theory.

Poisson-Boltzmann cell model for heterogeneously charged colloids

We introduce the Poisson-Boltzmann cell model for spherical colloidal particles with a heterogeneous surface charge distribution. This model is obtained by generalizing existing cell models for mixtures of homogeneously charged colloidal spheres. Our model has similar features as Onsager's second-virial theory for liquid crystals, but it predicts no orientational ordering if there is no positional ordering. This implies that all phases of heterogeneously charged colloids that are liquidlike with respect to translational degrees of freedom are also isotropic with respect to particle orientation.

Enhanced crystal stability in a binary mixture of charged colloidal spheres

We report on the phase behavior of a binary charged sphere mixture of size ratio Gamma=0.68 and charge ratio Lambda approximately 1 as a function of composition p and number density n . For p=0.1-0.3 we observe freezing at densities well below the freezing densities of the pure components. At all compositions our data indicate the formation of substitutional alloy crystals of body centered cubic structure. No indications for compound formation were observed. Rather, our findings point at the first observation of an upper azeotrope. Measurements of the crystallization kinetics reveal a combined density and composition dependence of growth velocities and nucleation rate densities, with small but significant anomalies at p*=0.2 . These correlated deviations can be rationalized within classical theories of solidification and suggest an increased similarity between melt and solid for this particular composition.