Damped reaction field method and the accelerated convergence of the real space Ewald summation

J Chem Phys. 2014 Oct 28;141(16):164108. doi: 10.1063/1.4898147.

Abstract

In this paper, we study a general theoretical framework which allows us to approximate the real space Ewald sum by means of effective force shifted screened potentials, together with a self term. Using this strategy it is possible to generalize the reaction field method, as a means to approximate the real space Ewald sum. We show that this method exhibits faster convergence of the Coulomb energy than several schemes proposed recently in the literature while enjoying a much more sound and clear electrostatic significance. In terms of the damping parameter of the screened potential, we are able to identify two clearly distinct regimes of convergence. First, a reaction field regime corresponding to the limit of small screening, where effective pair potentials converge faster than the Ewald sum. Second, an Ewald regime, where the plain real space Ewald sum converges faster. Tuning the screening parameter for optimal convergence occurs essentially at the crossover. The implication is that effective pair potentials are an alternative to the Ewald sum only in those cases where optimization of the convergence error is not possible.