meta-GGA算DOS图的问题

起就起

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

sobereva

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

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


这是PBE下的


根本没什么区别

起就起

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

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

起就起

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

起就起

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

起就起

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

卡开发发

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

起就起

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

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

卡开发发

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

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

起就起

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

好的，我去试试，谢谢。