Real space electrostatics for multipoles. II. Comparisons with the Ewald sum.

Abstract

We report on tests of the shifted potential (SP), gradient shifted force (GSF), and Taylor shifted force (TSF) real-space methods for multipole interactions developed in Paper I of this series, using the multipolar Ewald sum as a reference method. The tests were carried out in a variety of condensed-phase environments designed to test up to quadrupole-quadrupole interactions. Comparisons of the energy differences between configurations, molecular forces, and torques were used to analyze how well the real-space models perform relative to the more computationally expensive Ewald treatment. We have also investigated the energy conservation, structural, and dynamical properties of the new methods in molecular dynamics simulations. The SP method shows excellent agreement with configurational energy differences, forces, and torques, and would be suitable for use in Monte Carlo calculations. Of the two new shifted-force methods, the GSF approach shows the best agreement with Ewald-derived energies, forces, and torques and also exhibits energy conservation properties that make it an excellent choice for efficient computation of electrostatic interactions in molecular dynamics simulations. Both SP and GSF are able to reproduce structural and dynamical properties in the liquid models with excellent fidelity.