A BOUNDARY ELEMENT METHOD FOR MOLECULAR ELECTROSTATICS WITH ELECTROLYTE EFFECTS

A boundary element method is developed to compute the electrostatic potential inside and around molecules in an electrolyte solution. A set of boundary integral equations are derived based on the integral formulations of the Poisson equation and the linearized Poisson‐Boltzmann equation. The boundary integral equations are then solved numerically after discretizing the molecular surface into a number of flat triangular elements. The method is applied to a spherical molecule for which analytical solutions are available. Use is made of both constant and linearly varying unknowns over the boundary elements, and the method is tested for various values of parameters such as the dielectric constant of the molecule, ionic strength, and the location of the interior point charge. The use of the boundary integral method incorporating the nonlinear Poisson‐Boltzmann equation is also briefly discussed. Copyright © 1990 John Wiley & Sons, Inc.