The wall and multivalent counterion effects on the electrostatic force between like-charged spherical particles confined in a charged pore

The effect of wall confinement (wall charge and wall-sphere separation distance) on the electrostatic force between two charged spheres confined in a long charged pore in symmetric and asymmetric electrolytes have been quantified by solving the nonlinear Poisson-Boltzmann equation (PBE), using adaptive finite elements combined with error minimization techniques. The computed force indicated the strong effect of the wall potential on the reduction of the repulsive force for all type of electrolytes. The influence of the wall effect was reduced when the valence of the electrolyte was increased. A significant reduction in the repulsive force between the two spheres was also observed when the distance between the pore wall and the sphere surface was reduced. A smaller long-range repulsive interaction was observed between spheres when the solutions contained multivalent counterions as compared with a monovalent solution. However, at short ranges of separation distances multivalent counterions increase the electrostatic repulsive force between the spheres. The effect of the dimensionless radius of the spheres on the electrostatic force between them has been determined and a significant reduction observed as the dimensionless radius was reduced.     

The interaction between charged aggregates in electrolyte solution. A Monte Carlo simulation study

The Monte Carlo simulation technique has been used to calculate the electrostatic force acting upon a charged aggregate outside a similarly charged wall. Contrary to intuition and existing electrostatic theories, the force is found to be attractive for some realistic values of the parameters determining the system. High surface charge density, low temperature, low relative permitivity and polyvalent neutralizing counterions are all factors that favour a net attraction between the wall and the aggregate. In some cases the resulting electrostatic attractive force is found to be an order of magnitude larger than the ordinary van der Waals attraction applied in the DLVO theory of colloidal stability. The attractive interaction is interpreted as being due to currelations between the counterions in the electric double layers.     

Electrostatic double-layer interaction between spherical particles inside a rough capillary

Predictions of electrostatic double-layer interaction forces between two similarly charged spherical colloidal particles inside an infinitely long "rough" capillary are presented. A simple model of a rough cylindrical surface is proposed, which assumes the capillary wall to be a periodic function of axial position. The periodic roughness of the wall is characterized by the wavelength and amplitude of the undulations. The electrostatic double-layer interaction force between two spherical particles located axially inside this rough capillary is determined by solving the nonlinear Poisson-Boltzmann equation employing finite element analysis. The effect of surface roughness of the cylindrical enclosure on the interaction force between two particles is extensively studied on the basis of this model. The simulations are carried out for dimensionless amplitudes (amplitude/particle radii) ranging from 0.05 to 0.15 and scaled wavelengths (wavelength/particle radii) ranging from 0.4 to 4.0. The interaction force between the particles is significantly modified by the proximity of the rough capillary wall. Generally, the interaction force for rough capillaries oscillates around the corresponding interaction force in a smooth capillary depending on the magnitudes of the scaled amplitude and wavelength of the roughness. The influence of roughness on the electrostatic interactions becomes more pronounced when the surface potential of the cylinder wall is different from the sphere surface potentials. When the cylinder and the particle surfaces have large potential differences, the axial force experienced by a particle is dominated by the capillary roughness. There are dramatic oscillations of the force, which alternately becomes repulsive and attractive as the particle moves from the crest to the trough of the rough capillary wall. These results suggest that manipulation of colloidal particles in narrow microchannels may be subject to significant force variations owing to the roughness inherent in microfabricated channels etched on metal films.     

Analytical solution of Poisson–Boltzmann equation for interacting plates of arbitrary potentials and same sign

Efficient calculation of electrostatic interactions in colloidal systems is becoming more important with the advent of such probing techniques as atomic force microscopy. Such practice requires solving the nonlinear Poisson–Boltzmann equation (PBE). Unfortunately, explicit analytical solutions are available only for the weakly charged surfaces. Analysis of arbitrarily charged surfaces is possible only through cumbersome numerical computations. A compact analytical solution of the one-dimensional PBE is presented for two plates interacting in symmetrical electrolytes. The plates can have arbitrary surface potentials at infinite separation as long they have the same sign. Such a condition covers a majority of the colloidal systems encountered. The solution leads to a simple relationship which permits determination of surface potentials, surface charge densities, and electrostatic pressures as a function of plate separation H for different charging scenarios. An analytical expression is also presented for the potential profile between the plates for a given separation. Comparison of these potential profiles with those obtained by numerical analysis shows the validity of the proposed solution.     

Electrostatic interaction between two spherical colloidal particles

Explicit exact analytic expressions are obtained in the form of infinite series for the potential distribution and the potential energy of the electrostatic interaction for the system of two dissimilar spheres in an electrolyte solution on the basis of the linearized Poisson—Boltzmann equation without recourse to Derjaguin's approximation. The leading term of the expression for the interaction energy (the zeroth order approximation) corresponds to the interaction energy that would be obtained if both spheres were ion-penetrable spheres (“soft” spheres). This term is a screened Coulomb interaction due to a simple linear superposition of the unperturbed potentials of the respective spheres, which is proportional to the product of their unperturbed surface potentials. The first-order approximation corresponds to the interaction energy that would be obtained if either sphere were a soft particle (the other being hard). The first-order correction term consists of two sub-terms, each of which is proportional to the square of the unperturbed surface potential of either sphere and does not depend on the unperturbed surface potential of the other sphere, can be interpreted as the interaction between the soft sphere and its image with respect to the hard sphere. This image interaction is attractive if the surface potential of the hard sphere is constant and repulsive if the surface charge density of the sphere is constant. It is shown that Derjaguin's method as well as its extension to the interaction of unequal spheres by Hogg, Healy and Fuerstenau (HHF) is quite a good approximation.     

A Fast and Robust Poisson-Boltzmann Solver Based on Adaptive Cartesian Grids

An adaptive Cartesian grid (ACG) concept is presented for the fast and robust numerical solution of the 3D Poisson-Boltzmann Equation (PBE) governing the electrostatic interactions of large-scale biomolecules and highly charged multi-biomolecular assemblies such as ribosomes and viruses. The ACG offers numerous advantages over competing grid topologies such as regular 3D lattices and unstructured grids. For very large biological molecules and multi-biomolecule assemblies, the total number of grid-points is several orders of magnitude less than that required in a conventional lattice grid used in the current PBE solvers thus allowing the end user to obtain accurate and stable nonlinear PBE solutions on a desktop computer. Compared to tetrahedral-based unstructured grids, ACG offers a simpler hierarchical grid structure, which is naturally suited to multigrid, relieves indirect addressing requirements and uses fewer neighboring nodes in the finite difference stencils. Construction of the ACG and determination of the dielectric/ionic maps are straightforward, fast and require minimal user intervention. Charge singularities are eliminated by reformulating the problem to produce the reaction field potential in the molecular interior and the total electrostatic potential in the exterior ionic solvent region. This approach minimizes grid-dependency and alleviates the need for fine grid spacing near atomic charge sites. The technical portion of this paper contains three parts. First, the ACG and its construction for general biomolecular geometries are described. Next, a discrete approximation to the PBE upon this mesh is derived. Finally, the overall solution procedure and multigrid implementation are summarized. Results obtained with the ACG-based PBE solver are presented for: (i) a low dielectric spherical cavity, containing interior point charges, embedded in a high dielectric ionic solvent - analytical solutions are available for this case, thus allowing rigorous assessment of the solution accuracy; (ii) a pair of low dielectric charged spheres embedded in a ionic solvent to compute electrostatic interaction free energies as a function of the distance between sphere centers; (iii) surface potentials of proteins, nucleic acids and their larger-scale assemblies such as ribosomes; and (iv) electrostatic solvation free energies and their salt sensitivities - obtained with both linear and nonlinear Poisson-Boltzmann equation - for a large set of proteins. These latter results along with timings can serve as benchmarks for comparing the performance of different PBE solvers.     

Force measurements between Teflon AF and colloidal silica particles in electrolyte solutions

The interaction force between a very hydrophobic polymer surface and colloidal silica particles with a roughness of 10–15 nm has been measured in aqueous solutions of KOH and KCl using an atomic force microscope. The interaction can be described according to the DLVO theory by an electrical double-layer force that is repulsive at long distances and attractive at short distances and an attractive van der Waals force. The electrical double-layer potentials are compared to the zeta potentials of Teflon AF and the silica spheres. The roughness of the silica particles leads to an underestimation of the short-range attraction and the surface potential. Both KCl and KOH solutions affect the potential of the interacting surfaces. OH⁻ ions that adsorb preferentially to the Teflon AF surface create higher potentials than Cl⁻ ions. Range and strength of the attractive interaction are not affected by KCl solutions but reduced by addition of KOH. This can be explained by decreasing potential differences between the silica sphere and Teflon AF with increasing KOH concentration. In addition, the preferential adsorption of OH⁻ ions may lead to a reduction of the van der Waals interaction. The presence of nanobubbles, too, might play a role.     

Ion-mediated interactions between charged and neutral nanoparticles

Monte-Carlo simulations are used to study the ion-mediated effective interaction between weakly charged and highly charged nanoparticles in an implicit solvent. Three models of nanoparticles are successively studied, from crude charged hard spheres to dipolar and non-spherical nanoparticles. The analysis of the effective potential revealed that in an electrolyte solution, even a neutral nanoparticle feels an important repulsive force in the presence of a charged nanoparticle, with a typical range similar to the Debye length. When the two nanoparticles carry charges of opposite sign, we have shown that this repulsion can reverse the effect of the direct attractive electrostatic potential at short distances. This also yields the change of sign of the effective potential as a function of the relative orientations of two anisotropic nanoparticles. Moreover, we found that the 3-body terms of the effective potentials can overcome the 2-body terms, which is not observed in the case of symmetrically charged nanoparticles.     

Theory of capillary-induced interactions beyond the superposition approximation

Within a general theoretical framework, we study the effective, deformation-induced interaction between two colloidal particles trapped at a fluid interface in the regime of small deformations. In many studies, this interaction has been computed with the ansatz that the actual interface configuration for the pair is given by the linear superposition of the interface deformations around the single particles. Here, we assess the validity of this approach and compute the leading term of the effective interaction for a large interparticle separation beyond this so-called superposition approximation. As an application, we consider the experimentally relevant case of interface deformations owing to the electrostatic field emanating from charged colloidal particles. In mechanical isolation, i.e., if the net force acting on the total system consisting of the particles plus the interface vanishes, the superposition approximation is actually invalid. The effective capillary interaction is governed by contributions beyond this approximation and turns out to be attractive. For sufficiently small surface charges on the colloids, such that linearization is strictly valid, and at asymptotically large separations, the effective interaction does not overcome the direct electrostatic repulsion between the colloidal particles.     

Sedimentation of a charged colloidal sphere in a charged cavity

An analytical study is presented for the quasisteady sedimentation of a charged spherical particle located at the center of a charged spherical cavity. The overlap of the electric double layers is allowed, and the polarization (relaxation) effect in the double layers is considered. The electrokinetic equations that govern the ionic concentration distributions, electric potential profile, and fluid flow field in the electrolyte solution are linearized assuming that the system is only slightly distorted from equilibrium. Using a perturbation method, these linearized equations are solved for a symmetric electrolyte with the surface charge densities of the particle and cavity as the small perturbation parameters. An analytical expression for the settling velocity of the charged sphere is obtained from a balance among the gravitational, electrostatic, and hydrodynamic forces acting on it. Our results indicate that the presence of the particle charge reduces the magnitude of the sedimentation velocity of the particle in an uncharged cavity and the presence of the fixed charge at the cavity surface increases the magnitude of the sedimentation velocity of an uncharged particle in a charged cavity. For the case of a charged sphere settling in a charged cavity with equivalent surface charge densities, the net effect of the fixed charges will increase the sedimentation velocity of the particle. For the case of a charged sphere settling in a charged cavity with their surface charge densities in opposite signs, the net effect of the fixed charges in general reduces/increases the sedimentation velocity of the particle if the surface charge density of the particle has a greater/smaller magnitude than that of the cavity. The effect of the surface charge at the cavity wall on the sedimentation of a colloidal particle is found to increase with a decrease in the particle-to-cavity size ratio and can be significant in appropriate situations.     

Strong and weak adsorptions of polyelectrolyte chains onto oppositely charged spheres

We investigate the complexation of long thin polyelectrolyte (PE) chains with oppositely charged spheres. In the limit of strong adsorption, when strongly charged PE chains adapt a definite wrapped conformation on the sphere surface, we analytically solve the linear Poisson-Boltzmann equation and calculate the electrostatic potential and the energy of the complex. We discuss some biological applications of the obtained results. For weak adsorption, when a flexible weakly charged PE chain is localized next to the sphere in solution, we solve the Edwards equation for PE conformations in the Hulthen potential, which is used as an approximation for the screened Debye-Huckel potential of the sphere. We predict the critical conditions for PE adsorption. We find that the critical sphere charge density exhibits a distinctively different dependence on the Debye screening length than for PE adsorption onto a flat surface. We compare our findings with experimental measurements on complexation of various PEs with oppositely charged colloidal particles. We also present some numerical results of the coupled Poisson-Boltzmann and self-consistent field equation for PE adsorption in an assembly of oppositely charged spheres.     

Depletion potentials induced by charged colloidal rods

We present direct depletion potential measurements for a single colloidal sphere close to a wall in suspensions of charged colloidal rods. In contrast to earlier studies of purely entropic systems (Helden et al. Phys. Rev. Lett. 2003, 90, 048301), here electrostatic interactions are important. These enhance the depletion attraction and lead to repulsive parts in the interaction potentials, indicating correlation effects between the rods.     

DLVO interactions of tungsten oxide and cobalt oxide surfaces measured with the colloidal probe technique

We have investigated the DLVO surface forces of oxidized tungsten and cobalt surfaces using the atomic force microscope (AFM) colloidal probe technique. It was shown by X-ray photoelectron spectroscopy (XPS) and electrokinetic measurements that this model system is representative of industrial tungsten carbide (WC) and cobalt powders used in the production of hard metals. We found that the attractive van der Waals forces are well described by Hamaker constants, calculated from optical data for WO(3) and CoOOH. The repulsive electrostatic double layer forces between WO(3) surfaces increase with increasing pH due to an increasingly negative surface potential. This surface potential decreases with increasing ionic strength at pH 7.5. The electrostatic interaction between WO(3) and CoOOH is attractive at pH 10, suggesting a positively charged CoOOH surface.     

Like‐Charge Attraction: A Combined Electrostatic–Hydrodynamic Approach

Abstract

This article discusses the like‐charge attraction of colloidal spheres close to a charged plate and compares results produced by an electrostatic and a hydrodynamic model with experimental data. Hydrodynamic coupling is shown to be the dominating effect, while the electrostatic influence may often be neglected. Some observations, however, can be explained only by means of a combined electrostatic–hydrodynamic model, which is derived in this work. The combined model is able to predict not only the attractive force between particles of similar charge close to a charged plate but also the change to a purely repulsive force once the sphere‐plate distance is further reduced. This prediction matches qualitatively results of experiments reported in the literature.     

The influence of discrete surface charges on the force between charged surfaces

The force between two parallel charged flat surfaces, with discrete surface charges, has been calculated with Monte Carlo simulations for different values of the electrostatic coupling. For low electrostatic coupling (small counterion valence, small surface charge, high dielectric constant, and high temperature) the total force is dominated by the entropic contribution and can be described by mean field theory, independent of the character of the surface charges. For moderate electrostatic coupling, counterion correlation effects lead to a smaller repulsion than predicted by mean field theory. This correlation effect is strengthened by discrete surface charges and the repulsive force is further reduced. For large electrostatic coupling the total force for smeared out surface charges is known to be attractive due to counterion correlations. If discrete surface charges are considered the attractive force is weakened and can even be turned into a repulsive force. This is due to the counterions being strongly correlated to the discrete surface charges forming effective, oppositely directed, dipoles on the two walls.     

Interaction of similarly charged surfaces mediated by nanoparticles

We study the interaction between two like charged surfaces embedded in a solution of oppositely charged multivalent rod-like counterions.The counterions consist of two rigidly bonded point charges,each of valency Z.The strength of the electrostatic coupling increases with increasing surface charge density or valency of the charges.The system is analyzed by employing a self-consistent field theory,which treats the short and long range interactions of the counterions within different approximations.We find that in the weak coupling limit,the interactions are only repulsive.In the intermediate coupling regime,the multivalent rod-like counterions can mediate attractive interactions between the surfaces. For sufficiently long rods,bridging contributes to the attractive interaction.In the strong coupling limit,the charge correlations can contribute to the attractive interactions at short separations between the charged surfaces.Two minima can then appear in the force curve between surfaces.     

Electrostatic surface interactions in mixtures of symmetric and asymmetric electrolytes: a Monte Carlo study

Canonical Monte Carlo simulations of the interaction between a uniformly charged spherical particle and a discretely charged planar surface in solutions of symmetric and asymmetric electrolytes were performed. To assess the nature of the interactions, the force exerted on the colloidal particle perpendicular to the planar surface was calculated. Attractive minima in the interaction force between the similarly charged surfaces reveal the occurrence of two phenomena: long-range attraction of electrostatic origin and short-range attraction due to depletion effects. The degree of electrostatic coupling determines the magnitude and range of like-charge attraction between the two surfaces.     

带同种电荷胶体颗粒间的相互吸引

带同种电荷颗粒间存在静电吸引是当前胶体科学领域的一大研究热点, 也是颇有争议的一个课题。经典的DLVO 理论认为, 电解质溶液中一对孤立、带同种电荷的球型胶粒将表现出纯粹的库仑排斥。以往的实验事实也证实了DLVO 理论的正确性。然而, 近十几年来的实验研究表明, 在一定条件下, 同性颗粒间的有效相互作用势能可能表现出静电吸引。对此, 至今还没有一个令人信服的理论描述, 甚至连其产生根源也未有定论。带同种电荷颗粒的静电吸引将是一个很值得探讨的问题。     

Effective interaction of nonuniformly charged colloid spheres in a bulk electrolyte

The interaction energy between two like-charged colloid spheres with nonuniformly distributed surface charges immersed in a bulk electrolyte is calculated under the linearized Poisson-Boltzmann equation. The calculated results are sensitive to the relative orientations of the spheres due to the anisotropy of the distribution of charges on the sphere surface. The effective repulsive interaction after thermal average is weaker than the case when charges are uniformly distributed on the spheres, and with the increase of the nonuniformity, an attractive interaction between the two spheres emerges.