使用CP2K计算Fe2O3表面吸附结合能与文献差距较大

喝酸奶的巧乐兹 · 发表于 Post on 2025-12-11 23:03:57

各位老师好！我使用CP2K计算NH3在Fe2O3 001表面的吸附能。计算过程如下：
1.构建模型大小为15A-15A-25A，使用OT算法，每一步均开启D3(BJ)色散矫正，cutoff为1050，55。结构如下。
2.使用DZVP-MOLOPT-SR-GTH基组，OT、D3(BJ)、DFT+U=2eV，进行Fe2O3 001表面的几何优化、Fe2O3-NH3的几何优化
3.使用TZVP-MOLOPT(-SR)-GTH基组，OT、D3(BJ)、关闭DFT+U、分别计算Fe2O3、NH3、Fe2O3-NH3的能量，再利用能量相减得到结合能E1=-1.4076eV
4.使用TZVP-MOLOPT(-SR)-GTH基组，OT、D3(BJ)、关闭DFT+U、使用BSSE计算Fe2O3-NH3的结合能，得到结合能E2=-1.3463eV

多篇文献使用CASTEP软件计算得到吸附能在-1.25eV左右，而我使用CP2K不管是否使用BSSE，结合能均比文献大0.1eV以上，请问我的计算方法或参数设置有问题吗？该如何改正呢？
不使用BSSE的三个能量计算输入文件和使用BSSE的能量输入文件均已附上，希望各位老师不吝赐教！谢谢大家！

sobereva · 发表于 Post on 2025-12-12 01:52:40

CASTEP的结果又不是金标准，干嘛非要和它的结果一致。本来两个程序用的赝势、基函数都十分不一样。-1.3463eV和-1.25 eV哪算得上相差较大。
并且文献里具体用的什么计算级别要交代清楚，否则别人都不知道在理论方法层面是否有可比性。

还要注意相互作用能和吸附能不是一码事，下文说了。counterpoise任务给你的只是BSSE校正后的相互作用能
谈谈分子间结合能的构成以及分解分析思想
http://sobereva.com/733（http://bbs.keinsci.com/thread-50890-1-1.html）

喝酸奶的巧乐兹 · 发表于 Post on 2025-12-12 02:11:38

sobereva 发表于 2025-12-12 01:52
CASTEP的结果又不是金标准，干嘛非要和它的结果一致。本来两个程序用的赝势、基函数都十分不一样。-1.3463e ...

谢谢sob老师的耐心回复！
文献使用了CASTEP，超软赝势，截断能400eV，布里渊区采样使用2×2×1的Monkhorst-Pack网格，有0.2eV的费米展宽，Fe 3d电子的U值取4.0 eV。不知道我的设置与文献设置在是否具有可比性呢？

sobereva · 发表于 Post on 2025-12-12 02:23:27

喝酸奶的巧乐兹发表于 2025-12-12 02:11
谢谢sob老师的耐心回复！
文献使用了CASTEP，超软赝势，截断能400eV，布里渊区采样使用2×2×1的Monkhor ...

用了什么泛函，带没带色散校正都没说
倘若理论方法是精确一样的、模型是一样的、计算吸附能的流程是一样的，这俩数值方法截然不同的程序结果相差这么点完全是在正常范围内

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

[CP2K] 使用CP2K计算Fe2O3表面吸附结合能与文献差距较大