A simple yet accurate boundary element method for continuum dielectric calculations

Abstract

AbstractA simple yet accurate method for calculating electrostatic potentials using the boundary element continuum dielectric method is presented. It is shown that the limiting factor in accuracy is not the evaluation of integrals involving the interaction between boundary elements but rather a proper estimation of the self‐polarization of a patch upon itself. We derive a sum rule that allows us to calculate this important self‐polarization term in a self‐consistent and simple way. Intricate integration schemes used in previous treatments are consequently rendered unnecessary while concurrently achieving at least comparable accuracy over earlier methods. In some model systems for which analytic solutions are available, the computed surface polarization charge and reaction field energy are correct to better than six significant figures. An application of the method to the calculation of hydration free energies is presented. Good agreement with experimental values is obtained.