Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions

Citations Over TimeTop 1% of 2001 papers

Abstract

A new version of the DelPhi program, which provides numerical solutions to the nonlinear Poisson−Boltzmann (PB) equation, is reported. The program can divide space into multiple regions containing different dielectric constants and can treat systems containing mixed salt solutions where the valence and concentration of each ion is different. The electrostatic free energy is calculated by decomposing the various energy terms into Coulombic interactions so that that the calculated free energies are independent of the lattice used to solve the PB equation. This, together with algorithms that optimally position polarization charges on the molecular surface, leads to a significant decrease in the dependence of the electrostatic free energy on the resolution of the lattice used to solve the PB equation and, hence, to a remarkable improvement in the precision of the calculated values. The Gauss−Seidel algorithm used in the current version of DelPhi is retained so that the new program retains many of the optimization features of the old one. The program uses dynamic memory allocation and can easily handle systems requiring large grid dimensionsfor example a 3003 system can be conveniently treated on a single SGI R12000 processor. An algorithm that estimates the best relaxation parameter to solve the nonlinear equation for a given system is described, and is implemented in the program at run time. A number of applications of the program are presented.

Related Papers

• → Numerical solution of a modified Poisson-Boltzmann equation in electric double layer theory(1979)22 cited
• → On the electrical double layer theory. I. A numerical method for solving a generalized Poisson—Boltzmann equation(1978)40 cited
• → Solution of the Poisson‐Boltzmann Equation about a Cylindrical Particle: Functional Theoretical Approach(2005)2 cited
• → Solution to the Linearized Poisson–Boltzmann Equation for a Spheroidal Surface under a General Surface Condition(1996)13 cited
• → A fast semi-numerical technique for the solution of the poisson-boltzmann equation in a cylindrical nanowire(2007)1 cited