计算化学公社

标题: 求助VASP计算Ag(111)后用Multiwfn处理电子密度和外势信息和文献差别较大 [打印本页]

作者
Author: Jack 时间: 2022-5-25 21:32
标题: 求助VASP计算Ag(111)后用Multiwfn处理电子密度和外势信息和文献差别较大
大家好，我用VASP计算了Ag(111)表面的电子密度电势分布，采用Multiwfn分别处理包含电子密度和外势信息的LOCPOT。但是与文献差别较大，还请大家指导！
我计算的结果：
(, 下载次数 Times of downloads: 35)
文献的结果：
(, 下载次数 Times of downloads: 21)
同时附上文献、LOCPOT： (, 下载次数 Times of downloads: 19) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 6)

作者
Author: sobereva 时间: 2022-5-25 21:54
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“Multiwfn电子密度和电势分布”改了，以后务必注意

作者
Author: Jack 时间: 2022-5-25 22:26

sobereva 发表于 2022-5-25 21:54
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

谢谢Sob老师，下次注意

作者
Author: sobereva 时间: 2022-5-26 04:30
我不知道你怎么操作的
用Multiwfn官网上最新版本，载入CHGCAR，依次输入
y
13
18
Z
a
3
得到下面的平面平均电子密度曲线
(, 下载次数 Times of downloads: 23)

可见峰位置大约0.183 a.u.，具体来说是Bohr^(-3)。转换成Angstrom^(-3)就是0.183*0.5292^(-3)=1.23，可见和文献里的DFT结果基本一致（文献大概1.3）

作者
Author: Jack 时间: 2022-5-26 09:50
本帖最后由 Jack 于 2022-5-26 10:12 编辑

sobereva 发表于 2022-5-26 04:30
我不知道你怎么操作的
用Multiwfn官网上最新版本，载入CHGCAR，依次输入
y

谢谢Sob老师，是我操作问题，Multiwfn可以很好的处理这个。想请问下，关于这里的电势问题。当电势单位为原子单位a.u.时与V之间的换算关系是1a.u. = 27.2114 V么，是我忽略了哪一点或者换算错了导致数据与文献不对应。另外我用Multiwfn算出有负值，且随着真空层厚度增加而变小。真空层厚度为20埃的时候是-5 a.u.，厚度为40埃时大概为-3 a.u.
(, 下载次数 Times of downloads: 20)

作者
Author: sobereva 时间: 2022-5-27 04:00

Jack 发表于 2022-5-26 09:50
谢谢Sob老师，是我操作问题，Multiwfn可以很好的处理这个。想请问下，关于这里的电势问题。当电势单位为 ...

也许文献里做了平移，以令远离界面的曲线水平处的数值为0。我没细看文章，只是个猜测

作者
Author: Jack 时间: 2022-5-27 08:21

sobereva 发表于 2022-5-27 04:00
也许文献里做了平移，以令远离界面的曲线水平处的数值为0。我没细看文章，只是个猜测

谢谢Sob老师，我也是这样猜测的，而且也是合理的，真空层厚度无穷大时应该是0。但是这个原子单位a.u.换算过来的话，不知道怎么处理能与文献数值吻合。

作者
Author: sobereva 时间: 2022-5-27 08:42

Jack 发表于 2022-5-27 08:21
谢谢Sob老师，我也是这样猜测的，而且也是合理的，真空层厚度无穷大时应该是0。但是这个原子单位a.u.换算 ...

https://www.vasp.at/wiki/index.php/LOCPOT
LOCPOT是用eV作为单位记录的，和一般格点文件习俗不同
但Multiwfn假定格点数据是a.u.作为单位，所以你就别管Multiwfn此时显示的a.u.了，当前单位直接就和文献里的图对应。

作者
Author: Jack 时间: 2022-5-27 10:10

sobereva 发表于 2022-5-27 08:42
https://www.vasp.at/wiki/index.php/LOCPOT
LOCPOT是用eV作为单位记录的，和一般格点文件习俗不同
但M ...

确实，我处理后发现二者数值差不多。谢谢Sob老师耐心解答！

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