Monte Carlo simulation for the potential of mean force between ionic colloids in solutions of asymmetric salts

A new technique for Monte Carlo sampling of the hard-sphere collision force has been applied to study the interaction between a pair of spherical macr oions in primitive-model electrolyte solutions with valences 1:2, 2:1, and 2:2. Macroions of the same charge can attract each other in the presence of divalent counterions, in analogy with earlier observations for planar and cylindrical geometries. The attraction is most significant at intermediate counterion concentrations. In contrast to the entropic depletion force betw een neutral particles, attraction between macroions is of energetic origin. The entropic contribution to the potential of mean force is generally repu lsive at conditions corresponding to aqueous colloids with or without salt. For systems with divalent counterions, the potentials of mean force predic ted by mean-field approximations like the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory or the Sogami-Ise (SI) theory are qualitatively different fro m those observed in the simulations. However, for systems with monovalent c ounterions, predictions of DLVO theory are in fair agreement with simulatio n results. (C) 1999 American Institute of Physics. [S0021-9606(99)50639-9].