Polarizable TIP7P Water Model with Perturbation Charges Evaluated from ABEEM

• Published in: The journal of physical chemistry. B 2020-03, Vol.124 (12), p.2450-2464
• Main Authors: Zhao, Chong-Li, Zhao, Dong-Xia, Jiang, Qian-Ying, Zhang, Hai-Xia, Li, Shenmin, Yang, Zhong-Zhi
• Format: Article
• Language: English
• ISSN: 1520-6106; 1520-5207
• DOI: 10.1021/acs.jpcb.9b11775

A polarizable version of the rigid seven-site (TIP7P) water model with the atom-bond electronegativity equalization method (ABEEM) is proposed. The model uses direct polarization, where an isolated water monomer in the equilibrium geometry is assumed as a reference state and the polarization of the monomer arises from interacting with other molecules as a perturbation of the reference state. The charge on each site of the monomer splits into reference charge and perturbation charge. The perturbation charge arises only because of other reference charges. The interaction of the perturbation charge with other perturbation charges is replaced using polarization scaling to enhance the interaction of perturbation charge with the reference charges of the sites from other molecules. The perturbation charges are updated by evaluating explicit expressions once. This direct polarization is time-reversible because the charge update is independent of the charges in previous simulation steps. A Slater-type damping function moderates the short-range electrostatics to treat charge diffusion. The Ewald method corrects the long-range electrostatics both in the nuclei movement and in electronegativity equalization to diminish the size effect. The water model is parameterized by fitting the ab initio results of water clusters and the experimental results of water monomers and thermodynamic properties for liquid water. Owing to polarizability, the model performs better than the TIP7P model in terms of vaporization enthalpy, isothermal compressibility, and shear viscosity of the liquid phase. It performs better at the melting point of ice but slightly worse under critical conditions than the TIP7P model. Direct polarization has a low time complexity of O­(N) and is a good choice for ABEEM to improve its computational efficiency.