关于cp2k计算晶体态密度的一些基本概念和方法请教

WEILANDREW

大家晚上好，我是第一性原理计算新手，想做一些关于态密度的计算，但是对于一些基本概念还不太确定。比如，我用cp2k在gamma点对超胞进行单点计算，输出的gamma点各原子的PDOS和DOS。那么
（1）这时总的DOS是不是严格等于其组成原子各轨道的PDOS加和？比如在整个能级范围内是否 TDOS(Fe2O3) = PDOS(Fe) + PDOS(O)?   是否 Fe原子的分波态密度PDOS(total) = Fe原子的各电子轨道的PDOS(s+p+d+f+...) 求和?
（2）在用cp2k计算态密度时怎样选取泛函呢？用PBE泛函可以吗？
（3）用cp2k输出的DOS有没有后处理的脚本呢？
（4）用cp2k计算态密度，推荐是用扩胞计算还是 用单胞结合K点来计算？两者的结果是一致的吗？
（5）下面的α-Fe2O3是不是没有能隙？为啥在费米能级附近偏上方的DOS不为0呢？这个是materials project的结果，这个结果可靠吗？

非常感谢大家的解答！

sobereva

1 是。是
2 看体系。该用杂化泛函的用杂化泛函
3 按下文产生molden文件，作为Multiwfn输入文件用Multiwfn绘制DOS是最方便灵活的
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651（http://bbs.keinsci.com/thread-31827-1-1.html）
使用Multiwfn绘制态密度(DOS)图考察电子结构
http://sobereva.com/482（http://bbs.keinsci.com/thread-13242-1-1.html）

4 目前版本的CP2K直接产生DOS的功能无法考虑k点。而导出molden文件的功能也不支持k点。也即只能用超胞做gamma点计算

5 你可以Google查查此体系。从当前的图来说是没有能隙