计算化学公社

标题: 功函数计算的LOCPOT怎么分析呢 [打印本页]

作者
Author: qwe1832518773 时间: 2021-1-11 13:42
标题: 功函数计算的LOCPOT怎么分析呢
设置LVTOT= .TRUE. 得到了LOCPOT，打开一看里面有五列数据，分别是什么意思呢？ (, 下载次数 Times of downloads: 48)
小*虫上有帖子说是用脚本vtotav.f处理，“用ifort -o vtotav.x vtotav.f，然后用./vtov.x命令，选取你要沿的某一个方向”（http://mu****g.com/html/201403/7141305.html）
但是我没找到这个脚本，并且看这数据格式挺简单的，是否能自己处理呢？

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

(, 下载次数 Times of downloads: 51)
这里的哪一个数据是费米能级呢？是1.989042 eV嘛？
小弟VASP菜鸟，希望大佬们可以指导我一下，不胜感激

作者
Author: ene 时间: 2021-1-11 15:26
这种东西上谷狗搜一下就能找到
(, 下载次数 Times of downloads: 428)

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

作者
Author: qwe1832518773 时间: 2021-1-12 11:20

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

还是当了回伸手党吗

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

作者
Author: qwe1832518773 时间: 2021-1-12 11:23

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

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

作者
Author: hakuna 时间: 2021-1-12 23:39
grep E-fermi
1.989042应该是Thomas-Fermi screening length，不是Fermi能级，单位也不对

作者
Author: qwe1832518773 时间: 2021-1-15 14:43

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

感谢回复！

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