24小时热门版块排行榜    

查看: 1370  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +8 妮妮ninicgb 2026-03-15 12/600 2026-03-19 05:55 by anny19840123
[考研] 321求调剂 +3 何润采123 2026-03-18 3/150 2026-03-18 21:27 by li123456789.
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
请填处理意见