计算化学公社

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

[VASP] meta-GGA算DOS图的问题

[复制链接 Copy URL]

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位大神好:
       我用meta-GGA算CH2CH2分子的DOS图,跟用别的PBE完全不一样,而且感觉跟CH2CH2理论上的轨道根本对不上,觉得是错误的。有木有人用meta-GGA算过DOS??或者请指教哪里可能出问题了?
       meta-GGA是RTPSS。
       谢谢。

dos.pdf

60.86 KB, 下载次数 Times of downloads: 16

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2016-7-18 15:10:19 | 只看该作者 Only view this author
卡开发发 发表于 2016-7-9 16:03
PBE能否正常计算金属以及吸附?你试试其他Meta-GGA,如TPSS,我不排除是泛函(没准程序写的有问题)的问 ...

好的,我去试试,谢谢。

3621

帖子

3

威望

1万

eV
积分
18426

Level 6 (一方通行)

第一原理惨品小作坊

9#
发表于 Post on 2016-7-9 16:03:33 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2016-7-9 16:06 编辑
起就起 发表于 2016-7-9 16:00
算过了,DOS正常。而且,有个问题,用meta-GGA算CH2CH2气相分子能量是负值,但是算金属以及吸附,绝对值 ...

PBE能否正常计算金属以及吸附?你试试其他Meta-GGA,如TPSS,我不排除是泛函(没准程序写的有问题)的问题。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2016-7-9 16:00:15 | 只看该作者 Only view this author
卡开发发 发表于 2016-7-9 15:19
有没有试过有动能密度定义的赝势用PBE去算算看?DOS的区间如何?

算过了,DOS正常。而且,有个问题,用meta-GGA算CH2CH2气相分子能量是负值,但是算金属以及吸附,绝对值都是正值。

3621

帖子

3

威望

1万

eV
积分
18426

Level 6 (一方通行)

第一原理惨品小作坊

7#
发表于 Post on 2016-7-9 15:19:05 | 只看该作者 Only view this author
起就起 发表于 2016-7-8 21:09
我用meta-GGA以及不同的ISMEAR和SIGMA算出来的DOS图不太一样,但是DOS图的X轴的值都是-7到2 eV, 但PBE算出 ...
只是POTCAR有原来的不含kinetic energy density改成了含有kinetic energy density的POTCAR

有没有试过有动能密度定义的赝势用PBE去算算看?DOS的区间如何?
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2016-7-8 22:30:29 | 只看该作者 Only view this author
本帖最后由 起就起 于 2016-7-8 22:35 编辑

我试了以下参数,所有的所有都不变,只在INCAR中设置或者去掉METAGGA=RTPSS这个参数,得到的DOS的区间完全不一样,去掉METAGGA=RTPSS之后,就正常了。这是怎么回事?

Self-consistent meta-generalized gradient approximation within the projector-aug.pdf

265.36 KB, 下载次数 Times of downloads: 25

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2016-7-8 21:09:51 | 只看该作者 Only view this author
我用meta-GGA以及不同的ISMEAR和SIGMA算出来的DOS图不太一样,但是DOS图的X轴的值都是-7到2 eV, 但PBE算出来的X轴坐标范围都是-17到10 eV。这是怎么回事呢?

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2016-7-8 20:07:48 | 只看该作者 Only view this author
起就起 发表于 2016-7-8 20:03
哇塞,果然是大神啊。谢谢你的回答,我来解释一下。
第一,我用的VASP程序,算出vasprun.xml,然后用p4v ...

我觉得PBE算的是对的,可meta-GGA不知道怎么算成这样。我算meta-GGA的DOS跟PBE是一样的,只是POTCAR有原来的不含kinetic energy density改成了含有kinetic energy density的POTCAR,但是H的含有kinetic energy density的POTCAR没有找到,所以H的POTCAR依然用的不含有kinetic energy density的POTCAR。INCAR中的参数只是多了METAGGA=RTPSS,LASPH=.TRUE..

26

帖子

0

威望

55

eV
积分
81

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2016-7-8 20:03:52 | 只看该作者 Only view this author
sobereva 发表于 2016-7-8 19:32
要把情况交代清楚,用什么程序、怎么做的乙烯的DOS图。附件里是MgO的,和乙烯有什么关系?和理论上的对不上 ...

哇塞,果然是大神啊。谢谢你的回答,我来解释一下。
第一,我用的VASP程序,算出vasprun.xml,然后用p4v打开,show DOS, 然后导出DOS图的数据,之后用origin导入数据,做的图。

dos-ch2ch2-totaldos.jpg (83.84 KB, 下载次数 Times of downloads: 57)

PBE的

PBE的

ch2ch2-metaGGA.jpg (139.49 KB, 下载次数 Times of downloads: 52)

这个是meta-GGA算的

这个是meta-GGA算的

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

2#
发表于 Post on 2016-7-8 19:32:17 | 只看该作者 Only view this author
要把情况交代清楚,用什么程序、怎么做的乙烯的DOS图。附件里是MgO的,和乙烯有什么关系?和理论上的对不上,理论的具体指什么?
meta-GGA算的轨道能,以及DOS,并不会有什么问题。

这是g09在TPSS/6-31G*下产生fch,然后由multiwfn绘制的DOS图


这是PBE下的


根本没什么区别
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

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

GMT+8, 2024-11-23 17:12 , Processed in 0.195700 second(s), 25 queries , Gzip On.

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