模拟水表面张力，与实验值差别较大

五好先森 · 发表于 Post on 2023-9-6 20:24:11

本帖最后由五好先森于 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了。

thou · 发表于 Post on 2023-9-6 22:39:20

http://bbs.keinsci.com/thread-25479-1-1.html
不同水模型对表面张力的重现性有差异

五好先森 · 发表于 Post on 2023-9-7 08:34:16

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

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

sobereva · 发表于 Post on 2023-9-8 01:33:16

直接用OPC就完了，算各种属性都很好
跟程序本身毫无关系，而GROMACS效率远高于LAMMPS，毫无换LAMMPS的意义

五好先森 · 发表于 Post on 2023-9-8 08:44:51

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会导致界面张力偏差较大，能不能分情况使用两种模型啊？这样做又害怕会被质疑。

sobereva · 发表于 Post on 2023-9-8 14:37:27

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

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

五好先森 · 发表于 Post on 2023-9-8 19:36:24

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

明白明白，多谢社长。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 模拟水表面张力，与实验值差别较大

浏览过的版块