Multi-multigrid solution of modified Poisson-Boltzmann equation for arbitrarily shaped molecules

A new multi-multigrid method is presented for solving the modified Poisson–Boltzmann equation based on the Kirkwood Hierarchy of equations, with Loeb's closure, on a three-dimensional grid. The results are compared with standard Poisson–Boltzmann calculations, which are known to underestimate the local concentration of counterions near charged parts of molecules, mainly due to neglect of fluctuations in the ionic concentrations. In the present study, the Kirkwood hierarchy of equations is discretized with the finite volume method and solved using multigrid techniques. The new possibility for solution of the three-dimensional modified Poisson–Boltzmann equation, for the first time within a model including a dielectric discontinuity, and within reasonable computational time, enables the calculation of higher valence ion distributions around arbitrarily shaped biological macromolecules.