Electrostatics Ewald

Belhadj,M., Alper, H.A., Levy, R.M. Molecular dynamics simulations of wa ter with Ewald summation for the long-range electrostatic interactions. Chem. Phys. Lett. 179 (1991) 13-20. 

U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. The smooth particle mesh ewald method. J. Chem. Phys., 103 8577, 1995. Brock A. Luty, Ilario G. Tironi, and Wilfried F. van Gunsteren. Lattice-sum methods for calculating electrostatic interactions in molecular simulations. J. Chem. Phys., 103 3014-3021, 1995. 

There are three different algorithms for the calculation of the electrostatic forces in systems with periodic boundary conditions (a) the (optimized) Ewald method, which scales like (b) the Particle Mesh... 

In periodic boimdary conditions, one possible way to avoid truncation of electrostatic interaction is to apply the so-called Particle Mesh Ewald (PME) method, which follows the Ewald summation method of calculating the electrostatic energy for a number of charges [27]. It was first devised by Ewald in 1921 to study the energetics of ionic crystals [28]. PME has been widely used for highly polar or charged systems. York and Darden applied the PME method already in 1994 to simulate a crystal of the bovine pancreatic trypsin inhibitor (BPTI) by molecular dynamics [29]. 

Luty B A, M E David, I G Tironi and W F van Gunsteren 1994. A Comparison of Particle-Particle, Particle-Mesh and Ewald Methods for Calculating Electrostatics Interactions in Periodic Molecular Systems. Molecular Simulation 14 11-20. 

Thus, for a point r in the central cell that does not coincide with any atomic position Vi, i = 1,. . ., N, the electrostatic potential ( )(r) in Eq. (19) can be rewritten in the Ewald formulation as ( )i(r) + ( )2(r). The electrostatic potential at atom i is the potential due to all other atoms j together with their images as well as all nontrivial periodic images of atom i itself. This is like the potential ( )(r) except that the (infinite) potential due to i itself is missing. Thus the potential at i can be obtained by removing the potential qj ri — r ... 

Calculated Electrostatic Energy of the Site(eV) Ewald Summation 977 Pt. Ion Cluster... 

Variational electrostatic projection method. In some instances, the calculation of PMF profiles in multiple dimensions for complex chemical reactions might not be feasible using full periodic simulation with explicit waters and ions even with the linear-scaling QM/MM-Ewald method [67], To remedy this, we have developed a variational electrostatic projection (VEP) method [75] to use as a generalized solvent boundary potential in QM/MM simulations with stochastic boundaries. The method is similar in spirit to that of Roux and co-workers [76-78], which has been recently... 

---