Pressure calculation in polar and charged systems using Ewald summation: Results for the extended simple point charge model of water
- Institut fuer Experimentalphysik, Universitaet Wien, A-1090 Vienna (Austria)
Ewald summation and physically equivalent methods such as particle-mesh Ewald, kubic-harmonic expansions, or Lekner sums are commonly used to calculate long-range electrostatic interactions in computer simulations of polar and charged substances. The calculation of pressures in such systems is investigated. We find that the virial and thermodynamic pressures differ because of the explicit volume dependence of the effective, resummed Ewald potential. The thermodynamic pressure, obtained from the volume derivative of the Helmholtz free energy, can be expressed easily for both ionic and rigid molecular systems. For a system of rigid molecules, the electrostatic energy and the forces at the atom positions are required, both of which are readily available in molecular dynamics codes. We then calculate the virial and thermodynamic pressures for the extended simple point charge (SPC/E) water model at standard conditions. We find that the thermodynamic pressure exhibits considerably less system size dependence than the virial pressure. From an analysis of the cross correlation between the virial and thermodynamic pressure, we conclude that the thermodynamic pressure should be used to drive volume fluctuations in constant-pressure simulations. {copyright} {ital 1998 American Institute of Physics.}
- OSTI ID:
- 639163
- Journal Information:
- Journal of Chemical Physics, Vol. 109, Issue 7; Other Information: PBD: Aug 1998
- Country of Publication:
- United States
- Language:
- English
Similar Records
Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets
Temperature dependence of interactions of an ion pair in water: A molecular dynamics study