计算化学公社

标题: CP2K，利用cp2k计算同一吸附模型得到的吸附能相对于文章的计算值偏小该，该怎么办？ [打印本页]

作者
Author: tongxinsuo8s 时间: 2024-4-22 14:42
标题: CP2K，利用cp2k计算同一吸附模型得到的吸附能相对于文章的计算值偏小该，该怎么办？
本帖最后由 tongxinsuo8s 于 2024-4-22 14:45 编辑

各位老师好，我用CP2K计算单层Ti3C2O2吸附NO小分子，设置色散矫正为D3BJ，XY周期，MT，smear， DZVP-MOLOPT-SR-GTH，计算的吸附能是0.03eV，但是文章用castep计算得到的是0.4eV。我计算的NO分子与Ti3C2O2的O原子距离为2.4埃，文章是2.1埃。这是为什么呢？我有什么设置没有加上吗？

作者
Author: tongxinsuo8s 时间: 2024-4-22 16:56
嗷嗷，cp2k能量单位是hartree。。。我换算成eV以后，是0.64eV，和文献的0.41还是有差距。我通过vasp计算得到的吸附能是0.45eV，和文献差不多。

作者
Author: leeru 时间: 2024-4-22 17:28
我一般在做了收敛测试以后，我都会认为是文献的数据有问题。毕竟文献也是人写的，只要自己做了收敛性测试，就相信自己的计算吧

作者
Author: sobereva 时间: 2024-4-23 09:04
起码把DZVP-MOLOPT-SR-GTH基组先换了再说，这么小的基组只适合opt freq MD之类目的，算能量岂能用这么小的基组。尤其是算物理吸附吸附问题，如果还不考虑counterpoise或gCP，这基组BSSE问题大到姥姥家了。
顺带一提，很多人无端冤枉CP2K计算精度都是缺乏基组选用的最基本常识。

其它方面参考
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread-38911-1-1.html）

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