Numerical Studies on Electrical Interaction Forces and Free Energy between Colloidal Plates of Finite Size.

Abstract

By solving the nonlinear Poisson-Boltzmann (PB) equation with a finite element method (FEM), three-dimensional (3D) spatial distributions of the electric potential (ψ, scaled) in electrolyte solutions having two charged parallel finite plates (including cubes and prismatic rods) are determined for various separations (d, scaled by the Debye length, κ-1), surface potentials (ψs), and plate dimensions (length × width × thickness, each scaled by κ-1). The total interaction force between two plates, F, is the sum of the electrostatic double-layer (EDL) repulsion (the osmotic pressure, Fosm) and the Maxwell electrostatic stress (Fes). The EDL repulsion is estimated using the distribution of ψ not only between the facing surfaces of two parallel plates but also around the other extremities of the plates. The Maxwell stress (Fes) is localized near the extremities to act as a repulsive force on the midplane between the two plates. The ratio Fes/F is 0.07-0.5, depending on d, ψs, and dimensions. It is found that, with increasing dimensions, the total F values per unit area calculated for finite plates, F̃, decreasingly approach the exact ones for parallel infinite plates, F̃inf; for example, at d = 1 and ψs = 5, the ratio F̃/F̃inf is 2.83 for plates with dimensions of 1 × 1 × 1 and 1.18 for plates of 10 × 10 × 1. The repulsions arising from the extremities cannot be neglected for plates with dimensions <10 × 10 × 1. Furthermore, the total interaction forces (F) are calculated at a series of discrete d values, respectively, for parallel plates. We introduce a force fitting function, Ff(d), with parameters that can be determined so that Ff(d) fits well to the calculated serial F values. By integrating the Ff(d), we obtain the interaction free energy, G(d), for finite parallel plates that consists of two Γ functions.