计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 波函数分析与Multiwfn (Wavefunction Analysis & Multiwfn) › 库仑相互作用计算求助

库仑相互作用计算求助

查看数: 1612 | 评论数: 7 | 收藏 Add to favorites 0

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

帖子模式

zhangxx

发布时间: 2023-6-3 17:41

正文摘要:

本帖最后由 zhangxx 于 2023-6-3 20:46 编辑各位老师我想问一下，我用文献的公式计算我自己的体系的Ecoul，用NBO电荷写的chg文件，用Multiwfn的7功能，-2功能计算的值为-1.10eV，这个值算的是文献的Ecoul值吗， ...

回复 Reply

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: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