| 查看: 3273 | 回复: 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 |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有48人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有29人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
12楼2014-07-05 20:44:43
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 (幼儿园)
- 金币: 4546.6
- 散金: 723
- 红花: 5
- 帖子: 720
- 在线: 144.7小时
- 虫号: 2418920
- 注册: 2013-04-15
- 专业: 宇宙学
3楼2014-06-28 10:44:50
紫海心蝶
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4546.6
- 散金: 723
- 红花: 5
- 帖子: 720
- 在线: 144.7小时
- 虫号: 2418920
- 注册: 2013-04-15
- 专业: 宇宙学
5楼2014-06-30 21:05:17
