计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 求助VASP计算Ag(111)后用Multiwfn处理电子密度和外势信 ...

求助VASP计算Ag(111)后用Multiwfn处理电子密度和外势信息和文献差别较大

查看数: 4026 | 评论数: 8 | 收藏 Add to favorites 4

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Jack

发布时间: 2022-5-25 21:32

正文摘要:

大家好，我用VASP计算了Ag(111)表面的电子密度电势分布，采用Multiwfn分别处理包含电子密度和外势信息的LOCPOT。但是与文献差别较大，还请大家指导！我计算的结果：

回复 Reply

Jack 发表于 Post on 2022-5-27 10:10:09

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

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

sobereva 发表于 Post on 2022-5-27 08:42:02

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.了，当前单位直接就和文献里的图对应。

Jack 发表于 Post on 2022-5-27 08:21:07

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

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

sobereva 发表于 Post on 2022-5-27 04:00:57

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

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

sobereva 发表于 Post on 2022-5-26 04:30:35

我不知道你怎么操作的
用Multiwfn官网上最新版本，载入CHGCAR，依次输入
y
13
18
Z
a
3
得到下面的平面平均电子密度曲线

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

Jack 发表于 Post on 2022-5-25 22:26:49

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

谢谢Sob老师，下次注意

sobereva 发表于 Post on 2022-5-25 21:54:04

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

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