大体系主客体复合物相互作用能计算和BSSE校正问题

saaamon

近期在计算一个分子笼与富勒烯的相互作用能，分子笼较大，加上富勒烯总共超过300个原子数。首先用主客体复合物的单晶结构作为初猜，使用wB97M-V（def2-SVP）对主客体复合物的结构进行结构优化（orca计算耗时一两天，还算可以接受），想尝试使用def2-QZVP这样的4-zeta档次基组计算单点，并通过E_interaction = E_AB - E(A) - E(B)计算相互作用能，然而由于原子数过多，在orca上都很难算得动，因此想降低到def2-TZVP来，同时考虑BSSE校正（orca上使用gCP进行校正），来缩短计算的周期，但是得到的结果并不理想（不符合实验结果），并有以下疑问：
1. 主客体复合物结构优化时未使用溶剂模型，但是单点计算时使用了SMD模型，但是看到Sob的博文中提到“溶剂模型下考虑BSSE问题没有严格的办法”，是不是单点最好也在真空下计算？
2.主客体复合物结构优化后并进行单点计算可以得到E_AB，而对于单独的主客体单点E(A)和E(B)的计算存在疑问，应该使用如下哪种方法：a 通过把之前优化好的主客体复合物结构直接拆开，分别计算单点（加上gCP校正）得到E(A)和E(B)；b 通过把之前优化好的主客体复合物结构直接拆开，使用wB97M-V（def2-SVP）进行结构优化后，再分别计算单点（加上gCP校正）得到E(A)和E(B)？
3.主客体复合物结构优化时是否需要加上gCP校正？
请大神赐教！

sobereva

wB97M-V结合def2-TZVP结合gCP妥妥的，下文就是典型这么做的例子，可以作为参考文献引用证明合理性。结果与期望不符，不要盲目赖计算级别（胡乱赖计算级别是初学者最严重、常见问题之一），先搞清楚可能的误差来源都有哪些
全面揭示各种碳环与富勒烯之间独特的pi-pi相互作用！
http://sobereva.com/727（http://bbs.keinsci.com/thread-49376-1-1.html）

不要把wB97M-V写成wB97-MV

不推荐wB97M-V用于opt freq目的，优化慢，而且ORCA里都没有解析Hessian。B3LYP-D3(BJ)通常是又好又快的理想选择，离域性误差问题明显时改用wB97XD，看下文
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

def2-SVP优化弱相互作用复合物必须带gCP

问题2：仔细把下文看了便知
谈谈分子间结合能的构成以及分解分析思想
http://sobereva.com/733（http://bbs.keinsci.com/thread-50890-1-1.html）

saaamon

感谢sob老师赐教