Calculations of Double-Layer Electrostatic Interactions for the Sphere/Plane Geometry

Area of Science:

• Colloid and Surface Science
• Computational Electrochemistry
• Physical Chemistry

Background:

• The nonlinear Poisson-Boltzmann equation is crucial for understanding electrostatic interactions in electrolyte solutions.
• Accurate numerical solutions are needed for complex geometries like sphere/sphere and plane/sphere.
• Existing methods may lack efficiency or precision in resolving particle vicinity interactions.

Purpose of the Study:

• To develop and validate a novel numerical scheme for the nonlinear Poisson-Boltzmann equation.
• To investigate electrostatic interactions between charged particles in electrolyte solutions.
• To compare numerical results with analytical approximations for interaction energy.

Main Methods:

• Developed a numerical scheme using alternating direction overrelaxation and Newton-Raphson iteration.
• Employed grid transforming functions for improved mesh point distribution near particles.
• Calculated electric potential, force, and interaction energy for various boundary conditions (constant potential, constant charge, mixed).

Main Results:

• The numerical scheme successfully solved the nonlinear Poisson-Boltzmann equation for specified geometries.
• Electric potential distributions, interaction forces, and energies were computed.
• Calculated energy profiles were compared with analytical methods like Hogg-Healy-Fuerstenau and linear superposition approach (LSA).

Conclusions:

• The developed numerical method provides accurate solutions for electrostatic interactions.
• The linear superposition approach (LSA) is a reliable approximation for interaction energy at distances greater than kappa-1.
• The findings are pertinent to understanding colloid particle adsorption phenomena.