请教BSSE校正如何在隐式溶剂使用

dinghongming

我想算两个分子间的相互作用（主要是pi-pi堆积作用），参考sobereva老师相关博文http://sobereva.com/46，需要加入BSSE校正。在真空中，我在gaussian16中用如下关键词#  b3lyp/6-311+g(d,p) counterpoise=2 em=gd3bj  可以正常运行。
考虑到我的两个分子实际是在氯仿溶液中，所以我加入了隐式模型，#  b3lyp/6-311+g(d,p)  scrf=(solvent=chloroform) counterpoise=2 em=gd3bj
不过这样一开始就报错，把“counterpoise=2”去掉就可以跑了，不过这样只能算复合体的能量了，而且更关键的是无法加入BSSE校正了，不知该如何解决？谢谢！

ShuangfeiZhu

仔细看博文：谈谈BSSE校正与Gaussian对它的处理 http://sobereva.com/46：
“3 溶剂模型下考虑BSSE问题没有严格的办法，一般做法是在气相下做counterpoise计算得到E_BSSE，然后加到溶剂模型下以常规方式E(AB)-E(A)-E(B)计算的相互作用能上。”

dinghongming

ShuangfeiZhu 发表于 2020-1-15 12:56
仔细看博文：谈谈BSSE校正与Gaussian对它的处理 http://sobereva.com/46：
“3 溶剂模型下考虑BSSE问题没 ...

谢谢！