计算化学公社

标题: VASP计算H2O2的吸附能如何正确处理？ [打印本页]

作者
Author: XIEH 时间: 2025-6-14 17:02
标题: VASP计算H2O2的吸附能如何正确处理？
今天看到了sob老师关于“基组重叠误差（Basis Set Superposition Error, BSSE）”的博文http://sobereva.com/46，想到了之前使用VASP计算H2O2吸附能的疑惑：
①：E(*H2O2) = E(Slab*H2O2) - E(Slab) - E(H2O2)
之前按照上述式子处理得到的H2O2吸附能约为1 eV
后来尝试优化出一个box，其中包含slab与H2O2，但是H2O2没有吸附在slab上（至少间隔3埃，可能有物理吸附）
按照 ②：E(*H2O2) = E(Slab+H2O2) - E(Slab) - E(H2O2) 得到吸附能为0.5 eV

这0.5 eV的差值是来自BSSE还是由于剔除了物理吸附还是两者因素都有呢？

作者
Author: liuxp 时间: 2025-6-14 17:05
vasp 应该不存在BSSE吧，他用的是平面波

作者
Author: sobereva 时间: 2025-6-15 05:33
原子中心基函数才有BSSE的事，第一性原理程序里诸如CP2K（且不是平面波的SIRIUS模块）才需要考虑，QE、VASP等平面波程序都没BSSE的事。

作者
Author: XIEH 时间: 2025-6-16 15:52
“②：E(*H2O2) = E(Slab+H2O2) - E(Slab) - E(H2O2) ”写错了，应该是 E(*H2O2) = E(Slab*H2O2) - E(Slab+H2O2)

作者
Author: XIEH 时间: 2025-6-16 15:53

liuxp 发表于 2025-6-14 17:05
vasp 应该不存在BSSE吧，他用的是平面波

是这样啊，理论功底不行啊我

作者
Author: XIEH 时间: 2025-6-16 15:53

sobereva 发表于 2025-6-15 05:33
原子中心基函数才有BSSE的事，第一性原理程序里诸如CP2K（且不是平面波的SIRIUS模块）才需要考虑，QE、VASP ...

明白了，感谢社长的指点

作者
Author: liuxp 时间: 2025-6-16 21:35

XIEH 发表于 2025-6-16 15:52
“②：E(*H2O2) = E(Slab+H2O2) - E(Slab) - E(H2O2) ”写错了，应该是 E(*H2O2) = E(Slab*H2O2) - E(Slab+ ...

我觉得2不算吸附能

作者
Author: XIEH 时间: 2025-6-17 15:39

liuxp 发表于 2025-6-16 21:35
我觉得2不算吸附能

所以您认为就是应该拿box1（吸附物吸附在slab表面）-box2（slab）-box3（吸附物）来计算吸附能吗？
我做了一年计算，就是有种感觉：两部分反应物放在一个box里的能量就是比分别放在两个box里的能量之和要低，即使两部分隔得很远（2.7埃），是由于范德华作用导致的吗？那这一部分是否应该算入吸附能呢？

作者
Author: liuxp 时间: 2025-6-17 18:23
本帖最后由 liuxp 于 2025-6-17 18:25 编辑

XIEH 发表于 2025-6-17 15:39
所以您认为就是应该拿box1（吸附物吸附在slab表面）-box2（slab）-box3（吸附物）来计算吸附能吗？
我做 ...

我看那些视频吸附能是不能放到一个盒子里去的，我看sob老师的这个帖子http://sobereva.com/733，吸附能本身就包含吸附以后，slab和小分子吸附前后，由于形变引起的能量变化以及吸附两个相互作用能量之和。我个人觉得放到一个盒子里，哪怕有距离，也有一部分吸引作用，这才是导致两种情况不一样的结果，我个人看法，不一定对。

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