计算化学公社

标题: 库仑相互作用计算求助 [打印本页]

作者
Author: zhangxx 时间: 2023-6-3 17:41
标题: 库仑相互作用计算求助
本帖最后由 zhangxx 于 2023-6-3 20:46 编辑

各位老师我想问一下，我用文献的公式计算我自己的体系的Ecoul，用NBO电荷写的chg文件，用Multiwfn的7功能，-2功能计算的值为-1.10eV，这个值算的是文献的Ecoul值吗，我看文献算的体系的值才-0.3eV，我的值好像有点太大了

作者
Author: wzkchem5 时间: 2023-6-3 19:12
贴文献原文

作者
Author: zhangxx 时间: 2023-6-3 20:17

wzkchem5 发表于 2023-6-3 19:12
贴文献原文

DOI: 10.1021/acs.jpcc.6b09927

作者
Author: zhangxx 时间: 2023-6-3 20:47
这是文献原文的图片

作者
Author: wzkchem5 时间: 2023-6-4 00:49

zhangxx 发表于 2023-6-3 13:47
这是文献原文的图片

文献用的是Hirshfeld电荷，你用的是NPA电荷（注意那个不叫NBO电荷，尽管是NBO程序算出来的），结果当然不可比。
再者，你设置介电常数了吗（文献里设的4，Multiwfn这个功能我不清楚但是默认很可能是1）

作者
Author: sobereva 时间: 2023-6-4 03:07
NPA电荷的数值整体比Hirshfeld大得多，《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）里的对比充分体现了，所以基于NPA电荷算的静电相互作用能肯定比Hirshfeld电荷也大很多。

对于基于原子电荷计算静电相互作用能，最适合基于拟合静电势电荷，因此下文的方法也都是这么用的
使用Multiwfn做基于分子力场的能量分解分析
http://sobereva.com/442（http://bbs.keinsci.com/thread-10907-1-1.html）

Multiwfn的那个功能是当真空环境计算的，注意介电常数问题。

作者
Author: zhangxx 时间: 2023-6-4 14:52

sobereva 发表于 2023-6-4 03:07
NPA电荷的数值整体比Hirshfeld大得多，《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstra ...

那请问老师，我的介电常数应该怎么处理

作者
Author: sobereva 时间: 2023-6-6 09:08

zhangxx 发表于 2023-6-4 14:52
那请问老师，我的介电常数应该怎么处理

看实际研究的环境，诸如真空条件相对介电常数就是1

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