计算化学公社

标题: 求助，体系的总电荷不为零 [打印本页]

作者
Author: wjx2022 时间: 2023-5-30 16:03
标题: 求助，体系的总电荷不为零
在做NVT模拟时，体系的总电荷不为0，请问这个警告可以忽略吗？

作者
Author: Frozen-Penguin 时间: 2023-5-31 00:03
小数点后两位不为0，可能是力场参数中的电荷量不对，不能忽略

作者
Author: sobereva 时间: 2023-5-31 03:43
明显不能忽略。
建议按下文计算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）

作者
Author: wjx2022 时间: 2023-5-31 13:54

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

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

作者
Author: sobereva 时间: 2023-6-1 03:44

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

语义不明

作者
Author: slxc920113 时间: 2023-6-1 08:38

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

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

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