功函数计算结果与开尔文探针结果不符

qwe1832518773

我用VASP计算了一下石墨烯的功函（没有真空层），然后用vaspkit-426-3处理，得到z方向曲线如下
功函数的数值都要-30eV以上了，但是开尔文探针结果是4.53eV，一个是负值一个是正值，而且数值查这么多，这个算的对吗？
--------------------------------------------------------------------
今天下午又新算了一个又真空层的，曲线如下：

看了一下真空层那边的电势是4.65~5.14左右波动，然后在这个文件夹里用grep E-fermi OUTCAR查到 E-fermi :  -0.3687，根据功函数的定义是电子从费米能级逸散到真空能级需要的最小能量，所以我需要用真空能级减去费米能级等于5.0187~5.5087，这个就是计算得到的功函数，我这个理解对吗？

那么我用vaspkit得到的这个z方向曲线究竟是什么意义呢？

卡开发发

模型发来看看，初步怀疑真空层不够或者静电势平均的方向搞错了。

qwe1832518773

卡开发发 发表于 2021-1-16 13:18
模型发来看看，初步怀疑真空层不够或者静电势平均的方向搞错了。

卡老师，我更新了提问，模型也加在附件里了，请您再帮我解答一下疑惑，十分感谢！

卡开发发

1、块体模型没有真空，所以没有功函数（模型意义上的），所以graphite那个不要去折腾了。
2、理解基本正确。指的是电子从Fermi能级电离到真空静止能级所需要的能量（与非周期边界下的Koopmans定理是对应的），但是周期边界下没有严格意义上的真空静止能级，我们一般取远离表面且电势平稳的区域（如你第二张图大致横坐标25～35的区域）作为真空静止能级。严格来说，非周期边界条件下，真空静止能级按照电动力学的定义取两个带电粒子无穷远的势能为零点，周期性边界条件下，因为静电作用需要特殊处理导致能级被整体平移，所以只能通过这种方式来进行定义。

卡开发发

qwe1832518773 发表于 2021-1-16 20:58
卡老师，我更新了提问，模型也加在附件里了，请您再帮我解答一下疑惑，十分感谢！

1、块体模型没有真空，所以没有功函数（模型意义上的），所以graphite那个不要去折腾了。
2、理解基本正确。指的是电子从Fermi能级电离到真空静止能级所需要的能量（与非周期边界下的Koopmans定理是对应的），但是周期边界下没有严格意义上的真空静止能级，我们一般取远离表面且电势平稳的区域（如你第二张图大致横坐标25～35的区域）作为真空静止能级。严格来说，非周期边界条件下，真空静止能级按照电动力学的定义取两个带电粒子无穷远的势能为零点，周期性边界条件下，因为静电作用需要特殊处理导致能级被整体平移，所以只能通过这种方式来进行定义。

qwe1832518773

卡开发发 发表于 2021-1-16 21:17
1、块体模型没有真空，所以没有功函数（模型意义上的），所以graphite那个不要去折腾了。
2、理解基本正 ...

谢谢卡老师回复，我还有一个疑问，那我得到的第二个图里有用的信息只有那个25-30处真空能级的能量吗？前面哪些石墨烯层的能量没用吗

卡开发发

qwe1832518773 发表于 2021-1-16 21:31
谢谢卡老师回复，我还有一个疑问，那我得到的第二个图里有用的信息只有那个25-30处真空能级的能量吗？前 ...

对于功函数计算无用。

qwe1832518773

卡开发发 发表于 2021-1-16 21:40
对于功函数计算无用。

好的感谢卡老师！！

qwe1832518773

卡开发发 发表于 2021-1-16 21:40
对于功函数计算无用。

卡老师，能再麻烦您帮我看一下我另一个帖子吗？一直没得到回复
http://bbs.keinsci.com/thread-21341-1-1.html