求助：CP2K计算矿物表面吸附能结果过大

不要椰奶

本人使用cp2k计算小分子在矿物表面吸附能时，计算两组结果分别为Eads1=-3.817825367 Hartree，Eads2=-3.861938923 Hartree，与一般化学吸附结果相差太多，不知是否是inp文件设置有误，请大神们指点，inp文件中有何改进之处，两组文件均在相同精度条件下计算，计算级别为PBE/DZVP-MOLOPT-SR-GTH, 使用OT算法，开启DFT-D3(BJ)色散校正。

spinel

目前的任务是结构优化，要计算吸附能，起码要在更昂贵的基组条件下计算单点能，DZVP用来计算能量不合适

sobereva

看清楚版块再发帖，CP2K根本不是量子化学程序，必须发到第一性原理版块的CP2K分类，以后注意！

吸附能有几Hartree明显严重算错了。光给复合物的输入文件没给各个片段的没法完整判断情况。并且注意检查SCF收敛等细节避免存在硬伤。
DZVP-MOLOPT-SR-GTH这么小、BSSE问题那么大的基组算吸附能必定导致明显高估吸附能（虽然不会导致你当前结果那么离谱，结果离谱另有原因）。这基组只适合优化，单点至少也应该用TZVP-MOLOPT-GTH档次

PS：第5届起的北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里专门有一节“吸附能的计算”，详细讲了吸附能计算需要注意的各种细节，并给了十分完整的例子演示怎么正确计算并得到合理结果，若学过一遍肯定不会犯任何错误。

yuzc

没加BSSE