Molecular dynamics simulations of pure water employing two different empirical water models have been used to study the effects of different methods for truncation of long-range interactions in molecular mechanics calculations. As has been observed previously in integral equation studies, "shifting"
Molecular dynamics simulations of water with Ewald summation for the long range electrostatic interactions
โ Scribed by Mahfoud Belhadj; Howard E. Alper; Ronald M. Levy
- Publisher
- Elsevier Science
- Year
- 1991
- Tongue
- English
- Weight
- 713 KB
- Volume
- 179
- Category
- Article
- ISSN
- 0009-2614
No coin nor oath required. For personal study only.
โฆ Synopsis
The use of a partial charge version of the Ewald sum method to treat long range electrostatic interactions in molecular dynamics simulations of water has been investigated. The orientational structure and energetic properties of the liquid have been studied for several choices of Ewald sum parameters. The results, including the dielectric constant c,,, the radial-dependent Kirkwood factor G,(r), the orientational projection Hd( I), and the average dipole interaction energy, are compared with those from molecular dynamics simulations of water using the reaction field method. While the orientational structure of the liquid obtained from simulations which treat long range interactions by Ewald summation is different from corresponding results obtained from simulations with a reaction field, consistent values are calculated for the dielectric constant. The sensitivity ofthe results to variations in the convergence parameter (Y, the real-space cutoff 4, and the number of vectors used in the reciprocal space sum was also analyzed.
๐ SIMILAR VOLUMES
We have used molecular dynamics simulations to study the physical properties of modified TIP3P water model included in the CHARMM program, using four different methods-the Ewald summation technique, and three different spherical truncation methods-for the treatment of the long-range interactions. Bo
## Method of solution University of Belfast, N. Ireland (see application form in this A system of a few hundred molecules is simulated. The equaissue) tions of motion are solved in a series of time steps, at each of which the force on each mole~ule, due to its interaction with Computer: CDC 7600;I
The deuterium quadrupole coupling constant and asymmetry parameter in heavy water were determined using ab initio SCF calculations. Snapshots from a molecular dynamics simulation were used to give liquid water cluster configurations and the influence of simulation parameters on the quadrupole coupli