24小时热门版块排行榜    

查看: 1906  |  回复: 2

ewan

新虫 (正式写手)

[交流] LiF能带结构计算的困惑

大家好!小弟之前没有接触过能带计算,最近使用siesta计算64个原子LiF的能带结构(2*2*2的Supercell),对交换关联泛函是采用LDA,计算沿布里渊区的高对称轴Γ->X->W->Γ计算能带结构,k sampling只有一个Г点,并且MeshCutoff =300.0 Ry。我有点不能理解算出来的结果(如下图),从E(k)似乎有两个带隙(LiF带隙的实验值是14 eV),请问哪个是真正的带隙?最好是指明这两个带隙是如何产生的?(Fermi能级已经移至0 eV处)。如果哪位需要图的源文件的话,可以私聊啊!
        附有能带图两张(第一张是的能量范围是-22~6 eV,第二张是-2~-6 eV)和siesta计算的输入文件。

#######################################################################################
SystemName          bulk.LiF
SystemLabel         bulk_LiF

NumberOfAtoms       8
NumberOfSpecies     2

%block ChemicalSpeciesLabel
1    3    Li
2    9    F
%endblock ChemicalSpeciesLabel

# basis
PAO.BasisSize       DZP
PAO.EnergyShift     200 meV

LatticeConstant    4.04 Ang   # Lattice constant alat
%block LatticeVectors               # Lattice vectors, in units of latticeconstant
2.0      0.0       0.0
0.0      2.0       0.0
0.0      0.0       2.0        
%endblock LatticeVectors


%block SuperCell
2 0 0   0.0
0 2 0   0.0
0 0 2   0.0
%endblock SuperCell


%block kgrid_Monkhorst_Pack
1 0 0   0.0
0 1 0   0.0
0 0 1   0.0
%endblock kgrid_Monkhorst_Pack


BandLinesScale  ReciprocalLatticeVectors
%block BandLines
1     0.0   0.0   0.0     # Gamma-point
20    0.5   0.0   0.5     # X-point
30    0.5   0.25  0.75    # W-point
30    0.0   0.0   0.0     # Gamma-point
%endblock BandLines


%block  ProjectedDensityOfStates
-20 10 0.1  300    eV
%endblock ProjectedDensityOfStates

MeshCutoff          300.0 Ry

# SCF options
MaxSCFIterations      100           # Maximum number of SCF iter
DM.NumberPulay         3            # One Pulay every 3 iter
DM.MixingWeight       0.35         # New DM amount for next SCF cycle
DM.Tolerance          5.d-4         # Tolerance in maximum difference

SolutionMethod        diagon        # OrderN or Diagon
OccupationFunction    FD
ElectronicTemperature  25 meV       # Temp. for Fermi smearing

# exchange-correlation functional
XC.functional LDA
XC.authors CA



# structural infomation
AtomCoorFormatOut Ang

# output information
WriteCoorXmol T
WriteMDXmol T
WriteForces T

WriteKpoints            .true.
WriteEigenvalues        .true.
WriteKbands             .true.
WriteBands              .true.
WriteMullikenPop        1

# options for MD
MD.MaxForceTol          0.01 eV/Ang
MD.TypeOfRun CG
MD.NumCGSteps 400
MD.Broyden.History.Steps 6
MD.Broyden.Initial.Inverse.Jacobian 1.0

# options for continue
# MD.UseSaveCG T                        # neccessary!
# MD.UseSaveXV T                        # neccessary!
# DM.UseSaveDM          true            # to use continuae..ion files


AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
  -5.963747      -2.022693      -2.026392               1
  -6.051639      -5.949852      -2.013945               2
  -2.147364      -2.056847      -2.009531               2
  -2.039835      -5.953723      -2.024337               1
  -6.085698      -1.925561       2.008577               2
  -5.941415      -5.910826       2.052174               1
  -2.040599      -2.032355       2.028040               1
  -2.129653      -6.096888       2.008118               2
%endblock AtomicCoordinatesAndAtomicSpecies

LiF能带结构计算的困惑
Graph1.png


LiF能带结构计算的困惑-1
Graph2.png
回复此楼

» 猜你喜欢

没有明天,只有今天!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沐汐1998

铜虫 (正式写手)

用MS还是VASP?
2楼2015-08-14 16:36:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
franch: 金币+2, 谢谢回帖交流。。。。 2015-08-18 13:59:42
显然看Fermi能级上下两个能带,而且LDA低估带隙也不稀奇。

需要注意的是能带计算k最好充分点,印象path好像也不是这样的。
不一定挂在论坛,计算问题问题欢迎留言。
3楼2015-08-14 22:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ewan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +11 tcxiaoxx 2026-03-20 13/650 2026-03-25 15:38 by barlinike
[考研] 0854电子信息求调剂 +7 α____ 2026-03-22 9/450 2026-03-25 13:37 by α____
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +6 niko- 2026-03-22 6/300 2026-03-25 09:14 by aa331100
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
信息提示
请填处理意见