Phase separation in charge-stabilized colloidal suspensions: influence of nonlinear screening

Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion-microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations.

Concentration-dependent swelling and structure of ionic microgels: simulation and theory of a coarse-grained model

We study swelling and structural properties of ionic microgel suspensions within a comprehensive coarse-grained model that combines the polymeric and colloidal natures of microgels as permeable, compressible, charged spheres governed by effective interparticle interactions. The model synthesizes the Flory-Rehner theory of cross-linked polymer gels, the Hertz continuum theory of effective elastic interactions, and a theory of density-dependent effective electrostatic interactions. Implementing the model using Monte Carlo simulation and thermodynamic perturbation theory, we compute equilibrium particle size distributions, swelling ratios, volume fractions, net valences, radial distribution functions, and static structure factors as functions of concentration. Trial Monte Carlo moves comprising particle displacements and size variations are accepted or rejected based on the total change in elastic and electrostatic energies. The theory combines first-order thermodynamic perturbation and variational free energy approximations. For illustrative system parameters, theory and simulation agree closely at concentrations ranging from dilute to beyond particle overlap. With increasing concentration, as microgels deswell, we predict a decrease in the net valence and an unusual saturation of pair correlations. Comparison with experimental data for deionized, aqueous suspensions of PNIPAM particles demonstrates the capacity of the coarse-grained model to predict and interpret measured swelling behavior.

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.

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.

Colloidal gas-liquid condensation of polystyrene latex particles with intermediate kappa a values (5 to 160, a >> kappa(-1))

Polystyrene latex particles showed gas-liquid condensation under the conditions of large particle radius (a >> kappa(-1)) and intermediate kappa a, where kappa is the Debye-Hückel parameter and a is the particle radius. The particles were dissolved in deionized water containing ethanol from 0 to 77 vol %, settled to the bottom of the glass plate within 1 h, and then laterally moved toward the center of a cell over a 20 h period in reaching a state of equilibrium condensation. All of the suspensions that were 1 and 3 microm in diameter and 0.01-0.20 vol % in concentration realized similar gas-liquid condensation with clear gas-liquid boundaries. In 50 vol % ethanol solvent, additional ethanol was added to enhance the sedimentation force so as to restrict the particles in a monoparticle layer thickness. The coexistence of gas-liquid-solid (crystalline solid) was microscopically recognized from the periphery to the center of the condensates. A phase diagram of the gas-liquid condensation was created as a function of KCl concentration at a particle diameter of 3 microm, 0.10 vol % concentration, and 50:50 water/ethanol solvent at room temperature. The miscibility gap was observed in the concentration range from 1 to 250 microM. There was an upper limit of salt concentration where the phase separation disappeared, showing nearly critical behavior of macroscopic density fluctuation from 250 microM to 1 mM. These results add new experimental evidence to the existence of colloidal gas-liquid condensation and specify conditions of like-charge attraction between particles.

Phase diagrams of charged colloids from thermodynamic integration

We present full phase diagrams (including solid phases) of spherical charged colloids, using Monte Carlo sampling and thermodynamic integration of the Helmholtz free energy. Colloids and their co- and counterions are described by the primitive model for ionic systems that consists of hard-spheres with central point charges, while the solvent is taken into account solely through its dielectric constant. Two systems are considered: (i) a size-asymmetric system of oppositely charged spheres with size ratios q = 0.3 and 0.5 and (ii) a charge- and size-asymmetric system with colloid charge Q = 10 and counterions of charge -1 in the presence of monovalent added salt. In system (i), for both size ratios, the stable solid phase is equivalent to the NaCl crystal where the oppositely charged spheres take the lattice positions of Na and Cl ions. In system (ii), the phase diagram consists of gas-liquid and fluid-solid coexistence regions. We show that added salt stabilizes the fluid phase and shrinks the fluid-solid coexistence region, in agreement with experimental and theoretical results.

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.

Phase separation of charge-stabilized colloids: a Gibbs ensemble Monte Carlo simulation study

The fluid phase behavior of charge-stabilized colloidal suspensions is explored by applying a variant of the Gibbs ensemble Monte Carlo simulation method to a coarse-grained one-component model with implicit microions and solvent. The simulations take as input linear-response approximations for the effective electrostatic interactions--a hard-sphere-Yukawa pair potential and a one-body volume energy. The conventional Gibbs ensemble trial moves are supplemented by exchange of (implicit) salt between coexisting phases, with acceptance probabilities influenced by the state dependence of the effective interactions. Compared with large-scale simulations of the primitive model, with explicit microions, our computationally practical simulations of the one-component model closely match the pressures and pair distribution functions at moderate electrostatic couplings. For macroion valences and couplings within the linear-response regime, deionized aqueous suspensions with monovalent microions exhibit separation into macroion-rich and macroion-poor fluid phases below a critical salt concentration. The resulting pressures and phase diagrams are in excellent agreement with predictions of a variational free energy theory based on the same model.

Screening effects on structure and diffusion in confined charged colloids

Using molecular dynamics computer simulations we investigate structural and dynamic (diffusion) properties of charged colloidal suspension confined to narrow slit pores with structureless, uncharged walls. The system is modeled on an effective level involving only the macroions, which interact via a combination of a soft-sphere and a screened Coulomb potential. The aim of our study is to identify the role of the range of the macroion-macroion interaction controlled by the inverse Debye screening length, kappa. We also compare to bulk properties at the same chemical potential as determined in parallel grand canonical Monte Carlo simulations. Our results reveal a significant influence of the interaction range which competes, however, with the influence of density. At liquidlike densities a decrease of range yields a decreasing mobility (and a corresponding enhancement of local structure) in the bulk system, whereas the reverse effect occurs in narrow slits with thickness of a few particle diameter. These differences can be traced back to the confinement-induced, and kappa-dependent, reduction of overall density compared to the bulk reservoir. We also show that an increase of kappa softens the oscillations in the normal pressure as function of the wall separation, which is consistent with experimental observations concerning the influence of addition of salt.

Nonlinear screening and gas-liquid separation in suspensions of charged colloids

We calculate phase diagrams of charged colloidal spheres (valency Z and radius a) in a 1:1 electrolyte from multicentered nonlinear Poisson-Boltzmann theory. Our theory takes into account charge renormalization of the colloidal interactions and volume terms due to many-body effects. For valencies as small as Z = 1 and as large as 10(4) we find a gas-liquid spinodal instability in the colloid-salt phase diagram provided Z lambdaB/a > or similar 24+/-1, where lambdaB is the Bjerrum length.

Is perturbation DFT approach applicable to purely repulsive fluids?

A recently proposed third order + second order perturbation density functional theory (DFT) approach is tested for the validity and applicability to purely repulsive model fluids subjected to various external fields. Hard core repulsive Yukawa potential, point particle Yukawa potential, and inverse power potential are employed as sample models. Theoretical DFT results are compared with the corresponding simulation data obtained by grand canonical ensemble Monte Carlo simulation. This comparison indicates that the third order + second order perturbation DFT approach is suitable for these purely repulsive fluids only on condition of high accuracy of the imported bulk second order direct correlation function (DCF). However, in this case the origin of the successful performance somewhat differs from that observed for the mean field approximation applied to van der Waals fluids. In the present case it originates from the observation that the bulk second order DCF is strongly dependent on the density argument for the hard-core part, while for the distances exceeding the core dimension this dependence is considerably weaker.