24小时热门版块排行榜    

查看: 1367  |  回复: 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的回帖

lemonare

新虫 (小有名气)

是这一行:
AtomicCoordinatesFormat Fractional
每一天都是一个成长的过程!
3楼2011-10-17 13:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601求调剂 +3 Du.11 2026-03-16 3/150 2026-03-16 20:42 by 无际的草原
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见