计算化学公社

标题: 模拟水表面张力，与实验值差别较大 [打印本页]

作者
Author: 五好先森 时间: 2023-9-6 20:24
标题: 模拟水表面张力，与实验值差别较大
本帖最后由五好先森于 2023-9-6 20:31 编辑

我用Gromcas模拟298K，1atm下水（spc/e模型）的表面张力，为60 mN/m，与社长的文章（https://www.whxb.pku.edu.cn/CN/abstract/abstract29166.shtml）中结果相符，但是都明显小于实际测量值（~70 mN/m）。

但是有些文献中用LAMMPS，能得到非常接近实际值的水表面张力，如：https://doi.org/10.1021/jp048773n，该文献用的是F3C水模型。

请问，利用Gromacs能获得更准确的表面张力结果吗？是不是我只能改学LAMMPS了。

作者
Author: thou 时间: 2023-9-6 22:39
http://bbs.keinsci.com/thread-25479-1-1.html
不同水模型对表面张力的重现性有差异

作者
Author: 五好先森 时间: 2023-9-7 08:34

thou 发表于 2023-9-6 22:39
http://bbs.keinsci.com/thread-25479-1-1.html
不同水模型对表面张力的重现性有差异

多谢提点。我主要关注油水界面张力，看来得多尝试几种水模型了。

作者
Author: sobereva 时间: 2023-9-8 01:33
直接用OPC就完了，算各种属性都很好
跟程序本身毫无关系，而GROMACS效率远高于LAMMPS，毫无换LAMMPS的意义

作者
Author: 五好先森 时间: 2023-9-8 08:44

sobereva 发表于 2023-9-8 01:33
直接用OPC就完了，算各种属性都很好
跟程序本身毫无关系，而GROMACS效率远高于LAMMPS，毫无换LAMMPS的意义

社长好，我换了OPC水模型计算水的表面张力，298 K、1 atm条件下结果为73.49 mN/m，基本符合实验值（~70 mN/m），比SPC/E更准确。但是用OPC水模型计算的298 K、1 atm条件下水-正癸烷界面张力为58.00 mN/m，与实验值（51.98 mN/m）的偏差大于SPC/E水模型计算结果（52.21 mN/m）。

这让我感到有些纠结，当既要关注油-水界面张力又要关注水表面张力的情况下，担心统一都用OPC会导致界面张力偏差较大，能不能分情况使用两种模型啊？这样做又害怕会被质疑。

作者
Author: sobereva 时间: 2023-9-8 14:37

五好先森发表于 2023-9-8 08:44
社长好，我换了OPC水模型计算水的表面张力，298 K、1 atm条件下结果为73.49 mN/m，基本符合实验值（~70 m ...

如果确信误差就是来自于水模型而不是其它因素，也可以尝试OPC3或TIP4P-FB，看是否有更均衡的描述

作者
Author: 五好先森 时间: 2023-9-8 19:36

sobereva 发表于 2023-9-8 14:37
如果确信误差就是来自于水模型而不是其它因素，也可以尝试OPC3或TIP4P-FB，看是否有更均衡的描述

明白明白，多谢社长。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)