| 查看: 2854 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
紫海心蝶木虫 (正式写手)
|
[求助]
计算石墨烯超胞能带居然出现很大的带隙！求指教 已有1人参与
|
||
|
计算石墨烯的能带时，用超胞和原胞分别进行计算，INCAR参数相同，高对称点路径一致，为什么用原胞计算的能带没有带隙，用超胞计算的能带却出现很大的带隙呢？求指教！！！ 这是INVAR: SYSTEM=GRAPHENE PREC= Accurate ENCUT= 500 ISTART= 1 ICHARG= 11 ISMEAR= 0 SIGMA= 0.2 EDIFF= 1e-5 LORBIT= 10 这是KPOINTS: k-points 20 Line-mode Rec 0.000000000000 0.000000000000 0.000000000000 !G 0.000000000000 0.500000000000 0.000000000000 !M 0.000000000000 0.500000000000 0.000000000000 !M -0.333333333333 0.666666666667 0.000000000000 !K -0.333333333333 0.666666666667 0.000000000000 !K 0.000000000000 0.000000000000 0.000000000000 !G 原胞POSCAR: New structure 1.00000000000000 2.4599988461000000 0.0000000000000000 0.0000000000000000 -1.2299994230000000 2.1304214940000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 20.0000000000000000 C 2 Direct 0.3333333429999996 0.6666666269999979 0.0000000000000000 0.6666666269999979 0.3333333429999996 0.0000000000000000 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 超胞POSCAR: raphene 1.00000000000000 0.0000000000000000 4.9199978407307583 0.0000000000000000 4.2608431166374219 2.4599989203653791 0.0000000000000000 0.0000000000000000 0.0000000000000000 -20.0000000000000000 C 8 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.1666666666666643 0.1666666666666643 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.6666666666666643 0.6666666666666643 0.0000000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.6666666666666643 0.1666666666666643 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.0000000000000000 0.1666666666666643 0.6666666666666643 0.0000000000000000 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 请大家帮我分析下是什么原因  超胞.jpg  原胞.jpg |
回复此楼» 猜你喜欢
307求调剂
已经有13人回复
-
071010 323 分求调剂
已经有3人回复
-
材料调剂
已经有9人回复
-
材料专硕调剂
已经有10人回复
-
一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历
已经有8人回复
-
求调剂
已经有8人回复
-
284求调剂
已经有15人回复
-
085602化工求调剂（331分）
已经有3人回复
-
一志愿南昌大学324求调剂
已经有8人回复
-
材料科学与工程调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
|
可能是k-mesh网格不足，见下面的这篇文章也犯了同样的问题 http://journals.aps.org/prb/abstract/10.1103/PhysRevB.78.075435 被comment了之后，作者的reply提到了"k-mesh网格不足，导致他们发现了应变能打开石墨烯的带隙这样错误的结论” http://journals.aps.org/prb/abstract/10.1103/PhysRevB.80.167401 http://journals.aps.org/prb/abstract/10.1103/PhysRevB.80.167402 |
14楼2014-10-13 09:54:27
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
紫海心蝶: 金币+5, ★★★★★最佳答案 2014-06-28 10:35:21
感谢参与，应助指数 +1
紫海心蝶: 金币+5, ★★★★★最佳答案 2014-06-28 10:35:21
|
你采用的是2x2的超原胞吧. 超原胞的基矢写法优点奇怪, 特别是基矢a3, 用的是负值,另外 你的超原胞的另外两个基矢a1取在y轴, a2基矢取在第一卦相空间(x>0, y>0),这样这两个基矢的夹角是60度. 而你原胞的基矢a1取在x轴上,a2基矢取在第二卦相空间(x<0, y>0).这两个基矢的夹角是120度. 所以它们不一致. 另外超原胞的布里渊区会折叠, 对石墨烯原胞的K和K'(Dirac点所在的k点)并不会折叠到 2x2的超原胞的布里渊区中心.对2x2超原胞的石墨烯来计算能带时,高对称k点之间的间隔可能要取密些. |
2楼2014-06-28 09:50:41
紫海心蝶
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4540.8
- 散金: 723
- 红花: 5
- 帖子: 718
- 在线: 144.7小时
- 虫号: 2418920
- 注册: 2013-04-15
- 专业: 宇宙学
3楼2014-06-28 10:44:50
紫海心蝶
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4540.8
- 散金: 723
- 红花: 5
- 帖子: 718
- 在线: 144.7小时
- 虫号: 2418920
- 注册: 2013-04-15
- 专业: 宇宙学
5楼2014-06-30 21:05:17
