24小时热门版块排行榜    

查看: 1404  |  回复: 6
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 xijinyang 的 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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的回帖
相关版块跳转 我要订阅楼主 xijinyang 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知 +7 中南大学-金属材 2024-11-14 18/900 2024-11-15 20:22 by 中南大学-金属材
[论文投稿] angew三个小修,返修过去已经20多天了没动静,有朋友遇到这种情况吗? +5 迟迟未到场 2024-11-14 5/250 2024-11-15 20:16 by 秦时明月zy
[考博] 随缘读博 一篇中科院一区Top 一篇中科院二区Top 两篇专利 +6 Ricoch4t 2024-11-13 15/750 2024-11-15 19:24 by 青青之洲
[论文投稿] 找不到审稿人 +6 吉祥如意141913 2024-11-13 6/300 2024-11-15 16:48 by LEyblyt
[论文投稿] 职称论文 +10 xiaoqi_66 2024-11-09 13/650 2024-11-15 16:45 by sugarmei
[论文投稿] 核心期刊那个阶段最难 +4 公共茅斯 2024-11-11 4/200 2024-11-15 16:44 by sugarmei
[论文投稿] 小论文影响大论文查重吗? +3 阿道夫深度 2024-11-14 4/200 2024-11-15 16:19 by sugarmei
[基金申请] 博后面上太难了 +14 不安定因素 2024-11-14 17/850 2024-11-15 16:05 by lupbao
[教师之家] 大学老师 +11 考研一路顺风 2024-11-13 14/700 2024-11-15 15:24 by LNP@mRNA
[基金申请] 博后资助名单出来了 +12 Shxjjxjkx 2024-11-14 14/700 2024-11-15 14:40 by fanfenggui
[基金申请] 去年七月底入站的还能申请下一批吗? +4 brightwo 2024-11-14 4/200 2024-11-15 14:18 by gazi1111
[硕博家园] 实验进度 +9 麻辣烫。 2024-11-11 14/700 2024-11-15 13:14 by hanweiguo
[基金申请] 广东省自然科学基金-面上项目消息有了么? +4 pachang 2024-11-14 6/300 2024-11-15 09:39 by zhoulei7170
[硕博家园] 课题和专业一点都不相关怎么办 +7 imt4n9 2024-11-11 10/500 2024-11-15 01:55 by 成长虫
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[基金申请] 浙江省青基怎么看A还是B? +8 ToallGoof 2024-11-11 14/700 2024-11-13 09:11 by Praybest
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[论文投稿] 论文接受后,编辑部发邮件提醒缺少附件 15+4 上善若水明泪 2024-11-10 4/200 2024-11-12 14:24 by 北京莱茵润色
信息提示
请填处理意见