计算化学公社

标题: 使用Gaussian和CP2K分别计算得到的同一化学键解离能，结果差异很大 [打印本页]

作者
Author: liyaqin 时间: 2026-4-23 22:19
标题: 使用Gaussian和CP2K分别计算得到的同一化学键解离能，结果差异很大
各位老师好，我想请教关于使用Gaussian和CP2K做势能面扫描的区别。
下图显示了用Gaussian和CP2K分别计算得到的化学键解离能，其中在Gaussian里基组赝势选择是b3lyp/6-31g，在CP2K里分别用b3lyp/6-31g和PBE/DZVP-MOLOPT-SR-GTH进行了计算，但是可以看出Gaussian计算得到的曲线就很合理和lammps计算得到的曲线重合程度很高，而CP2K计算得到的曲线就很奇怪，甚至曲线趋势也不合理。
想请问老师们这是什么原因呢？

作者
Author: sobereva 时间: 2026-4-24 09:53
确保CP2K做杂化泛函计算时没有犯基本常识性错误，进而不仅浪费耗时还导致SCF收敛到错误结果，要注意的下文都说了
CP2K做杂化泛函计算的关键要点和简单例子
http://sobereva.com/690（http://bbs.keinsci.com/thread-41087-1-1.html）
当前CP2K的PBE曲线也不平滑，确保SCF收敛足够充分、CUTOFF等影响数值精度的设置足够合理

描述要准确，“b3lyp/6-31g” 不叫 “基组赝势”，这叫计算级别

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