计算化学公社

标题: meta-GGA算DOS图的问题 [打印本页]

作者
Author: 起就起 时间: 2016-7-8 18:23
标题: meta-GGA算DOS图的问题
各位大神好：
   我用meta-GGA算CH2CH2分子的DOS图，跟用别的PBE完全不一样，而且感觉跟CH2CH2理论上的轨道根本对不上，觉得是错误的。有木有人用meta-GGA算过DOS？？或者请指教哪里可能出问题了？
   meta-GGA是RTPSS。
   谢谢。

作者
Author: sobereva 时间: 2016-7-8 19:32
要把情况交代清楚，用什么程序、怎么做的乙烯的DOS图。附件里是MgO的，和乙烯有什么关系？和理论上的对不上，理论的具体指什么？
meta-GGA算的轨道能，以及DOS，并不会有什么问题。

这是g09在TPSS/6-31G*下产生fch，然后由multiwfn绘制的DOS图
(, 下载次数 Times of downloads: 41)

这是PBE下的
(, 下载次数 Times of downloads: 65)

根本没什么区别

作者
Author: 起就起 时间: 2016-7-8 20:03

sobereva 发表于 2016-7-8 19:32
要把情况交代清楚，用什么程序、怎么做的乙烯的DOS图。附件里是MgO的，和乙烯有什么关系？和理论上的对不上 ...

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

作者
Author: 起就起 时间: 2016-7-8 20:07

起就起发表于 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..

作者
Author: 起就起 时间: 2016-7-8 21:09
我用meta-GGA以及不同的ISMEAR和SIGMA算出来的DOS图不太一样，但是DOS图的X轴的值都是-7到2 eV, 但PBE算出来的X轴坐标范围都是-17到10 eV。这是怎么回事呢？

作者
Author: 起就起 时间: 2016-7-8 22:30
本帖最后由起就起于 2016-7-8 22:35 编辑

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

作者
Author: 卡开发发 时间: 2016-7-9 15:19

起就起发表于 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的区间如何？

作者
Author: 起就起 时间: 2016-7-9 16:00

卡开发发发表于 2016-7-9 15:19
有没有试过有动能密度定义的赝势用PBE去算算看？DOS的区间如何？

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

作者
Author: 卡开发发 时间: 2016-7-9 16:03
本帖最后由卡开发发于 2016-7-9 16:06 编辑

起就起发表于 2016-7-9 16:00
算过了，DOS正常。而且，有个问题，用meta-GGA算CH2CH2气相分子能量是负值，但是算金属以及吸附，绝对值 ...

PBE能否正常计算金属以及吸附？你试试其他Meta-GGA，如TPSS，我不排除是泛函（没准程序写的有问题）的问题。

作者
Author: 起就起 时间: 2016-7-18 15:10

卡开发发发表于 2016-7-9 16:03
PBE能否正常计算金属以及吸附？你试试其他Meta-GGA，如TPSS，我不排除是泛函（没准程序写的有问题）的问 ...

好的，我去试试，谢谢。

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