计算化学公社

标题: 利用cp2k计算同一吸附模型得到的吸附能相对于VASP的计算值偏大该如何校正呢？ [打印本页]

作者
Author: zbh 时间: 2021-12-18 10:59
标题: 利用cp2k计算同一吸附模型得到的吸附能相对于VASP的计算值偏大该如何校正呢？
各位老师好，我想问一下利用cp2k计算同一吸附模型得到的吸附能相对于VASP的计算值偏大该如何校正呢？（比如用vasp算出来是-6，而cp2k的计算值为-10）

作者
Author: Aridea 时间: 2021-12-18 18:46
首先保证能量计算精度足够高，让所有参数收敛了；其次考虑BSSE可以再论坛搜一下怎么操作

作者
Author: ghifi37 时间: 2021-12-18 19:36
用CP2K做结构优化，再用VASP做单点是最方便的。

作者
Author: sobereva 时间: 2021-12-18 22:40
没有任何前提下，本来就不存在VASP比CP2K算得更准一说，至于要将CP2K结果向VASP的结果做校正完全是在不当地贬损CP2K。
至少先说清楚俩程序都具体怎么算的。而且什么类型吸附说清楚，以判断CP2K考虑BSSE的必要性

仔细看
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

作者
Author: 卡开发发 时间: 2021-12-19 03:13
本帖最后由卡开发发于 2021-12-19 04:00 编辑

其实可以把体系和计算细节拿出来看看，便于排查。-6和-10这样的偏差也得看，如果是eV这样的单位可能要怀疑是不是算错了，要是kcal/mol的单位，那可能LCAO和PW就是会有这样的偏差，也在容忍的范围内，有时候提升基组质量或修正BSSE也许能够改善类似问题。另外赝势本身的差异也会造成两者结果差异，这点较为难排查。

我个人经验而言，赝势的角度说，同样是平面波从Delta计算反映来看，VASP稍微新一些的或者那些带GW尾缀的赝势精度相对于GTH会高一些（https://molmod.ugent.be/deltacodesdft），此处的图是2007的数据。
(, 下载次数 Times of downloads: 20)
不过就同样都使用GTH的情况应该也有相应的测试（https://www.cp2k.org/_media/events:2019_ghent:automation.pdf）：
(, 下载次数 Times of downloads: 29)
Delta值的数据虽然也只是一个参考，我就是从客观数据讨论一下。

作者
Author: zbh 时间: 2021-12-19 10:01

Aridea 发表于 2021-12-18 18:46
首先保证能量计算精度足够高，让所有参数收敛了；其次考虑BSSE可以再论坛搜一下怎么操作

谢谢您的指导

作者
Author: zbh 时间: 2021-12-19 10:07

sobereva 发表于 2021-12-18 22:40
没有任何前提下，本来就不存在VASP比CP2K算得更准一说，至于要将CP2K结果向VASP的结果做校正完全是在不当地 ...

谢谢老师的指导

作者
Author: zbh 时间: 2021-12-19 10:08

卡开发发发表于 2021-12-19 03:13
其实可以把体系和计算细节拿出来看看，便于排查。-6和-10这样的偏差也得看，如果是eV这样的单位可能要怀疑 ...

感谢您的指导

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