求助，体系的总电荷不为零

查看数: 1474 | 评论数: 5 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2023-5-30 16:03

在做NVT模拟时，体系的总电荷不为0，请问这个警告可以忽略吗？

回复 Reply

slxc920113 发表于 Post on 2023-6-1 08:38:43

wjx2022 发表于 2023-5-31 13:54
老师，我的itp文件中正负电荷只相差0.0002，我平均一下补充上去可以吗？

这种很可能是数值误差引起的，找个电荷较大的原子，手动改一改就行。

sobereva 发表于 Post on 2023-6-1 03:44:43

wjx2022 发表于 2023-5-31 13:54
老师，我的itp文件中正负电荷只相差0.0002，我平均一下补充上去可以吗？

语义不明

wjx2022 发表于 Post on 2023-5-31 13:54:15

sobereva 发表于 2023-5-31 03:43
明显不能忽略。
建议按下文计算1.2*CM5原子电荷，最适合OPLS-AA。Multiwfn输出的原子电荷位数很多，能确保 ...

老师，我的itp文件中正负电荷只相差0.0002，我平均一下补充上去可以吗？

sobereva 发表于 Post on 2023-5-31 03:43:44

明显不能忽略。
建议按下文计算1.2*CM5原子电荷，最适合OPLS-AA。Multiwfn输出的原子电荷位数很多，能确保加准确为0
计算适用于OPLS-AA力场做模拟的1.2*CM5原子电荷的懒人脚本
http://sobereva.com/585（http://bbs.keinsci.com/thread-21462-1-1.html）
ORCA结合Multiwfn计算RESP、RESP2和1.2*CM5原子电荷的懒人脚本
http://sobereva.com/637（http://bbs.keinsci.com/thread-28178-1-1.html）

Frozen-Penguin 发表于 Post on 2023-5-31 00:03:13

小数点后两位不为0，可能是力场参数中的电荷量不对，不能忽略

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