Multifwfn基于力场能量分解和counterpoise计算分子间弱相互作用的差异问题

lily1970

老师好，我再分别利用Multifwn和Counterpoise计算分子间弱相互作用时发现，两种方法得出的结果差距较大，
以下是我的计算步骤，烦请老师指点下哪里出了问题：

背景：计算多氯联苯二聚体中多氯联苯单体之间的相互作用

1、使用Multiwfn基于分子力场的能量分解分析，利用B3LYP-D3(BJ)/6-311G**优化，按照博文中http://sobereva.com/442计算水二聚体的例子进行计算，
得出的结果为：
Interaction energy components between all fragments:
                         Electrostatic   Repulsion   Dispersion     Total
Frag   1 -- Frag   2:         3.80       155.51      -136.96        22.35 kJ/mol
2、利用Gaussian进行计算，并进行BSSE校正
（1）优化二聚体的结构，关键字为：# opt freq b3lyp/6-311g(d,p) em=gd3bj
（2）利用优化后的结构进行Counterpoise任务，关键字为：# counterpoise=2 b3lyp/6-311+g(2d,p) em=gd3bj
得出的结果为：
Counterpoise corrected energy =   -4603.976740867549
                   BSSE energy =       0.005356537390
               sum of monomers =   -4603.953138876017
           complexation energy =     -18.17 kcal/mole (raw)
           complexation energy =     -14.81 kcal/mole (corrected) = -61.99 kJ/mol


想问一下老师，我的计算过程是哪里出了问题呢，还有以上两种方法哪种的精确度较高呢？谢谢老师

向往阳光

1.用不同方法计算的能量的绝对值一般不能直接比较。
2.当然是高斯计算的准确，Multiwfn只是基于分子力场分解，有时候参数会不合理，那样自然会定性错误。

lily1970

向往阳光 发表于 2020-7-29 17:24
1.用不同方法计算的能量的绝对值一般不能直接比较。
2.当然是高斯计算的准确，Multiwfn只是基于分子力场分 ...

谢谢老师~那么我的操作过程是没问题的对吧

sobereva

没有你的做EDA-FF分解的具体细节，也没有结构图，没法回答是什么原因导致的显著差异
仔细检查原子电荷、原子类型的设置合理性。
既有可能EDA-FF操作或设置不当，也有可能当前力场模型、静电作用考虑方式对当前体系不理想。
最好先把博文里的碱基对的例子完整重复一遍。