功函数计算的LOCPOT怎么分析呢

qwe1832518773 · 发表于 Post on 2021-1-11 13:42:14

设置LVTOT= .TRUE. 得到了LOCPOT，打开一看里面有五列数据，分别是什么意思呢？

小*虫上有帖子说是用脚本vtotav.f处理，“用ifort -o vtotav.x vtotav.f，然后用./vtov.x命令，选取你要沿的某一个方向”（http://mu****g.com/html/201403/7141305.html）
但是我没找到这个脚本，并且看这数据格式挺简单的，是否能自己处理呢？

还有一个关于Fermi能级计算的问题，我在功函数这个文件夹里直接grep Fermi OUTCAR，出现

这里的哪一个数据是费米能级呢？是1.989042 eV嘛？
小弟VASP菜鸟，希望大佬们可以指导我一下，不胜感激

qwe1832518773 · 发表于 Post on 2021-1-15 14:43:41

hakuna 发表于 2021-1-12 23:39
grep E-fermi
1.989042应该是Thomas-Fermi screening length，不是Fermi能级，单位也不对

感谢回复！

hakuna · 发表于 Post on 2021-1-12 23:39:53

grep E-fermi
1.989042应该是Thomas-Fermi screening length，不是Fermi能级，单位也不对

qwe1832518773 · 发表于 Post on 2021-1-12 11:23:36

hakuna 发表于 2021-1-12 00:25
LOCPOT是个静电势文件，你截的数据是slab里各点的静电势，是一个108x128x300的三维矩阵，也可以看成是a、b ...

好的谢谢大佬讲解。想再追问一下，我看网上搜费米能级找的是E_Fermi 这个词，但我用grep Fermi OUTCAR只能找到我贴出来的那两行，那我的费米能级是1.989042eV吗

qwe1832518773 · 发表于 Post on 2021-1-12 11:20:12

ene 发表于 2021-1-11 15:26
这种东西上谷狗搜一下就能找到

还是当了回伸手党吗

有点不好意思，非常感谢您的回复

hakuna · 发表于 Post on 2021-1-12 00:25:25

LOCPOT是个静电势文件，你截的数据是slab里各点的静电势，是一个108x128x300的三维矩阵，也可以看成是a、b平面上108x128两维矩阵，这样的二维矩阵沿c方向共300个，对每层a、b上的数据求平均，共300个数据，然后对c做图，真空区平坦处即真空电势，与费米能级的差值即功函数

ene · 发表于 Post on 2021-1-11 15:26:40

这种东西上谷狗搜一下就能找到

vtotav.f (3.96 KB, 下载次数 Times of downloads: 428)

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

[VASP] 功函数计算的LOCPOT怎么分析呢

评分 Rate

评分 Rate