利用cp2k计算同一吸附模型得到的吸附能相对于VASP的计算值偏大该如何校正呢？

zbh · 发表于 Post on 2021-12-18 10:59:34

各位老师好，我想问一下利用cp2k计算同一吸附模型得到的吸附能相对于VASP的计算值偏大该如何校正呢？（比如用vasp算出来是-6，而cp2k的计算值为-10）

Aridea · 发表于 Post on 2021-12-18 18:46:07

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

ghifi37 · 发表于 Post on 2021-12-18 19:36:58

用CP2K做结构优化，再用VASP做单点是最方便的。

sobereva · 发表于 Post on 2021-12-18 22:40:14

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

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

卡开发发 · 发表于 Post on 2021-12-19 03:13:26

本帖最后由卡开发发于 2021-12-19 04:00 编辑

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

我个人经验而言，赝势的角度说，同样是平面波从Delta计算反映来看，VASP稍微新一些的或者那些带GW尾缀的赝势精度相对于GTH会高一些（https://molmod.ugent.be/deltacodesdft），此处的图是2007的数据。

不过就同样都使用GTH的情况应该也有相应的测试（https://www.cp2k.org/_media/events:2019_ghent:automation.pdf）：

Delta值的数据虽然也只是一个参考，我就是从客观数据讨论一下。

zbh · 发表于 Post on 2021-12-19 10:01:27

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

谢谢您的指导

zbh · 发表于 Post on 2021-12-19 10:07:32

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

谢谢老师的指导

zbh · 发表于 Post on 2021-12-19 10:08:10

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

感谢您的指导

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 利用cp2k计算同一吸附模型得到的吸附能相对于VASP的计算值偏大该如何校正呢？

评分 Rate

浏览过的版块