Electrostatic interaction, calculation

The nonbonded terms comprise the Van der Waal s (VdW) interactions (calculated from the Lennard-Jones potential), the electrostatic interaction (calculated from Coulomb s law) and eventually the H-bond interaction (calculated by a modified Lennard-Jones potential). 

B. A. Luty, I. G. Tironi, and W. F. van Gunsteren. Lattice-sum methods for calculating electrostatic interactions in molecular simulations. J. Chem. Phys., 103 3014-3021, 1995. 

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. 

Brock A. Luty and Wilfried F. van Gunsteren. Calculating electrostatic interactions using the particle-particle particle-mesh method with nonperiodic long-range interactions. J. Phys. Chem., 100 2581-2587, 1996. 

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]. 

Electrostatic terms other than the simple charge interactions above are commonly included in molecular mechanics calculations. particularly dipole-dipole interactions. More recently, second-order electrostatic interactions like those describing polarizability have been added to some force fields. 

The first modification is to simply scale the dielectric permittivity of free space (T( ) by a scale factorD to rn ediate or dam pen thelong range electrostatic interactions. Its value was often set to be between 1.0 and 7H.0, the macroscopic value for water. A value of D=2..5, so that u=2..5Ug, wasoften used in early CIIARMM calculation s. 

A central multipole expansion therefore provides a way to calculate the electrostatic interaction between two molecules. The multipole moments can be obtained from the wave-function and can therefore be calculated using quantum mechanics (see Section 2.7.3) or can be determined from experiment. One example of the use of a multipole expansion is... 

VVe therefore return to the point-charge model for calculating electrostatic interactions. If sufficient point charges are used then all of the electric moments can be reproduced and the multipole interaction energy. Equation (4.30), is exactly equal to that calculated from the Coulomb summation. Equation (4.19). 

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. 

The conductor-like screening model (COSMO) is a continuum method designed to be fast and robust. This method uses a simpler, more approximate equation for the electrostatic interaction between the solvent and solute. Line the SMx methods, it is based on a solvent accessible surface. Because of this, COSMO calculations require less CPU time than PCM calculations and are less likely to fail to converge. COSMO can be used with a variety of semiempirical, ah initio, and DFT methods. There is also some loss of accuracy as a result of this approximation. 

The alkali metals tend to ionize thus, their modeling is dominated by electrostatic interactions. They can be described well by ah initio calculations, provided that diffuse, polarized basis sets are used. This allows the calculation to describe the very polarizable electron density distribution. Core potentials are used for ah initio calculations on the heavier elements. 

Another difference between the force fields is the calculation of electrostatic interactions. AMBER, BIO+, and OPLS use point charges to model electrostatic interactions. MM+ calculates electrostatic interactions using bond dipoles. The bond dipole method may not adequately simulate very polar or charged systems. 

---