### 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, r_{c}. However, the estimation of the electrostatic potential profile and dipole orientational function is poor. The Wolf method cannot estimate electrostatic configuration at r_{c} L_{z}2 (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 r_{c} L _{z}3. 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 r_{c} L _{z}3. In spite of this shortcoming, the IPSn method gives the most accurate result in estimating the surface tension at r_{c} L _{z}2. 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 r_{c} 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 language | English |
---|---|

Article number | 174112 |

Journal | Journal of Chemical Physics |

Volume | 134 |

Issue number | 17 |

DOIs | |

Publication status | Published - 2011 May 7 |

### Fingerprint

### ASJC Scopus subject areas

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

### Cite this

*Journal of Chemical Physics*,

*134*(17), [174112]. https://doi.org/10.1063/1.3578473

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

Research output: Contribution to journal › Article

*Journal of Chemical Physics*, vol. 134, no. 17, 174112. https://doi.org/10.1063/1.3578473

}

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

C2 - 21548678

AN - SCOPUS:79957462395

VL - 134

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 17

M1 - 174112

ER -