Science Institute of the University of Iceland, VR-III, 107, Reykjavík, Iceland.
Phys Chem Chem Phys. 2013 Oct 21;15(39):16542-56. doi: 10.1039/c3cp52097h. Epub 2013 Aug 15.
A transferable potential energy function for describing the interaction between water molecules is presented. The electrostatic interaction is described rigorously using a multipole expansion. Only one expansion center is used per molecule to avoid the introduction of monopoles. This single center approach turns out to converge and give close agreement with ab initio calculations when carried out up to and including the hexadecapole. Both dipole and quadrupole polarizability are included. All parameters in the electrostatic interaction as well as the dispersion interaction are taken from ab initio calculations or experimental measurements of a single water molecule. The repulsive part of the interaction is parametrized to fit ab initio calculations of small water clusters and experimental measurements of ice Ih. The parametrized potential function was then used to simulate liquid water and the results agree well with experiment, even better than simulations using some of the point charge potentials fitted to liquid water. The evaluation of the new interaction potential for condensed phases is fast because point charges are not present and the interaction can, to a good approximation, be truncated at a finite range.
本文提出了一个可用于描述水分子间相互作用的转移势能函数。静电相互作用通过使用多极展开式进行严格描述。每个分子仅使用一个展开中心,以避免引入单极。事实证明,这种单中心方法可以收敛,并在进行到包括十六极为止时与从头算计算结果吻合得很好。同时考虑了偶极子和四极子极化率。静电相互作用以及色散相互作用中的所有参数均来自于单个水分子的从头算计算或实验测量。相互作用的排斥部分被参数化以拟合小分子水团簇的从头算计算和冰 Ih 的实验测量结果。然后使用参数化的势能函数来模拟液态水,结果与实验吻合得很好,甚至优于使用一些拟合液态水的点电荷势的模拟结果。由于不存在点电荷,因此评估凝聚相的新相互作用势能的速度很快,并且可以很好地将相互作用截断到有限的范围内。