库仑相互作用计算求助

zhangxx · 发表于 Post on 2023-6-3 17:41:20

本帖最后由 zhangxx 于 2023-6-3 20:46 编辑

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

sobereva · 发表于 Post on 2023-6-6 09:08:07

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

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

zhangxx · 发表于 Post on 2023-6-4 14:52:45

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

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

sobereva · 发表于 Post on 2023-6-4 03:07:56

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的那个功能是当真空环境计算的，注意介电常数问题。

wzkchem5 · 发表于 Post on 2023-6-4 00:49:36

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

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

zhangxx · 发表于 Post on 2023-6-3 20:47:35

这是文献原文的图片

zhangxx · 发表于 Post on 2023-6-3 20:17:31

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

DOI: 10.1021/acs.jpcc.6b09927

wzkchem5 · 发表于 Post on 2023-6-3 19:12:10

贴文献原文

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

[Multiwfn使用咨询] 库仑相互作用计算求助