摘要

The linear Poisson-Boltzmann equation (LPBE) is one well-known implicit solvent continuum model for computing the electrostatic potential of biomolecules in ionic solvent. To overcome its singular difficulty caused by Dirac delta distributions of point charges and to further improve its solution accuracy, we develop in this paper a new scheme for solving the current LPBE model, a new LPBE model, and a new LPBE finite element program package based on our previously proposed PBE solution decomposition. Numerical tests on biomolecules and a nonlinear Born ball model with an analytical solution validate the new LPBE solution decomposition schemes, demonstrate the effectiveness and efficiency of the new program package, and confirm that the new LPBE model can significantly improve the solution accuracy of the current LPBE model.