计算化学公社

标题: Multifwfn基于力场能量分解和counterpoise计算分子间弱相互作用的差异问题 [打印本页]

作者
Author: lily1970 时间: 2020-7-29 15:47
标题: Multifwfn基于力场能量分解和counterpoise计算分子间弱相互作用的差异问题
老师好，我再分别利用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

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

作者
Author: 向往阳光 时间: 2020-7-29 17:24
1.用不同方法计算的能量的绝对值一般不能直接比较。
2.当然是高斯计算的准确，Multiwfn只是基于分子力场分解，有时候参数会不合理，那样自然会定性错误。

作者
Author: lily1970 时间: 2020-7-29 20:00

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

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

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

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