Simple and robust solver for the Poisson-Boltzmann equation

Manuel Baptista, Roman Schmitz, B. Dünweg

Research output: Contribution to journalArticleResearchpeer-review

19 Citations (Scopus)


A variational approach is used to develop a robust numerical procedure for solving the nonlinear Poisson-Boltzmann equation. Following Maggs and Rossetto [Phys. Rev. Lett. 88, 196402 (2002)], we construct an appropriate constrained free energy functional such that its Euler-Lagrange equations are equivalent to the Poisson-Boltzmann equation. This is a formulation that searches for a true minimum in function space, in contrast to previous variational approaches that rather searched for a saddle point. We then develop, implement, and test an algorithm for its numerical minimization, which is quite simple and unconditionally stable. The analytic solution for planar geometry is used for validation. Some results are presented for a charged colloidal sphere surrounded by counterions and optimizations based upon fast Fourier transforms and hierarchical preconditioning are briefly discussed.

Original languageEnglish
Article number016705
JournalPhysical Review E
Issue number1
Publication statusPublished - 6 Aug 2009
Externally publishedYes

Cite this