计算化学公社

标题: 求助：CP2K计算矿物表面吸附能结果过大 [打印本页]

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

作者
Author: spinel 时间: 2026-1-14 15:40
目前的任务是结构优化，要计算吸附能，起码要在更昂贵的基组条件下计算单点能，DZVP用来计算能量不合适

作者
Author: sobereva 时间: 2026-1-15 02:06
看清楚版块再发帖，CP2K根本不是量子化学程序，必须发到第一性原理版块的CP2K分类，以后注意！

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

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

作者
Author: yuzc 时间: yesterday 16:34
没加BSSE

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