计算化学公社

标题: 求助RESP电荷分析和Mulliken电荷差异的问题 [打印本页]

作者
Author: zhangzh 时间: 2020-6-25 01:04
标题: 求助RESP电荷分析和Mulliken电荷差异的问题
我在计算乙烷分子的电荷分布的时候，发现使用高斯自带的电荷分布选项查看Mulliken电荷分布的结果如下
(, 下载次数 Times of downloads: 31)
但我在使用Multiwfn计算其RESP电荷时，发现RESP电荷的数值和上述值存在很大的差异，请问在什么体系中应该采用Mulliken电荷，什么体系中应该采用RESP电荷呢？
(, 下载次数 Times of downloads: 19)

作者
Author: zhangzh 时间: 2020-6-25 01:07
(, 下载次数 Times of downloads: 26)
抱歉一楼的RESP电荷有错误，这张图是正确操作后生成的，定性没有问题但是定量上差异较大

作者
Author: hebrewsnabla 时间: 2020-6-25 01:14
本帖最后由 hebrewsnabla 于 2020-6-25 01:16 编辑

先看这个
《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）
RESP拟合静电势电荷的原理以及在Multiwfn中的计算 http://bbs.keinsci.com/forum.php ... =10880&fromuid=5840

采用哪种电荷不是取决于体系，而是取决于计算目的。

作者
Author: zhangzh 时间: 2020-6-25 01:22

hebrewsnabla 发表于 2020-6-25 01:14
先看这个
《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）
...

谢谢，第二篇文章我看了，第一篇文章好像一直打不开，还有我发现文章范例中的RMSE和RRMSE值都很小，但是我的RRMSE值都在0.9以上，我的步骤是完全按照例子做的，这是说明拟合结果不好吗？应该如何降低RRMSE值呢

作者
Author: sobereva 时间: 2020-6-25 22:31

zhangzh 发表于 2020-6-25 01:22
谢谢，第二篇文章我看了，第一篇文章好像一直打不开，还有我发现文章范例中的RMSE和RRMSE值都很小，但是 ...

第一篇能正常打开

对当前体系这两种电荷根本就没多大差异，你看绝对差异才多少，都远没达到0.1
当前你的计算结果并没问题。RRMSE大是因为本身分子范德华表面附近静电势数值的绝对值就非常小，所以显得好像相对差大似的，但从绝对误差上来看根本没多大，正如RMSE所体现的。

作者
Author: xinxinup 时间: 2020-7-12 10:26
楼主的RESP电荷是写错了吗？怎么比Mulliken电荷少了个数量级？

作者
Author: sobereva 时间: 2020-7-13 11:19

xinxinup 发表于 2020-7-12 10:26
楼主的RESP电荷是写错了吗？怎么比Mulliken电荷少了个数量级？

没错
本来烃类的原子电荷数量级就非常小，绝对差异根本就是微乎其微。有的力场还直接把烃类所有原子电荷当成0

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