计算化学公社

 找回密码 Forget password
 注册 Register
Views: 4023|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

91

帖子

0

威望

6067

eV
积分
6158

Level 6 (一方通行)

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

文献的结果:

同时附上文献、LOCPOT: Grand-Canonical Model of Electrochemical Double Layers from a Hybrid Density–Po.pdf (972.04 KB, 下载次数 Times of downloads: 19) CHGCAR.rar (7.16 MB, 下载次数 Times of downloads: 15) LOCPOT.rar (6.17 MB, 下载次数 Times of downloads: 6)

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

2#
发表于 Post on 2022-5-25 21:54:04 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“Multiwfn电子密度和电势分布”改了,以后务必注意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

91

帖子

0

威望

6067

eV
积分
6158

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2022-5-25 22:26:49 | 只看该作者 Only view this author
sobereva 发表于 2022-5-25 21:54
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此 ...

谢谢Sob老师,下次注意

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

4#
发表于 Post on 2022-5-26 04:30:35 | 只看该作者 Only view this author
我不知道你怎么操作的
用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)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

91

帖子

0

威望

6067

eV
积分
6158

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2022-5-26 09:50:49 | 只看该作者 Only view this author
本帖最后由 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.



6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

6#
发表于 Post on 2022-5-27 04:00:57 | 只看该作者 Only view this author
Jack 发表于 2022-5-26 09:50
谢谢Sob老师,是我操作问题,Multiwfn可以很好的处理这个。想请问下,关于这里的电势问题。当电势单位为 ...

也许文献里做了平移,以令远离界面的曲线水平处的数值为0。我没细看文章,只是个猜测
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

91

帖子

0

威望

6067

eV
积分
6158

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2022-5-27 08:21:07 | 只看该作者 Only view this author
sobereva 发表于 2022-5-27 04:00
也许文献里做了平移,以令远离界面的曲线水平处的数值为0。我没细看文章,只是个猜测

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

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

8#
发表于 Post on 2022-5-27 08:42:02 | 只看该作者 Only view this author
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.了,当前单位直接就和文献里的图对应。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

91

帖子

0

威望

6067

eV
积分
6158

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2022-5-27 10:10:09 | 只看该作者 Only view this author
sobereva 发表于 2022-5-27 08:42
https://www.vasp.at/wiki/index.php/LOCPOT
LOCPOT是用eV作为单位记录的,和一般格点文件习俗不同
但M ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 06:43 , Processed in 0.193277 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list