Cutoff radius effect of the isotropic periodic sum and Wolf method in liquid-vapor interfaces of water

Kazuaki Z. Takahashi, Tetsu Narumi, Kenji Yasuoka

Research output: Contribution to journalArticle

36 Citations (Scopus)

Abstract

As a more economical but similarly accurate computation method than the Ewald sum, the isotropic periodic sum (IPS) method for nonpolar molecules (IPSn) and polar molecules (IPSp), along with the Wolf method are of interest, but the cutoff radius dependence is an important issue. To evaluate the cutoff radius effect of the three methods, a water-vapor interfacial system has been studied by molecular dynamics. The Wolf method can produce adequate results for surface tension compared to that of the Ewald sum (within 2.9) at a long enough cutoff radius, rc. However, the estimation of the electrostatic potential profile and dipole orientational function is poor. The Wolf method cannot estimate electrostatic configuration at rc Lz2 (L z is the longest lattice of the system). We have found that the convergence of the surface tension and the electrostatic configuration of the IPSn method is faster than that of the IPSp method. Moreover, the IPSn method is most accurate among the three methods for the same cutoff radius. Furthermore, the behavior of the surface tension against the cutoff radius shows a greater difference for the IPSn and IPSp method. The surface tension of the IPSp method fluctuates and presents a similar result to that of the Ewald sum, but the surface tension for the IPSn method greatly deviates near rc L z3. The cause of this deviation is the difference between the interfacial configuration of the water surface and the cutoff treatment of the IPS method. The deviation becomes insignificant far from rc L z3. In spite of this shortcoming, the IPSn method gives the most accurate result in estimating the surface tension at rc L z2. From all the results in this work, the IPSn and IPSp method have been found to be more accurate than the Wolf method. In conclusion, the surface tension and structure of water-vapor interface can be calculated by the IPSn method when rc is greater than or equal to the longest lattice of the system. The IPSp method and the Wolf method require a longer cutoff radius than the longest lattice of the system to estimate interfacial properties.

Original languageEnglish
Article number174112
JournalJournal of Chemical Physics
Volume134
Issue number17
DOIs
Publication statusPublished - 2011 May 7

Fingerprint

liquid-vapor interfaces
wolves
Surface tension
interfacial tension
cut-off
Vapors
Molecules
radii
Water
Liquids
water
molecules
Electrostatics
Steam
electrostatics
water vapor
configurations
deviation
estimates
surface water

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

Cutoff radius effect of the isotropic periodic sum and Wolf method in liquid-vapor interfaces of water. / Takahashi, Kazuaki Z.; Narumi, Tetsu; Yasuoka, Kenji.

In: Journal of Chemical Physics, Vol. 134, No. 17, 174112, 07.05.2011.

Research output: Contribution to journalArticle

@article{770d270953f74543b09a01d4f8732bea,
title = "Cutoff radius effect of the isotropic periodic sum and Wolf method in liquid-vapor interfaces of water",
abstract = "As a more economical but similarly accurate computation method than the Ewald sum, the isotropic periodic sum (IPS) method for nonpolar molecules (IPSn) and polar molecules (IPSp), along with the Wolf method are of interest, but the cutoff radius dependence is an important issue. To evaluate the cutoff radius effect of the three methods, a water-vapor interfacial system has been studied by molecular dynamics. The Wolf method can produce adequate results for surface tension compared to that of the Ewald sum (within 2.9) at a long enough cutoff radius, rc. However, the estimation of the electrostatic potential profile and dipole orientational function is poor. The Wolf method cannot estimate electrostatic configuration at rc Lz2 (L z is the longest lattice of the system). We have found that the convergence of the surface tension and the electrostatic configuration of the IPSn method is faster than that of the IPSp method. Moreover, the IPSn method is most accurate among the three methods for the same cutoff radius. Furthermore, the behavior of the surface tension against the cutoff radius shows a greater difference for the IPSn and IPSp method. The surface tension of the IPSp method fluctuates and presents a similar result to that of the Ewald sum, but the surface tension for the IPSn method greatly deviates near rc L z3. The cause of this deviation is the difference between the interfacial configuration of the water surface and the cutoff treatment of the IPS method. The deviation becomes insignificant far from rc L z3. In spite of this shortcoming, the IPSn method gives the most accurate result in estimating the surface tension at rc L z2. From all the results in this work, the IPSn and IPSp method have been found to be more accurate than the Wolf method. In conclusion, the surface tension and structure of water-vapor interface can be calculated by the IPSn method when rc is greater than or equal to the longest lattice of the system. The IPSp method and the Wolf method require a longer cutoff radius than the longest lattice of the system to estimate interfacial properties.",
author = "Takahashi, {Kazuaki Z.} and Tetsu Narumi and Kenji Yasuoka",
year = "2011",
month = "5",
day = "7",
doi = "10.1063/1.3578473",
language = "English",
volume = "134",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "17",

}

TY - JOUR

T1 - Cutoff radius effect of the isotropic periodic sum and Wolf method in liquid-vapor interfaces of water

AU - Takahashi, Kazuaki Z.

AU - Narumi, Tetsu

AU - Yasuoka, Kenji

PY - 2011/5/7

Y1 - 2011/5/7

N2 - As a more economical but similarly accurate computation method than the Ewald sum, the isotropic periodic sum (IPS) method for nonpolar molecules (IPSn) and polar molecules (IPSp), along with the Wolf method are of interest, but the cutoff radius dependence is an important issue. To evaluate the cutoff radius effect of the three methods, a water-vapor interfacial system has been studied by molecular dynamics. The Wolf method can produce adequate results for surface tension compared to that of the Ewald sum (within 2.9) at a long enough cutoff radius, rc. However, the estimation of the electrostatic potential profile and dipole orientational function is poor. The Wolf method cannot estimate electrostatic configuration at rc Lz2 (L z is the longest lattice of the system). We have found that the convergence of the surface tension and the electrostatic configuration of the IPSn method is faster than that of the IPSp method. Moreover, the IPSn method is most accurate among the three methods for the same cutoff radius. Furthermore, the behavior of the surface tension against the cutoff radius shows a greater difference for the IPSn and IPSp method. The surface tension of the IPSp method fluctuates and presents a similar result to that of the Ewald sum, but the surface tension for the IPSn method greatly deviates near rc L z3. The cause of this deviation is the difference between the interfacial configuration of the water surface and the cutoff treatment of the IPS method. The deviation becomes insignificant far from rc L z3. In spite of this shortcoming, the IPSn method gives the most accurate result in estimating the surface tension at rc L z2. From all the results in this work, the IPSn and IPSp method have been found to be more accurate than the Wolf method. In conclusion, the surface tension and structure of water-vapor interface can be calculated by the IPSn method when rc is greater than or equal to the longest lattice of the system. The IPSp method and the Wolf method require a longer cutoff radius than the longest lattice of the system to estimate interfacial properties.

AB - As a more economical but similarly accurate computation method than the Ewald sum, the isotropic periodic sum (IPS) method for nonpolar molecules (IPSn) and polar molecules (IPSp), along with the Wolf method are of interest, but the cutoff radius dependence is an important issue. To evaluate the cutoff radius effect of the three methods, a water-vapor interfacial system has been studied by molecular dynamics. The Wolf method can produce adequate results for surface tension compared to that of the Ewald sum (within 2.9) at a long enough cutoff radius, rc. However, the estimation of the electrostatic potential profile and dipole orientational function is poor. The Wolf method cannot estimate electrostatic configuration at rc Lz2 (L z is the longest lattice of the system). We have found that the convergence of the surface tension and the electrostatic configuration of the IPSn method is faster than that of the IPSp method. Moreover, the IPSn method is most accurate among the three methods for the same cutoff radius. Furthermore, the behavior of the surface tension against the cutoff radius shows a greater difference for the IPSn and IPSp method. The surface tension of the IPSp method fluctuates and presents a similar result to that of the Ewald sum, but the surface tension for the IPSn method greatly deviates near rc L z3. The cause of this deviation is the difference between the interfacial configuration of the water surface and the cutoff treatment of the IPS method. The deviation becomes insignificant far from rc L z3. In spite of this shortcoming, the IPSn method gives the most accurate result in estimating the surface tension at rc L z2. From all the results in this work, the IPSn and IPSp method have been found to be more accurate than the Wolf method. In conclusion, the surface tension and structure of water-vapor interface can be calculated by the IPSn method when rc is greater than or equal to the longest lattice of the system. The IPSp method and the Wolf method require a longer cutoff radius than the longest lattice of the system to estimate interfacial properties.

UR - http://www.scopus.com/inward/record.url?scp=79957462395&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=79957462395&partnerID=8YFLogxK

U2 - 10.1063/1.3578473

DO - 10.1063/1.3578473

M3 - Article

VL - 134

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 17

M1 - 174112

ER -