|
|
[交流]
【求助】金红石结构TiO2(110)表面体系能带结构计算得到的结果只有一半
计算金红石结构的TiO2(110)面的能带结构时有几个问题想交流一下:
表面第一布里渊区的选取(我这个是按定义推导的,还有文献说表面的布里渊区是体相布里渊区在该方向的投影,这个没弄明白,呵呵。到底怎么选,希望能有人指点一下):
表面可以看作是二维周期结构(A,B方向有周期性),C方向为真空层。那么表面的基矢为a1和a2(A,B方向),可以求得它们的倒格矢b1,b2从而得到第一布里渊区。例如金红石结构的TiO2(110)面的基矢是两个垂直的向量(可由MS导出),它们的倒格矢构成一个长方形网格,第一布里渊区也是长方形,高对称点选为1/4的长方形顶点(与参考文献的一样):例如我选取的:
KPOINTS
k-points along high symmetry lines
10
Line-mode
rec
0 0 !G
0 0.5 !X
0 0.5 ! X
0.5 0.5 ! M
0.5 0.5 ! M
0.5 0 !L
0.5 0 ! L
0 0 !G
正常情况下会算40个K点,但是我的结果只得到了G到X、X到M的20个点,后面的没算(我把G到X、X到M的能带结构画出来发现和文献中能带结构类似)。之后我把后面的 M到L、L到G的20个点又单独算了一下,只得到M到L的10个点的结果,很郁闷。
我把INCAR和EIGENVAL 粘过来,请大家看看哪出错了,请指教。
INCAR和EIGENVAL
SYSTEM = TiO2(110)
ISTART=0
ICHARG=11
ENCUT=400
EDIFF=0.1E-06
EDIFFG=0.1E-03
NSW=0
IBRION=-1
PREC=Accurate
EIGENVAL
TiO2(110)
256 20 156 (20个k点,156条能带) |
» 收录本帖的淘帖专辑推荐
» 猜你喜欢
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有3人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有4人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
» 本主题相关价值贴推荐,对您同样有帮助:
» 抢金币啦!回帖就可以得到:
查看全部散金贴
|