CP2K绘制石墨烯能带图如何得到合理的结果

flyingchow · 发表于 Post on 2022-6-11 04:32:13

各位老师好，我这两天学习如何绘制能带，我通过cp2k网站上的用QE绘制能带的教程顺利地完成了石墨烯能带的绘制。
我然后又用DFTB+绘制了石墨烯的能带，和QE的结果是相似的。
图如下：（左图QE绘制，右图DFTB+绘制）

但是当我想用CP2K绘制能带图的时候，得到的结果和上面两个有很大区别。我用的都是相同的k点对称路径，如下图：

我把我的cp2k输入文件附上，请老师们帮我看看有什么问题？谢谢。

graphene.inp (2.94 KB, 下载次数 Times of downloads: 45)

rabbitkiller · 发表于 Post on 2022-6-11 07:12:51

cp2k的图费米能级在哪？scf收敛了么？感觉好奇怪……

gog · 发表于 Post on 2022-6-11 07:39:52

主要是赝势库不同吧。算能带，最准确的，还是要纯粹的平面波基组。castep和qe的赝势也比较全了。

flyingchow · 发表于 Post on 2022-6-11 12:17:19

rabbitkiller 发表于 2022-6-11 07:12
cp2k的图费米能级在哪？scf收敛了么？感觉好奇怪……

收敛了，你可以自己拿我的输入文件算一个看看。

flyingchow · 发表于 Post on 2022-6-11 12:18:34

gog 发表于 2022-6-11 07:39
主要是赝势库不同吧。算能带，最准确的，还是要纯粹的平面波基组。castep和qe的赝势也比较全了。

除了基组和赝势的原因，还有别的原因吗？

gog · 发表于 Post on 2022-6-11 16:17:04

flyingchow 发表于 2022-6-11 12:18
除了基组和赝势的原因，还有别的原因吗？

其他的，我就不知道了，说不上来。功底不够深。

yoghurt117 · 发表于 Post on 2026-1-14 08:48:10

同问，我用网上推荐的DZVP基组优化晶格，结果算出来边长2.62，比2.46大了8%左右，能带图画出来也不太对，和楼主的差不多，不知道怎么回事

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] CP2K绘制石墨烯能带图如何得到合理的结果