Multiwfn计算的Coulomb相互作用和自己编程计算的不一样

破晓时的曙光 · 发表于 Post on 2020-8-16 10:46:03

本帖最后由破晓时的曙光于 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代码。

sobereva · 发表于 Post on 2020-8-16 11:34:16

你没注意距离的单位转换，用Bohr算出来的才是Hartree

liuyuje714 · 发表于 Post on 2020-8-16 11:39:11

你的单位换算错了，实际算出来的E的单位是静电势能(e2/Å)，换成kj/mol需要乘以1389.35457644677，换成eV需要乘以14.3996454784562

破晓时的曙光 · 发表于 Post on 2020-8-16 14:20:34

的确是单位错了，谢谢各位老师

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

[辅助/分析程序] Multiwfn计算的Coulomb相互作用和自己编程计算的不一样

浏览过的版块