Search

Donnan equilibrium between a salt-free colloidal dispersion and an electrolyte solution has been investigated by Monte Carlo simulations. The Donnan potential is directly calculated by considering two compartments separated by a semipermeable membrane. In order to understand the role played by colloid–ion interactions, the influences of colloidal characteristics, including particle size R, intrinsic particle charge Z, counterion valency z c , and concentration c p , on Donnan potential Ψ D and effective charge Z eff are examined. Our simulations show that the electroneutrality condition is not followed in both compartments and the Donnan potential arises due to the net charge density. The Donnan potential grows by increasing c p and Z eff and by decreasing dielectric constant e r , i.e., Ψ D ∼ Z eff c p / e r approximately. Note that the effective charge varies with R , Z , c p , e r and z c as well. When the salt concentration is increased, the net charge density is lowered and thus the Donnan potential decays accordingly. The validity of the classical theory based on the Nernst equation and the electroneutrality assumption is also examined. In general, the simulation results at high-salt condition can be well represented by such mean-field theory.

The phenomenon of counterion condensation around a flexible polyelectrolyte chain with N monomers is investigated by Monte Carlo simulations in terms of the degree of ionization alpha, which is proportional to the effective charge. It is operationally defined as the ratio of observed to intrinsic counterion concentration, alpha = co/ci. The observed counterion concentration in the dilute polyelectrolyte solution is equivalent to an electrolyte solution of concentration co with the same counterion chemical potential. It can be determined directly by thermodynamic experiments such as ion-selective electrode. With the polyelectrolyte fixed at the center of the spherical Wigner-Seitz cell, the polymer conformation, counterion distribution, and chemical potential can be obtained. Our simulation shows that the degree of ionization rises as the polymer concentration decreases. This behavior is opposite to that calculated from the infinitely long charged rod model, which is often used to study counterion condensation. Moreover, we find that, for a specified line charge density, alpha decreases with an increment in chain length and chain flexibility. In fact, the degree of ionization is found to decline with increasing polymer fractal dimension, which can be tuned by varying bending modulus and solvent quality. Those results can be qualitatively explained by a simple model of two-phase approximation.

The sedimentation profile of a dilute colloidal solution follows the barometric distribution owing to the balance between gravitational force and thermal fluctuation. However, the electrostatic interactions may lead to significant deviation even in the low volume fraction limit (e.g., 10(-5)). On the basis of Monte Carlo simulations for a dilute, salt-free colloidal dispersion, five regimes can be identified through the resulting colloidal sedimentation profile and the counterion distribution. The electrostatic interactions depends on the Coulomb strength E(c) defined as the ratio of the Bjerrum length to the colloid size. At weak colloid-ion attractions (small E(c)), counterions tend to distribute uniformly in the container. However, both barometric and inflated profiles of colloids can be observed. On the contrary, at strong colloid-ion attraction (large E(c)), counterions accumulate in the vicinity of the colloids. Significant counterion condensation effectively decreases the strength of colloid-colloid repulsion and barometric profile of colloids can be obtained as well. As a result, the sedimentation profile and counterion distribution are indicative of the strength of effective colloid-colloid and colloid-ion interactions. It is also found that local electroneutrality condition is generally not satisfied and charge separation (or internal electric field) is neither a sufficient nor necessary condition for nonbarometric distributions.

At strong electrostatic coupling, counterions are accumulated in the vicinity of the surface of the charged particle with intrinsic charge Z. In order to explain the behavior of highly charged particles, effective charge Z(*) is therefore invoked in the models based on Debye-Huckel approximation, such as the Derjaguin-Landau-Verwey-Overbeek potential. For a salt-free colloidal suspension, we perform Monte Carlo simulations to obtain various thermodynamic properties omega in a spherical Wigner-Seitz cell. The effect of dielectric discontinuity is examined. We show that at the same particle volume fraction, counterions around a highly charged sphere with Z may display the same value of omega as those around a weakly charged sphere with Z(*), i.e., omega(Z)=omega(Z(*)). There exists a maximally attainable value of omega at which Z=Z(*). Defining Z(*) as the effective charge, we find that the effective charge passes through a maximum and declines again due to ion-ion correlation as the number of counterions is increased. The effective charge is even smaller if one adopts the Debye-Huckel expression omega(DH). Our results suggest that charge renormalization can be performed by chemical potential, which may be observed in osmotic pressure measurements.