计算化学公社

标题: vaspkit得到的平均静电势文件的费米能级位置？ [打印本页]

作者
Author: pukkie 时间: 2022-1-23 21:18
标题: vaspkit得到的平均静电势文件的费米能级位置？
各位老师好。我在用vasp计算一个金属-半导体接触，得到体系的LOCPOT文件后，使用vaspkit 426对z方向静电势做平均，得到PLANAR_AVERAGE.dat文件。
我的问题是：直接用这个文件画图，体系的费米能级能否认为是在0处？为什么？
我个人理解是不行，因为从该文件的真空层处可以得到真空能级位置，如果0处是费米能级，则功函数值=真空能级能量，这不合理。
但师兄说vaspkit 426已经将费米能级固定为零。请各位老师指点迷津，谢谢。

作者
Author: wsz 时间: 2022-1-25 22:34
你的理解是对的

作者
Author: LI7M 时间: 2023-1-30 16:34
前辈你好！请问要画静电势分布图，是用vaspkit处理静态计算得到的LOCPOT文件呢，还是能带计算的呢

作者
Author: pukkie 时间: 2023-2-9 16:38

LI7M 发表于 2023-1-30 16:34
前辈你好！请问要画静电势分布图，是用vaspkit处理静态计算得到的LOCPOT文件呢，还是能带计算的呢

我用的静态计算的LOCPOT再vaspkit-426取均值

作者
Author: duchenggong 时间: 2023-2-19 19:54
本帖最后由 duchenggong 于 2023-2-19 20:14 编辑

是不是有几个公式可以把相关的几个都计算出来
计算静电势完成后，grep E-Fermi OUTCAR ,可以得到费米能级EF
功函数W=X-EF，X是将静电式文件拖到p4vasp里得到的average
Evbm=-W+0.5Eg
Ecbm=-W-0.5Eg
个人见解，如有不对的地方还请老师们批评指正

作者
Author: asd9998876 时间: 2023-2-20 11:33
VASPKIT426只是格点文件处理，费米能级还是要减去的。不过单个费米能级具体数值没有什么意义，把他取为零，那功函数就是真空能级的能量这个是很合理的

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