石墨烯计算问题

镇魂歌

计算石墨烯原胞的时候是零带隙，有狄拉克点，而计算石墨烯超胞的时候会出现带隙，没有在0上出现狄拉克点

卡开发发

以前遇到过，这个问题前面仔细研究了一下，结果还是让人瞠目结舌的：
1、单胞计算：对网格进行平移，使之过（-0.333 0.667 0），结果发觉DMol3的kponits offset只能[-0.5,0.5]，失败；
2、单胞计算：尝试将（-0.333 0.667 0）包含到SCF过程的KPOINTS中，手工修改了kpoints文件，算出来有带隙，失败；
3、单胞计算：修改bandstr_kpoints文件，将（-0.333 0.667 0 间隔）修改成（-0.333333333 0.666666667 0 间隔），带隙0，成功；
4、4x4超胞计算：过程同3，带隙0，成功。

镇魂歌

好的，我先试下，谢谢

镇魂歌

卡开发发 发表于 2016-5-26 21:48
以前遇到过，这个问题前面仔细研究了一下，结果还是让人瞠目结舌的：
1、单胞计算：对网格进行平移，使之 ...

你好，这个是我把石墨烯改成四方结构，那计算这个模型的能带出现了带隙，是取点的原因么

卡开发发

镇魂歌 发表于 2016-5-27 10:10
你好，这个是我把石墨烯改成四方结构，那计算这个模型的能带出现了带隙，是取点的原因么

如果带隙是0.0*这样的级别的，估计就是我说的那种原因，你可以试试。

镇魂歌

卡开发发 发表于 2016-5-27 11:57
如果带隙是0.0*这样的级别的，估计就是我说的那种原因，你可以试试。

六方的已经试出来了，有效

镇魂歌

卡开发发 发表于 2016-5-26 21:48
以前遇到过，这个问题前面仔细研究了一下，结果还是让人瞠目结舌的：
1、单胞计算：对网格进行平移，使之 ...

老师，如果把六方的石墨烯改为四方的，计算的带隙不为零，那band structure.kpoint 应该改哪一个点

卡开发发

镇魂歌 发表于 2016-6-6 14:50
老师，如果把六方的石墨烯改为四方的，计算的带隙不为零，那band structure.kpoint 应该改哪一个点

四方的高对称点应当不同于六方的，你看下G[0,0,0],X[0,0.5,0],M[0.5,0.5,0]。

镇魂歌

卡开发发 发表于 2016-6-6 15:05
四方的高对称点应当不同于六方的，你看下G[0,0,0],X[0,0.5,0],M[0.5,0.5,0]。

还是有带隙

卡开发发

镇魂歌 发表于 2016-6-6 15:12
还是有带隙

我算了个小点的体系没带隙。

你把结构（*.cif）和参数(*.input)发来看看。

镇魂歌

卡开发发 发表于 2016-6-6 15:23
我算了个小点的体系没带隙。

你把结构（*.cif）和参数(*.input)发来看看。

好的

卡开发发

镇魂歌 发表于 2016-6-6 15:37
好的

这个又得说说解决问题的思路，以之前的经验石墨烯做优化结构几乎是不变的，将你的结构和建立的结构直接计算单点能，并且考虑以下问题：
1、对你的结构进行单点能计算：
（1）参考你的input，计算单点能，有带隙，失败；
（2）分别调整多极展开，截断半径，k点，有带隙，失败；
2、按照你的尺寸直接建立的结构进行单点能计算：
（1）参考你的input，计算单点能，无带隙；
3、一样的参数算相同的体系的单点能，为什么结果不同？那可能就是结构不同了，虽然偏差可能不是很大；
（1）参考你的input优化直接建立的结构，有带隙，显然就是优化结构导致的；
（2）分析你的设置中不可靠的因素，其中我认为截断半径出现问题的可能最大，一般建议截断半径还是>4A吧
（3）以你的结构作为初始结构，对你的input仅修改cutoff进行计算，无带隙，成功。
另外还有一点，前面有个地方我也弄错了，这个结构不是正方的，所以高对称点应该要G-X-M-Y[0,0.5,0]。

学会解决问题的方法很重要。

镇魂歌

卡开发发 发表于 2016-6-6 17:45
这个又得说说解决问题的思路，以之前的经验石墨烯做优化结构几乎是不变的，将你的结构和建立的结构直接计 ...

麻烦您了，谢谢老师。老师，您说的根据我的结构尺寸直接建立的结构是根据晶胞的大小来构建的么

卡开发发

镇魂歌 发表于 2016-6-7 10:00
麻烦您了，谢谢老师。老师，您说的根据我的结构尺寸直接建立的结构是根据晶胞的大小来构建的么

恩，一模一样的超胞大小。

镇魂歌

卡开发发 发表于 2016-6-7 10:27
恩，一模一样的超胞大小。

好的，再次谢谢老师