计算化学公社

标题: CP2K计算dos的问题求助 [打印本页]

作者
Author: zhou9527 时间: 2021-12-8 10:41
标题: CP2K计算dos的问题求助
各位老师，同学大家好！据我所知，在CP2K中有两种方式可以得到带隙。第一种在MD模拟的时候，可以输出HOMO-LUMO Gap，计算的结果基本和VASP （2.2）是一致的,参数如下：
&DFT
&PRINT
      &MO_CUBES
         NHOMO 1
         NLUMO 1
      &END
   &END
&END
同理，在优化和单点同样也能输出所谓的HOMO-LUMO Gap, 但是如果采用PDOS计算得到相应的Gap,采用:
&DFT
&PRINT
   &PDOS
      NLUMO -1
   &END
&END
&END
通过converlution.py处理之后，得到的带隙（0.5 eV),请问为什么这两个Gap的结果相差这么大呢？
还有一个问题就是，这两个能级是不是一致的呢，类似VASP中的EIGENVL和DOSCAR？
希望，有清楚的老师、同学不吝赐教！十分感谢！

作者
Author: sobereva 时间: 2021-12-8 11:34
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“CP2K dos”改了，以后务必注意，下次将删帖+扣分处理。

作者
Author: zhou9527 时间: 2021-12-8 15:09

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

好的，记下了，麻烦sob老师了

作者
Author: sobereva 时间: 2021-12-9 06:24
如果晶胞足够大，是用gamma点算的，直接用最低空轨道和最高占据轨道能量差当做带隙就完了，不牵扯到计算PDOS当前问题和做MD也没有丝毫直接联系

作者
Author: zhou9527 时间: 2021-12-9 19:34
sob老师，您好！晶胞足够大也是用Gamma点算的,但是我想对MD中的特定一帧做电子结构分析，需要计算相应的DOS，因为用VASP（vasp_gam）还是相对慢一些，但是用CP2K算出来的dos和在MD过程中输出的HOMO-LUMO Gap的数值差距比较大，所以不甚清楚原因在哪

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