Chemical Charge Regulation and Charge Renormalization in Concentrated Colloidal Suspensions

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.

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.

Charge regulating macro-ions in salt solutions: screening properties and electrostatic interactions

We revisit the charge-regulation mechanism of macro-ions and apply it to mobile macro-ions in a bathing salt solution. In particular, we examine the effects of correlation between various adsorption/desorption sites and analyze the collective behavior in terms of the solution effective screening properties. We show that such a behavior can be quantified in terms of the charge asymmetry of the macro-ions, defined by their preference for a non-zero effective charge, and their donor/acceptor propensity for exchanging salt ions with the bathing solution. Asymmetric macro-ions tend to increase the screening, while symmetric macro-ions can in some cases decrease it. Macro-ions that are classified as donors display a rather regular behavior, while those that behave as acceptors exhibit an anomalous non-monotonic Debye length. The screening properties, in their turn, engender important modifications to the disjoining pressure between two charged surfaces. Our findings are in particular relevant for solutions of proteins, whose exposed amino acids can undergo charge dissociation/association processes to/from the bathing solution, and can be considered as a solution of charged regulated macro-ions, as analyzed here.

A simple model for electrical charge in globular macromolecules and linear polyelectrolytes in solution

We present a model for calculating the net and effective electrical charge of globular macromolecules and linear polyelectrolytes such as proteins and DNA, given the concentration of monovalent salt and pH in solution. The calculation is based on a numerical solution of the non-linear Poisson-Boltzmann equation using a finite element discretized continuum approach. The model simultaneously addresses the phenomena of charge regulation and renormalization, both of which underpin the electrostatics of biomolecules in solution. We show that while charge regulation addresses the true electrical charge of a molecule arising from the acid-base equilibria of its ionizable groups, charge renormalization finds relevance in the context of a molecule’s interaction with another charged entity. Writing this electrostatic interaction free energy in terms of a local electrical potential, we obtain an “interaction charge” for the molecule which we demonstrate agrees closely with the “effective charge” discussed in charge renormalization and counterion-condensation theories. The predictions of this model agree well with direct high-precision measurements of effective electrical charge of polyelectrolytes such as nucleic acids and disordered proteins in solution, without tunable parameters. Including the effective interior dielectric constant for compactly folded molecules as a tunable parameter, the model captures measurements of effective charge as well as published trends of pKa shifts in globular proteins. Our results suggest a straightforward general framework to model electrostatics in biomolecules in solution. In offering a platform that directly links theory and experiment, these calculations could foster a systematic understanding of the interrelationship between molecular 3D structure and conformation, electrical charge and electrostatic interactions in solution. The model could find particular relevance in situations where molecular crystal structures are not available or rapid, reliable predictions are desired.

Strong coupling electrostatics for randomly charged surfaces: antifragility and effective interactions

We study the effective interaction mediated by strongly coupled Coulomb fluids between dielectric surfaces carrying quenched, random monopolar charges with equal mean and variance, both when the Coulomb fluid consists only of mobile multivalent counterions and when it consists of an asymmetric ionic mixture containing multivalent and monovalent (salt) ions in equilibrium with an aqueous bulk reservoir. We analyze the consequences that follow from the interplay between surface charge disorder, dielectric and salt image effects, and the strong electrostatic coupling that results from multivalent counterions on the distribution of these ions and the effective interaction pressure they mediate between the surfaces. In a dielectrically homogeneous system, we show that the multivalent counterions are attracted towards the surfaces with a singular, disorder-induced potential that diverges logarithmically on approach to the surfaces, creating a singular but integrable counterion density profile that exhibits an algebraic divergence at the surfaces with an exponent that depends on the surface charge (disorder) variance. This effect drives the system towards a state of lower thermal 'disorder', one that can be described by a renormalized temperature, exhibiting thus a remarkable antifragility. In the presence of an interfacial dielectric discontinuity, the singular behavior of counterion density at the surfaces is removed but multivalent counterions are still accumulated much more strongly close to randomly charged surfaces as compared with uniformly charged ones. The interaction pressure acting on the surfaces displays in general a highly non-monotonic behavior as a function of the inter-surface separation with a prominent regime of attraction at small to intermediate separations. This attraction is caused directly by the combined effects from charge disorder and strong coupling electrostatics of multivalent counterions, which dominate the surface-surface repulsion due to the (equal) mean charges on the two surfaces and the osmotic pressure of monovalent ions residing between them. These effects can be quite significant even with a small degree of surface charge disorder relative to the mean surface charge. The strong coupling, disorder-induced attraction is typically much stronger than the van der Waals interaction between the surfaces, especially within a range of several nanometers for the inter-surface separation, where such effects are predicted to be most pronounced.

Charging of poly(methyl methacrylate) (PMMA) colloids in cyclohexyl bromide: locking, size dependence, and particle mixtures

We studied suspensions of sterically stabilized poly(methyl methacrylate) (PMMA) particles in the solvent cyclohexyl bromide (CHB; εr = 7.92). We performed microelectrophoresis measurements on suspensions containing a single particle species and on binary mixtures, using confocal microscopy to measure the velocity profiles of the particles. We measured the charge of so-called locked PMMA particles, for which the steric stabilizer, a comb-graft stabilizer of poly(12-hydroxystearic acid) (PHSA) grafted on a backbone of PMMA, was covalently bonded to the particle, and for unlocked particles, for which the stabilizer was adsorbed to the surface of the particle. We observed that locked particles had a significantly higher charge than unlocked particles. We found that the charge increase upon locking was due to chemical coupling of 2-(dimethylamino)ethanol to the PMMA particles, which was used as a catalyst for the locking reaction. For particles of different size we obtained the surface potential and charge from the electrophoretic mobility of the particles. For locked particles we found that the relatively high surface potential (∼ +5.1 kBT/e or +130 mV) was roughly constant for all particle diameters we investigated (1.2 μm < σ < 4.4 μm), and that the particle charge was proportional to the square of the diameter.

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.

Beyond standard Poisson-Boltzmann theory: ion-specific interactions in aqueous solutions

The Poisson-Boltzmann mean-field description of ionic solutions has been successfully used in predicting charge distributions and interactions between charged macromolecules. While the electrostatic model of charged fluids, on which the Poisson-Boltzmann description rests, and its statistical mechanical consequences have been scrutinized in great detail, much less is understood about its probable shortcomings when dealing with various aspects of real physical, chemical and biological systems. These shortcomings are not only a consequence of the limitations of the mean-field approximation per se, but perhaps are primarily due to the fact that the purely Coulombic model Hamiltonian does not take into account various additional interactions that are not electrostatic in their origin. We explore several possible non-electrostatic contributions to the free energy of ions in confined aqueous solutions and investigate their ramifications and consequences on ionic profiles and interactions between charged surfaces and macromolecules.

Surfaces with quenched and annealed disordered charge distributions

We consider surfaces with disordered charge distribution. The disorder can be caused by mobile charges, as for example in mixed lipid bilayers, or by weakly charged surfaces where charge regulation takes place (e.g. carboxyl groups). Using Monte-Carlo simulation methods we find for quenched as well as annealed disordered charge distributions counterion densities close to the surface that are significantly larger than for ordered regularly spaced surface ions. Our field-theoretic results agree well with results obtained from Monte-Carlo simulations of the system. Furthermore, we obtain expressions for the effective interaction between charged colloids and charged rods close to a charged surface and discuss the effect of the surface-ion mobility and polarization charges on the interaction. In general, polarization effects as well as surface-ion mobility lead to a weakening of the effective interaction between charged objects.

Adsorption of barium and calcium chloride onto negatively charged alpha-Fe(2)O(3) particles

Adsorption of cations (Na(+), Ca(2+), Ba(2+)) onto negatively charged (pH 10.4) hematite (alpha-Fe(2)O(3)) particles has been studied. The oxide material was carefully prepared in order to obtain monodisperse suspensions of well-crystallized, quasi-spherical particles (50 nm in diameter). The isoelectric point (IEP) is located at pH 8.5. Adsorption of barium ions onto oxide particles was carried out and the electrophoretic mobility was measured throughout the adsorption experiment. Comparison with calcium adsorption at full coverage reveals a higher uptake of Ba(2+). In both cases it shows also that chloride ions coadsorb with M(2) ions. Simultaneous uptake of the positive and negative ions explains why the electrophoretic mobility does not reverse to cationic migration. A theoretical study of the surface speciation has been carried out, using the MuSiC model. It reveals the presence of negative as well as positive sites on both sides of the point of zero charge (PZC) of the hematite particles, which may explain the coadsorption of Ba(2+) and Cl(-) at pH 10.4. The effective charge of the oxide particles, calculated from the electrophoretic mobility, is in very good agreement with the results found with the MuSiC modelization and the chloride/barium adsorption ratio. It also verifies the theory of ionic condensation. Calorimetric measurements gave a negative heat for the overall reaction occurring when Ba(2+)/Cl(-) ions adsorb onto hematite. Despite the fact that anions (Cl(-) and OH(-)) adsorption onto mineral oxides is an exothermic phenomenon, it is likely that barium and calcium adsorption is endothermic, denoting the formation of an inner-sphere complex as reported in the literature.

Interaction potentials, structural ordering and effective charges in dispersions of charged colloidal particles

As colloidal dispersions of charged particles exhibit a wide variety of commercial, technological and scientific applications, a considerable theoretical effort has been devoted to finding an effective interaction potential from primitive models. The forces derived from this potential should justify the spatial ordering experimentally observed under certain conditions. This paper reviews the advances in these theoretical studies as well as some experiments (based on the mentioned order) that try to corroborate them. Special attention has been paid to the Derjaguin-Landau-Verwey-Overbeek (DLVO) potential. Nowadays, many of these theoretical investigations suggest that it could be applied if some of its parameters are renormalized. Nevertheless, to achieve a renormalization procedure in a strict way (from a primitive model) is a difficult task as a result of the size and charge asymmetries between small ions and macroions. Thus, several procedures for computing renormalized charges in a simple way have been developed. However, the notion of effective charge has also been widely used (as a adjustable parameter) in order to justify results found for several kinds of colloids (like solid particle dispersions or micellar systems) by means of quite different experimental techniques. Renormalization (as well as ion condensation) approaches, experiments and the controversial relationship between theoretical and phenomenological effective charges are also reviewed in this work.

Effective interactions between electric double layers

This review summarizes and assesses recent theoretical and experimental advances, with special emphasis on the effective interaction between charge-stabilized colloids, in the bulk or in confined geometries, and on the ambiguities of defining an effective charge of the colloidal particles. Some consideration is given to the often neglected discrete solvent effects.