计算化学公社

标题: CP2K绘制石墨烯能带图如何得到合理的结果 [打印本页]

作者
Author: flyingchow 时间: 2022-6-11 04:32
标题: CP2K绘制石墨烯能带图如何得到合理的结果
各位老师好，我这两天学习如何绘制能带，我通过cp2k网站上的用QE绘制能带的教程顺利地完成了石墨烯能带的绘制。
我然后又用DFTB+绘制了石墨烯的能带，和QE的结果是相似的。
图如下：（左图QE绘制，右图DFTB+绘制）
(, 下载次数 Times of downloads: 16)
但是当我想用CP2K绘制能带图的时候，得到的结果和上面两个有很大区别。我用的都是相同的k点对称路径，如下图：
(, 下载次数 Times of downloads: 16)
我把我的cp2k输入文件附上，请老师们帮我看看有什么问题？谢谢。
(, 下载次数 Times of downloads: 45)

作者
Author: rabbitkiller 时间: 2022-6-11 07:12
cp2k的图费米能级在哪？scf收敛了么？感觉好奇怪……

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

作者
Author: flyingchow 时间: 2022-6-11 12:17

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

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

作者
Author: flyingchow 时间: 2022-6-11 12:18

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

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

作者
Author: gog 时间: 2022-6-11 16:17

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

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

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

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