24小时热门版块排行榜    

查看: 1365  |  回复: 2

lemonare

新虫 (小有名气)

[求助] Siesta能带计算,带隙总是不对

我是Siesta新手,现在在尝试计算ZnS的能带结构,计算出的结果看起来形状还可以,但在带隙上差的很大,希望高手帮忙看看,给予修改建议~
FDF文件:

SystemName          ZnS
SystemLabel         ZnS

NumberOfAtoms       2
NumberOfSpecies     2

%block ChemicalSpeciesLabel
1  30  Zn
2  16   S
%endblock ChemicalSpeciesLabel

PAO.BasisSize       DZP
PAO.EnergyShift     50 meV

AtomicCoordinatesFormat ScaledCartesian # Format for coordinates
AtomicCoorFormatOut     Ang

%block AtomicCoordinatesAndAtomicSpecies
   .000   .000   .000   1
   .250   .250   .250   2
%endblock AtomicCoordinatesAndAtomicSpecies

LatticeConstant     5.409 Ang
%block LatticeVectors
  0.000  0.500  0.500
  0.500  0.000  0.500
  0.500  0.500  0.000
%endblock LatticeVectors

kgrid_cutoff        20. Ang

%block BandLines
1 0.5000 0.2500 0.7500
20 0.5000 0.5000 0.5000
20 0.0000 0.0000 0.0000
20 0.5000 0.0000 0.5000
20 0.5000 0.2500 0.7500
20 0.3750 0.3750 0.7500
%endblock BandLines

%block ProjectedDensityOfStates
-15.00 20.00 0.05 500 eV
%endblock ProjectedDensityOfStates

# DFT, Grid, SCF

XC.functional           LDA         # Exchange-correlation functional type
XC.authors              CA          # Particular parametrization of xc func
SpinPolarized           .false.     # Spin unpolarized calculation
MeshCutoff              300. Ry     # Equivalent planewave cutoff for the grid
MaxSCFIterations        150         # Maximum number of SCF iterations per step
DM.MixingWeight         0.3         # New DM amount for next SCF cycle
DM.Tolerance            1.d-4       # Tolerance in maximum difference
                                    # between input and output DM
DM.NumberPulay          5           # Number of SCF steps between pulay mixing
         

SolutionMethod       diagon   
ElectronicTemperature  25 meV  

MD.TypeOfRun         CG
MD.NumCGsteps        300
MD.MaxCGDispl         0.05  Ang
MD.MaxForceTol        0.04 eV/Ang

MD.VariableCell        .true.

WriteCoorInitial     
WriteCoorStep      
WriteForces         
WriteKpoints            .true.
WriteEigenvalues        .true.
WriteKbands             .true.
WriteBands              .true.
WriteMullikenPop        1            # Write Mulliken Population Analysis
WriteCoorXmol           .false.
WriteMDCoorXmol         .false.
WriteMDhistory          .false.
WriteCoorXmol           .false.


# Options for saving/reading information

DM.UseSaveDM            .true.      # Use DM Continuation files
MD.UseSaveXV            .true.      # Use stored positions and velocities
MD.UseSaveCG            .true.      # Use stored positions and velocities
SaveRho                              # Write valence pseudocharge at the mesh
SaveDeltaRho                         # Write RHOscf-RHOatm at the mesh
SaveElectrostaticPotential .false.   # Write the total elect. pot. at the mesh
SaveTotalPotential      .false.      # Write the total pot. at the mesh
WriteSiestaDim          .false.      # Write minimum dim to siesta.h and stop
WriteDenchar            .true.       # Write information for DENCHAR
WriteWaveFunctions      .true.
UseSavedata             .true.
回复此楼

» 收录本帖的淘帖专辑推荐

siesta

» 猜你喜欢

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

每一天都是一个成长的过程!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)

【答案】应助回帖

★ ★
zzy870720z(金币+2): 鼓励交流 2011-10-15 09:26:31
lemonare(金币+2): 不是的,后来我找到问题了 2011-10-17 13:30:13
估计和赝势有关系吧。严格的说你的赝势估计没有做过验证,不够准确。我最近也发现自己用siesta算的能隙有点问题,但是用别的软件就算的好。结构参数都差不多,关键siesta的赝势既不是官方的,也没有经过验证,麻烦。
坠入无边的网~~
2楼2011-10-14 15:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemonare

新虫 (小有名气)

是这一行:
AtomicCoordinatesFormat Fractional
每一天都是一个成长的过程!
3楼2011-10-17 13:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lemonare 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 285化工学硕求调剂(081700) +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
信息提示
请填处理意见