THE FAST MULTIPOLE BOUNDARY-ELEMENT METHOD FOR MOLECULAR ELECTROSTATICS - AN OPTIMAL APPROACH FOR LARGE SYSTEMS

We propose a fast implementation of the boundary element method for solving the Poisson equation, which approximately determines the electrostatic field around solvated molecules of arbitrary shape. The method presented uses computational resources of order O(N) only, where N is the number of elements representing the dielectric boundary at the molecular surface. The method is based on the Fast Multipole Algorithm by Rokhlin and Greengard, which is used to calculate the Coulombic interaction between surface elements in linear time. We calculate the solvation energies of a sphere, a small polar molecule, and a moderately sized protein. The values obtained by the boundary element method agree well with results from finite difference calculations and show a higher degree of consistency due to the absence of grid dependencies. The boundary element method can be taken to a much higher accuracy than is possible with finite difference methods and can therefore be used to verify their validity. (C) 1995 by John Wiley and Sons, Inc.