24小时热门版块排行榜    

查看: 2272  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xijinyang

金虫 (小有名气)


[交流] 【求助】怎么确定abinit计算能带结构的费米能位置

我想计算graphene的能带结构,然后就按教程的做了下,具体的输入文件是:

ndtset 2

#Dataset 1 : usual self-consistent calculation
kptopt1 1          # Option for the automatic generation of k points,
                   # taking into account the symmetry
nshiftk1 1
shiftk1  0.0 0.0 0.5  # These shifts will be the same for all grids
ngkpt1  36 36 1  
prtden1  1         # Print the density, for use by dataset 2
tolvrs1  1.0d-20   # This value is way too large for most realistic studies of materials
nband1  16
occopt1 3
iscf1 7
#Dataset 2 : the band structure
iscf2    -2
getden2  -1
kptopt2  -3
nband2   16
ndivk2   15 15 15      # 10, 12 and 17 divisions of the 3 segments, delimited
                       # by 4 points.
kptbounds2  0.0  0.0  0.0 # gamma point
            0.5  0.0  0.0 # M point
            1/3  1/3  0.0 # K point
            0.0  0.0  0.0 # Gamma point
tolwfr2  1.0d-12
enunit2  1             # Will output the eigenenergies in eV


#Definition of the unit cell
acell  4.6209 4.6209 9.44863                  
rprim  1.0  0.0       0.0        # The hexagonal axes
      -0.5  0.8660    0.0
       0.0  0.0       1.0

#Definition of the atom types
ntypat 1          # There is only one type of atom
znucl 6           # The keyword "znucl" refers to the atomic number of the
                  # possible type(s) of atom. The pseudopotential(s)
                  # mentioned in the "files" file must correspond
                 
                        

#Definition of the atoms
natom 2           # There are two atoms
typat 1 1         
xred   2/3   1/3  0.0  # Triplet giving the REDUCED coordinate of atom 1.
       1/3   2/3  0.0  # Triplet giving the REDUCED coordinate of atom 2.
#Definition of the planewave basis set
ecut 30.0         # Maximal kinetic energy cut-off, in Hartree

#Definition of the SCF procedure
nstep 100          # Maximal number of SCF cycles
diemac 12.0      

计算完,我画的能带图如下:

可是我不知道fermi能在哪个位置?是在0点么?如果是0点的话,怎么我的dirac锥点位置比0点要高啊?是不是我参数设置有问题?
谢谢大家!非常感谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)



gzqdyouxia(金币+1):鼓励交流 2011-01-07 13:42:53
引用回帖:
Originally posted by mazuju028 at 2010-12-30 10:50:34:
计算结果中有Feimi energy啊

请教怎么确定费米能级和导带与价带的相对位置?
4楼2011-01-06 15:26:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

frank_zhan

金虫 (正式写手)


★ ★
youzhizhe(金币+1):鼓励交流。 2010-12-30 10:39:14
zhang668(金币+1):多谢交流 2010-12-30 19:53:59
我觉得会不会是因为在定义unit cell时acell  4.6209 4.6209 9.44863的问题。因为 c方向上9.44863bhor相当于0.5 nm,是不是不足够消除层和层之间的 相互作用?我觉得这个方向的长度至少要超过2nm才能说是单层的graphene吧。
2楼2010-12-30 10:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
zzy870720z(金币+1):谢谢提示 2010-12-30 17:37:00
zhang668(金币+1):多谢交流 2010-12-30 19:54:12
计算结果中有Feimi energy啊
3楼2010-12-30 10:50:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)


同问,我也想知道!
5楼2011-01-06 19:47:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见