计算化学公社

标题: Multiwfn计算的Coulomb相互作用和自己编程计算的不一样 [打印本页]

作者
Author: 破晓时的曙光 时间: 2020-8-16 10:46
标题: Multiwfn计算的Coulomb相互作用和自己编程计算的不一样
本帖最后由破晓时的曙光于 2020-8-16 14:20 编辑

我想计算两个分子间的库仑相互作用，自己编程计算了一下，但不确定对不对，发现Multiwfn手册4.A.9 计算excitonic coupling的公式与计算Coulomb相互作用的公式一致，只是所用电荷不同。因此，我用multiwfn也计算Coulomb相互作用，发现和自己编程计算结果不一样。仔细检查并未找到问题。
我的Multiwfn操作步骤如下：
1. 将包含分子坐标和电荷的test.chg文件拖入Multiwfn。
2. 输入 7 Population analysis and atomic charges
输入 -2 Calculate interaction energy between fragments based on atomic charges
输入1-6
输入7-18

最后Multiwfn计算得到的能量是 20.84 kJ/mol（0.216 eV），编程计算得到的是0.408 eV，这差的还是蛮多的。以下是我的test.chg文件和python代码。

作者
Author: sobereva 时间: 2020-8-16 11:34
你没注意距离的单位转换，用Bohr算出来的才是Hartree

作者
Author: liuyuje714 时间: 2020-8-16 11:39
你的单位换算错了，实际算出来的E的单位是静电势能(e2/Å)，换成kj/mol需要乘以1389.35457644677，换成eV需要乘以14.3996454784562

作者
Author: 破晓时的曙光 时间: 2020-8-16 14:20
的确是单位错了，谢谢各位老师

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