关于BSSE能量矫正的详细处理过程。

molx79

本人才学疏浅，恳请大家赐教。目前做的体系一直是弱相互作用。我验证一片文献，目的是检验我理解的是否正确。体系：SO2...HCHO。
文献内容：mp2/aug-cc-pvdz条件下opt ，freq，为了得到更为准确的能量，在此优化的基础上用CCSD（T）/aug-cc-pvtz计算单点。分别得到Eb和Eint
疑问：1、mp2条件下处理的Eint和文献结果相差不多。但是CCSD（T)条件下处理的Eint就会差很多。
         我的处理过程如下：（1）讨论Eint时，复合物opt是加上关键词counter=N。
                                          （2）mp2/DZ  opt  freq  counter=N  计算条件下得到的复合物的结构用CCSD（T）/TZ计算单点的EAB
                                           （3）mp2/DZ  opt  freq  counter=N  计算条件下得到的复合物的结构的两个片段分别用CCSD（T）/TZ计算单点EA*和EB*
                                            （4）Eint=EAB-EA*-EB*，如果考虑上BSSE，则把mp2/DZ  opt  freq  counter=N  计算条件下得到的BSSE加进去。
                                         和文献的结果严重不符。
请大家帮帮看看哪有问题？不胜感激。