Macroion correlation effects in electrostatic screening and thermodynamics of highly charged colloids

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.

Long-range forces and charge inversions in model charged colloidal dispersions at finite concentration

In charged colloidal dispersion systems the interest is in finding their stability conditions, phase transitions, and transport properties, either in bulk or confinement, among other physicochemical quantities, for which the knowledge of the dispersions' molecular structure and the associated macroion-macroion forces is crucial. To investigate these phenomena simple models have been proposed. Most of the theoretical and simulation studies on charged particles suspensions are at infinite dilution conditions. Hence, these studies have been focused on the electrolyte structure around one or two isolated central particle(s), where phenomena as charge reversal, charge inversion and surface charge amplification have been shown to be relevant. However, experimental studies at finite volume fraction exhibit interesting phenomenology which imply very long-range correlations. A simple, yet useful, model is the Colloidal Primitive Model, in which the colloidal dispersion is modeled as a mixture of size (and charge) asymmetrical hard spheres, at finite volume fraction. In this paper we review recent integral equations solutions for this model, where very long-range attractive-repulsive forces, as well as new long-range, giant charge inversions are reported. The calculated macroions radial distribution functions, charge distributions, and macroion-macroion forces are qualitatively consistent with existing experimental results, and Monte Carlo and molecular dynamics simulations.

Outsized Amplitude-Modulated Structure of Very-Long-Range Charge Inversions in Model Colloidal Dispersions

Most theoretical and simulation studies on charged suspensions are at infinite dilution and are focused on the electrolyte structure around one or two isolated particles. Some classic experimental studies with latex particle solutions exhibit interesting phenomenology which imply very-long-range correlations. Here, we apply an integral equation theory to a model charged macroion suspension, at finite volume fraction, and find an amplitude-modulated charge inversion structure, with outsized amplitudes and of very-long-range extension. These inversions are different from the standard charge inversions in that they occur at finite macroions' volume fraction, far away from the central macroion, are outsized, and increase, not decrease, with increasing particle charge and distance to the central particle, which is indicative of long-range correlations. We find our results to be in agreement with our Monte Carlo simulations and qualitatively consistent with existing experimental results.

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.

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.

Ultrafiltration of charge-stabilized dispersions at low salinity

We present a comprehensive study of cross-flow ultrafiltration (UF) of charge-stabilized suspensions, under low-salinity conditions of electrostatically strongly repelling colloidal particles. The axially varying permeate flux, near-membrane concentration-polarization (CP) layer and osmotic pressure profiles are calculated using a macroscopic diffusion-advection boundary layer method, and are compared with filtration experiments on aqueous suspensions of charge-stabilized silica particles. The theoretical description based on the one-component macroion fluid model (OCM) accounts for the strong influence of surface-released counterions on the renormalized colloid charge and suspension osmotic compressibility, and for the influence of the colloidal hydrodynamic interactions and electric double layer repulsion on the concentration-dependent suspension viscosity η, and collective diffusion coefficient Dc. A strong electro-hydrodynamic enhancement of Dc and η, and likewise of the osmotic pressure, is predicted theoretically, as compared with their values for a hard-sphere suspension. We also point to the failure of generalized Stokes-Einstein relations describing reciprocal relations between Dc and η. According to our filtration model, Dc is of dominant influence, giving rise to an only weakly developed CP layer having practically no effect on the permeate flux. This prediction is quantitatively confirmed by our UF measurements of the permeate flux using an aqueous suspension of charged silica spheres as the feed system. The experimentally detected fouling for the largest considered transmembrane pressure values is shown not to be due to filter cake formation by crystallization or vitrification.

Effective charges and virial pressure of concentrated macroion solutions

The stability of colloidal suspensions is crucial in a wide variety of processes, including the fabrication of photonic materials and scaffolds for biological assemblies. The ionic strength of the electrolyte that suspends charged colloids is widely used to control the physical properties of colloidal suspensions. The extensively used two-body Derjaguin-Landau-Verwey-Overbeek (DLVO) approach allows for a quantitative analysis of the effective electrostatic forces between colloidal particles. DLVO relates the ionic double layers, which enclose the particles, to their effective electrostatic repulsion. Nevertheless, the double layer is distorted at high macroion volume fractions. Therefore, DLVO cannot describe the many-body effects that arise in concentrated suspensions. We show that this problem can be largely resolved by identifying effective point charges for the macroions using cell theory. This extrapolated point charge (EPC) method assigns effective point charges in a consistent way, taking into account the excluded volume of highly charged macroions at any concentration, and thereby naturally accounting for high volume fractions in both salt-free and added-salt conditions. We provide an analytical expression for the effective pair potential and validate the EPC method by comparing molecular dynamics simulations of macroions and monovalent microions that interact via Coulombic potentials to simulations of macroions interacting via the derived EPC effective potential. The simulations reproduce the macroion-macroion spatial correlation and the virial pressure obtained with the EPC model. Our findings provide a route to relate the physical properties such as pressure in systems of screened Coulomb particles to experimental measurements.

Rotational self-diffusion in suspensions of charged particles: simulations and revised Beenakker-Mazur and pairwise additivity methods

We present a comprehensive joint theory-simulation study of rotational self-diffusion in suspensions of charged particles whose interactions are modeled by the generic hard-sphere plus repulsive Yukawa (HSY) pair potential. Elaborate, high-precision simulation results for the short-time rotational self-diffusion coefficient, D(r), are discussed covering a broad range of fluid-phase state points in the HSY model phase diagram. The salient trends in the behavior of D(r) as a function of reduced potential strength and range, and particle concentration, are systematically explored and physically explained. The simulation results are further used to assess the performance of two semi-analytic theoretical methods for calculating D(r). The first theoretical method is a revised version of the classical Beenakker-Mazur method (BM) adapted to rotational diffusion which includes a highly improved treatment of the salient many-particle hydrodynamic interactions. The second method is an easy-to-implement pairwise additivity (PA) method in which the hydrodynamic interactions are treated on a full two-body level with lubrication corrections included. The static pair correlation functions required as the only input to both theoretical methods are calculated using the accurate Rogers-Young integral equation scheme. While the revised BM method reproduces the general trends of the simulation results, it significantly underestimates D(r). In contrast, the PA method agrees well with the simulation results for D(r) even for intermediately concentrated systems. A simple improvement of the PA method is presented which is applicable for large concentrations.

Anomalous columnar order of charged colloidal platelets

Monte Carlo computer simulations are carried out for a model system of like-charged colloidal platelets in the isothermal-isobaric ensemble (NpT). The aim is to elucidate the role of electrostatic interactions on the structure of synthetic clay systems at high particle densities. Short-range repulsions between particles are described by a suitable hard-core model representing a discotic particle. This potential is supplemented with an electrostatic potential based on a Yukawa model for the screened Coulombic potential between infinitely thin disklike macro-ions. The particle aspect-ratio and electrostatic parameters were chosen to mimic an aqueous dispersion of thin, like-charged, rigid colloidal platelets at finite salt concentration. An examination of the fluid phase diagram reveals a marked shift in the isotropic-nematic transition compared to the hard cut-sphere reference system. Several statistical functions, such as the pair correlation function for the center-of-mass coordinates and structure factor, are obtained to characterize the structural organization of the platelets phases. At low salinity and high osmotic pressure we observe anomalous hexagonal columnar structures characterized by interpenetrating columns with a typical intercolumnar distance corresponding to about half of that of a regular columnar phase. Increasing the ionic strength leads to the formation of glassy, disordered structures consisting of compact clusters of platelets stacked into finite-sized columns. These so-called "nematic columnar" structures have been recently observed in systems of charge-stabilized gibbsite platelets. Our findings are corroborated by an analysis of the static structure factor from a simple density functional theory.

Coarse-graining intermolecular interactions in dispersions of highly charged colloids

Effective pair potentials between charged colloids, obtained from Monte Carlo simulations of two single colloids in a closed cell at the primitive model level, are shown to reproduce accurately the structure of aqueous salt-free colloidal dispersions, as determined from full primitive model simulations by Linse et al. (Linse, P.; Lobaskin, V. Electrostatic Attraction and Phase Separation in Solutions of Like-Charged Colloidal Particles. Phys. Rev. Lett.1999, 83, 4208). Excellent agreement is obtained even when ion-ion correlations are important and is in principle not limited to spherical particles, providing a potential route to coarse-grained colloidal interactions in more complex systems.

On the calculation of the structure of charge-stabilized colloidal dispersions using density-dependent potentials

The structure of charge-stabilized colloidal dispersions has been studied through a one-component model using a Yukawa potential with density-dependent parameters examined with integral equation theory and Monte Carlo simulations. Partial thermodynamic consistency was guaranteed by considering the osmotic pressure of the dispersion from the approximate mean-field renormalized jellium and Poisson-Boltzmann cell models. The colloidal structures could be accurately described by the Ornstein-Zernike equation with the Rogers-Young closure by using the osmotic pressure from the renormalized jellium model. Although we explicitly show that the correct effective pair-potential obtained from the inverse Monte Carlo method deviates from the Yukawa shape, the osmotic pressure constraint allows us to have a good description of the colloidal structure without losing information on the system thermodynamics. Our findings are corroborated by primitive model simulations of salt-free colloidal dispersions.

The metastable states of foam films containing electrically charged micelles or particles: experiment and quantitative interpretation

The stepwise thinning (stratification) of liquid films containing electrically charged colloidal particles (in our case - surfactant micelles) is investigated. Most of the results are applicable also to films from nanoparticle suspensions. The aim is to achieve agreement between theory and experiment, and to better understand the physical reasons for this phenomenon. To test different theoretical approaches, we obtained experimental data for free foam films from micellar solutions of three ionic surfactants. The theoretical problem is reduced to the interpretation of the experimental concentration dependencies of the step height and of the final film thickness. The surface charges of films and micelles are calculated by means of the charge-regulation model, with a counterion-binding (Stern) constant determined from the fit of surface tension isotherms. The applicability of three models was tested: the Poisson-Boltzmann (PB) model; the jellium-approximation (JA), and the cell model (CM). The best agreement theory/experiment was obtained with the JA model without using any adjustable parameters. Two theoretical approaches are considered. First, in the energy approach the step height is identified with the effective diameter of the charged micelles, which represents an integral of the electrostatic-repulsion energy calculated by the JA model. Second, in the osmotic approach the step height is equal to the inverse cubic root of micelle number density in the bulk of solution. Both approaches are in good agreement with the experiment if the suspension of charged particles (micelles) represents a jellium, i.e. if the particle concentration is uniform despite the field of the electric double layers. The results lead to a convenient method for determining the aggregation number of ionic surfactant micelles from the experimental heights of the steps.

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.

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.

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.

Density dependent interactions and structure of charged colloidal dispersions in the weak screening regime

We determine the structure of charge-stabilized colloidal suspensions at low ionic strength over an extended range of particle volume fractions using a combination of light and small angle neutron scattering experiments. The variation of the structure factor with concentration is analyzed within a one-component model of a colloidal suspension. We show that the observed structural behavior corresponds to a nonmonotonic density dependence of the colloid effective charge and the mean interparticle interaction energy. Our findings are corroborated by similar observations from primitive model computer simulations of salt-free colloidal suspensions.