Coupling between bulk- and surface chemistry in suspensions of charged colloids

Electrokinetic effects of ambient and excess carbonization of dielectric surfaces in aqueous environments

The charge state of surfaces in contact with aqueous electrolytes is crucial for the performance and stability of dielectric surfaces in general and lyophobic colloids in particular. Thus far the role of adsorbed molecular CO2 remained largely unexplored. The aim of the present investigation is to study the de-charging and re-charging for two model surfaces upon addition of CO2 and/or 1:1 electrolytes (NaCl, HCl) under precisely controlled boundary conditions up to millimolar concentrations of additives. Starting from the salt- and CO2-free state, the ζ-potential magnitudes drop linearly with the logarithm of the CO2-concentrations over several orders of magnitude in CO2-concentrations. Hydrophobic Polystyrene nearly fully discharges, hydrophilic SiO2 reveals a 60% charge reduction. From the surface specific effects of instead adding NaCl or HCl, we discriminate and parameterize empirically the relative contribution of three individual mechanisms for decreasing the ζ-potential magnitudes (screening, pH-driven charge regulation, dielectric charge regulation) combining during CO2-addition. Moreover, depending on the achieved CO2-induced de-charging, the behavior upon subsequent addition of NaCl and HCl switches between two limiting cases. Screening dominates for surfaces in the native state without CO2, but a significant re-charging is observed for surfaces conditioned under excess CO2-concentrations.

Charge fluctuations in charge-regulated systems: dependence on statistical ensemble

We investigate charge regulation of nanoparticles in concentrated suspensions, focusing on the effect of different statistical ensembles. We find that the choice of ensemble does not affect the mean charge of nanoparticles, but significantly alters the magnitude of its fluctuation. Specifically, we compared the behaviors of colloidal charge fluctuations in the semi-grand canonical and canonical ensembles and identified significant differences between the two. The choice of ensemble-whether the system is isolated or is in contact with a reservoir of acid and salt-will, therefore, affect the Kirkwood-Shumaker fluctuation-induced force inside concentrated suspensions. Our results emphasize the importance of selecting an appropriate ensemble that accurately reflects the experimental conditions when studying fluctuation-induced forces between polyelectrolytes, proteins, and colloidal particles in concentrated suspensions.

Surface charge regulation using classical density functional theory: the effect of divalent potential determining ions

The charge regulation approach has been used to describe the charge of surfaces susceptible to the presence of protons and other ions. Conventionally, this model is used with the Poisson-Boltzmann equation, which generally neglects the finite size of the ions and the electrostatic correlations. Recently, progress has been made by coupling charge regulation with classical density functional theory (DFT), which explicitly includes these correlations. However, little is known about charge regulation at surfaces with both acid-base equilibria and complexation with multivalent ions. The main purpose of this work is to investigate the role divalent ions play in charge regulation. Using DFT, we show that the size of the divalent ion has significant consequences on the surface charge density and it should not be neglected. For the surface reactions investigated, the larger the size of the divalent cation, the greater the charge on the surface due to higher divalent concentration there. At low divalent concentration, the ion correlations play a second-order but non-negligible role; using Poisson-Boltzmann theory with point ions cannot recover the DFT surface charge. At high concentrations, ion correlations play a dominant role by creating charge inversion.

Structural correlations in highly asymmetric binary charged colloidal mixtures

We explore structural correlations of strongly asymmetric mixtures of binary charged colloids within the primitive model of electrolytes considering large charge and size ratios of 10 and higher. Using computer simulations with explicit microions, we obtain the partial pair correlation functions between the like-charged colloidal macroions. Interestingly the big-small correlation peak amplitude is smaller than that of the big-big and small-small macroion correlation peaks, which is unfamiliar for additive repulsive interactions. Extracting optimal effective microion-averaged pair interactions between the macroions, we find that on top of non-additive Yukawa-like repulsions an additional shifted Gaussian attractive potential between the small macroions is needed to accurately reproduce their correct pair correlations. For small Coulomb couplings, the behavior is reproduced in a coarse-grained theory with microion-averaged effective interactions between the macroions. However, the accuracy of the theory deteriorates with increasing Coulomb coupling. We emphasize the relevance of entropic interactions exerted by the microions on the macroions. Our results are experimentally verifiable in binary mixtures of micron-sized colloids and like-charge nanoparticles.

Charge regulation of colloidal particles in aqueous solutions

We study the charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize a colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. On contact with water, these surface groups undergo proton transfer reactions, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. This theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.

Modeling Surface Charge Regulation of Colloidal Particles in Aqueous Solutions

Colloidal particles are mostly charged in an aqueous solution because of the protonation or deprotonation of ionizable groups on the surface. The surface charge density reflects a complex interplay of ion distributions within the electric double layer and the surface reaction equilibrium. In this work, we present a coarse-grained model to describe the charge regulation of various colloidal systems by an explicit consideration of the inhomogeneous ion distributions and surface reactions. With the primitive model for aqueous solutions and equilibrium constants for surface reactions as the inputs, the theoretical model is able to make quantitative predictions of the surface-charge densities and zeta potentials for diverse colloidal particles over a wide range of pH and ionic conditions. By accounting for the ionic size effects and electrostatic correlations, our model is applicable to systems with multivalent ions that exhibit charge inversion and provides a faithful description of the interfacial properties without evoking the empirical Stern capacitance or specific ion adsorptions.

Non-monotonic concentration dependence of the electro-phoretic mobility of charged spheres in realistic salt free suspensions

Using super-heterodyne Doppler velocimetry with multiple scattering correction, we extend the optically accessible range of concentrations in experiments on colloidal electro-kinetics. Here, we measured the electro-phoretic mobility and the DC conductivity of aqueous charged sphere suspensions covering about three orders of magnitude in particle concentrations and transmissions as low as 40%. The extended concentration range for the first time allows the demonstration of a non-monotonic concentration dependence of the mobility for a single particle species. Our observations reconcile previous experimental observations made on other species over restricted concentration ranges. We compare our results to the state-of-the-art theoretical calculations using a constant particle charge and the carefully determined experimental boundary conditions as input. In particular, we consider the so-called realistic salt free conditions, i.e., we respect the release of counterions by the particles, the solvent hydrolysis, and the formation of carbonic acid from dissolved neutral CO₂. We also compare our results to previous results obtained under similarly well-defined conditions. This allows identification of three distinct regions of differing density dependence. There is an ascent during the build-up of double layer overlap, which is not expected by theory, an extended plateau region in quantitative agreement with theoretical expectation based on a constant effective charge and a sudden decrease, which occurs way before the expected gradual decrease. Our observations suggest a relation of the non-monotonic behavior to a decrease in particle charge, and we tentatively discuss possibly underlying mechanisms.

Charge Regulation of Colloidal Particles: Theory and Simulations

To explore charge regulation (CR) in physicochemical and biophysical systems, we present a model of colloidal particles with sticky adsorption sites which account for the formation of covalent bonds between the hydronium ions and the surface functional groups. Using this model and Monte Carlo simulations, we find that the standard Ninham and Parsegian (NP) theory of CR leads to results which deviate significantly from computer simulations. The problem with the NP approach is traced back to the use of a bulk equilibrium constant to account for surface chemical reactions. To resolve this difficulty we present a new theory of CR. The fundamental ingredient of the new approach is the sticky length, which is nontrivially related to the bulk equilibrium constant. The theory is found to be in excellent agreement with computer simulations, without any adjustable parameters. As an application of the theory we calculate the effective charge of colloidal particles containing carboxyl groups, as a function of pH and salt concentration.

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.

Hydrodynamic simulations of charge-regulation effects in colloidal suspensions

Self-organization of charged soft matter is of crucial importance in biology. However, it is an extremely complex phenomenon due to dynamical couplings between the hydrodynamic flow, ions, and charges of soft matter. For colloidal suspensions, the coupling between the former two has already been studied by numerical simulations, with the colloid surface charge being fixed. However, the self-organization of colloids and/or the application of an external electric field make the electrostatic environment of each colloid inhomogeneous in both space and time. Thus, this leads to inhomogenisation of the surface charge of each colloid under ionisation equilibrium conditions. This effect is known as "charge regulation" and is of great importance in various electrostatic and electrokinetic phenomena not only in colloid suspensions but also in solutions of biomolecules. However, there has so far been no success in taking the charge regulation effect into account in numerical simulations of colloidal electrokinetics. Here, we extend the fluid particle dynamics (FPD) method to incorporate the charge regulation effect. We present a theoretical formulation of the method and its application to two types of problems, where charge regulation plays an important role: (i) cluster formation of colloid particles and (ii) a single colloid particle under an external field. By these examples, we show not only the importance of considering charge-regulation effects in the self-organization of charged systems but also the applicability of our simulation method to more complex problems of charged soft matter systems.

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.

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.

Ionization at a solid-water interface in an applied electric field: Charge regulation

We investigate ionization at a solid-water interface in an applied electric field. We attach an electrode to a dielectric film bearing silanol or carboxyl groups with an areal density Γ₀, where the degree of dissociation α is determined by the proton density in water close to the film. We show how α depends on the density n₀ of NaOH in water and the surface charge density σ_m on the electrode. For σ_m > 0, the protons are expelled away from the film, leading to an increase in α. In particular, in the range 0 < σ_m < eΓ₀, self-regulation occurs to realize α ≅ σ_m/eΓ₀ for n₀ ≪ n_c, where n_c is 0.01 mol/L for silica surfaces and is 2 × 10^-5 mol/L for carboxyl-bearing surfaces. We also examine the charge regulation with decreasing the cell thickness H below the Debye length κ^-1, where a crossover occurs at the Gouy-Chapman length. In particular, when σ_m ∼ eΓ₀ and H ≪ κ^-1, the surface charges remain only partially screened by ions, leading to a nonvanishing electric field in the interior.

Solvent Role in the Formation of Electric Double Layers with Surface Charge Regulation: A Bystander or a Key Participant?

The charge formation at interfaces involving electrolyte solutions is due to the chemical equilibrium between the surface reactive groups and the potential determining ions in the solution (i.e., charge regulation). In this Letter we report our findings that this equilibrium is strongly coupled to the precise molecular structure of the solution near the charged interface. The neutral solvent molecules dominate this structure due to their overwhelmingly large number. Treating the solvent as a structureless continuum leads to a fundamentally inadequate physical picture of charged interfaces. We show that a proper account of the solvent effect leads to an unexpected and complex system behavior that is affected by the molecular and ionic excluded volumes and van der Waals interactions.

Charge Regulation in the Electrical Double Layer: Ion Adsorption and Surface Interactions

Charge regulation in the electrical double layer has important implications for ion adsorption, interparticle forces, colloidal stability, and deposition phenomena. Although charge regulation generally receives little attention, its consequences can be major, especially when considering interactions between unequally charged surfaces. The present article discusses common approaches to quantify such phenomena, especially within classical Poisson-Boltzmann theory, and pinpoints numerous situations where a consideration of charge regulation is essential. For the interpretation of interaction energy profiles, we advocate the use of the constant regulation approximation, which summarizes the surface properties in terms of two quantities, namely, the diffuse layer potential and the regulation parameter. This description also captures some pronounced regulation effects observed in the presence of multivalent ions.

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.

Crystallization kinetics of colloidal model suspensions: recent achievements and new perspectives

Colloidal model systems allow studying crystallization kinetics under fairly ideal conditions, with rather well-characterized pair interactions and minimized external influences. In complementary approaches experiment, analytic theory and simulation have been employed to study colloidal solidification in great detail. These studies were based on advanced optical methods, careful system characterization and sophisticated numerical methods. Over the last decade, both the effects of the type, strength and range of the pair-interaction between the colloidal particles and those of the colloid-specific polydispersity have been addressed in a quantitative way. Key parameters of crystallization have been derived and compared to those of metal systems. These systematic investigations significantly contributed to an enhanced understanding of the crystallization processes in general. Further, new fundamental questions have arisen and (partially) been solved over the last decade: including, for example, a two-step nucleation mechanism in homogeneous nucleation, choice of the crystallization pathway, or the subtle interplay of boundary conditions in heterogeneous nucleation. On the other hand, via the application of both gradients and external fields the competition between different nucleation and growth modes can be controlled and the resulting microstructure be influenced. The present review attempts to cover the interesting developments that have occurred since the turn of the millennium and to identify important novel trends, with particular focus on experimental aspects.