Staggered Mesh Ewald: An Extension of the Smooth Particle-Mesh Ewald Method Adding Great Versatility

Abstract

We draw on an old technique for improving the accuracy of mesh-based field calculations to extend the popular Smooth Particle Mesh Ewald (SPME) algorithm as the Staggered Mesh Ewald (StME) algorithm. StME improves the accuracy of computed forces by up to 1.2 orders of magnitude and also reduces the drift in system momentum inherent in the SPME method by averaging the results of two separate reciprocal space calculations. StME can use charge mesh spacings roughly 1.5× larger than SPME to obtain comparable levels of accuracy; the one mesh in an SPME calculation can therefore be replaced with two separate meshes, each less than one third of the original size. Coarsening the charge mesh can be balanced with reductions in the direct space cutoff to optimize performance: the efficiency of StME rivals or exceeds that of SPME calculations with similarly optimized parameters. StME may also offer advantages for parallel molecular dynamics simulations because it permits the use of coarser meshes without requiring higher orders of charge interpolation and also because the two reciprocal space calculations can be run independently if that is most suitable for the machine architecture. We are planning other improvements to the standard SPME algorithm, and anticipate that StME will work synergistically will all of them to dramatically improve the efficiency and parallel scaling of molecular simulations.

An Iron Reservoir to the Catalytic Metal

“…Long-range electrostatic interactions were calculated using the Particle mesh Ewald summation method (31), and a cut-off of 9 Å was used for non-bonded interactions. The simulations were performed using the NPT ensemble at a constant pressure of 1 bar, and a time step of 2 fs was used to integrate the equation of motion.…”

A GIPC1-Palmitate Switch Modulates Dopamine Drd3 Receptor Trafficking and Signaling

“…A series of coarse-grained MD simulations (CGMD) was carried out for a system composed of the transmembrane 7 (TM7)-helix-8 segment from Drd3. The CGMD simulations were done with the Martini force field, following established protocols (47)(48)(49). Briefly, an integration step of 0.04 ps was used, and long-range electrostatic and van der Waals interactions were calculated with shift functions.…”

Recent Advances and Future Trend on the Emerging Role of GPUs as Platforms for Virtual Screening-Based Drug Discovery

“…Previous works have investigated this approach in a similar fashion but for long range interactions using Ewald-Mesh methods (Cerutti et al, 2009). Given the molecular sizes involved in protein-ligand interactions, we are only interested in short-range electrostatics.…”