请问VASP能否计算带电表面对中性气体物理吸附过程中的范德华力呢

yidawang

各位朋友大家好，我们组最近在研究带电石墨烯极板对甲烷等非极性气体分子的物理吸附过程，希望通过计算来获得一些参考数据。通过阅读VASP手册，我们知道VASP似乎不好处理电中性表面对离子的吸附。那么请问一下带电表面对中性气体的吸附，计算结果有可信度吗？

具体的计算方法如下：
构建一个在z方向上20Å以上的真空层，不进行结构优化，而是直接计算在不同位置处 （带电系统总能量）−（单独的带电极板能量）−（单独的电中性气体能量），并以这个结果作为范德华力的能量。通过能带图我们可以发现，在费米能级附近是石墨烯的狄拉克锥，因此少量的电荷（0.05e）总是位于石墨烯的轨道上，而气体分子仍然是电中性的。

计算结果遇到的主要问题是，在添加不同大小的电荷时，范德华力的能量在非常远的距离仍然有一定的区别。结果如图所示，图示是 带电情况下范德华力与不带电情况下范德华力的区别。



请问一下这个计算结果可信吗？是否需要添加例如IDIPOL等参数来修正呢？

我已经阅读了论坛内有关表面吸附的帖子，但是仍然没有解决我的问题，例如以下几篇帖子。

http://bbs.keinsci.com/thread-7016-1-1.html

http://bbs.keinsci.com/thread-14426-1-1.html

http://bbs.keinsci.com/thread-9422-1-1.html

希望各位能够不吝赐教，非常感谢！

卡开发发

不可信，相同的真空层情况可以比较相对能量。石墨烯比较难处理，如果是一般表面体系，我觉得可能的做法是，可以构造一个上下对称的表面来进行处理，然后用VASPSOL对FERMI_SHIFT进行修正。

yidawang

卡开发发 发表于 2021-1-9 22:11
不可信，相同的真空层情况可以比较相对能量。石墨烯比较难处理，如果是一般表面体系，我觉得可能的做法是， ...

非常感觉卡开发发前辈！

您这边提到的相同真空层可以比较相对能量，现在我想出来的一个方案是，构建了一个100Å的真空层，中间只有一个甲烷气体分子，我直接设定50Å处为范德华力的零点来修正。请问对于同一个NELECT的系统，比较甲烷气体在不同位置时的相对能量，这种比较方法可行吗？

再次感谢前辈的指导！

卡开发发

yidawang 发表于 2021-1-10 03:57
非常感觉卡开发发前辈！

您这边提到的相同真空层可以比较相对能量，现在我想出来的一个方案是，构建了 ...

应该仅可以讨论趋势。

yidawang

卡开发发 发表于 2021-1-10 09:15
应该仅可以讨论趋势。

好的！感谢前辈！祝新年快乐，工作顺利！