| 查看: 1377 | 回复: 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 |
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
有关石墨烯能带结构的计算
已经有11人回复
- 半导体价带顶和能带间隙计算？
已经有10人回复
- 最新计算能带的方法 得出的带隙和实验非常的吻合。
已经有9人回复
- SIESTA中能带文件*.bands的处理
已经有49人回复
- siesta计算dipole moment（电偶极矩）
已经有7人回复
- 计算一种氧化物，光吸收系数计算的光吸收边远大于能带带隙，求解释
已经有24人回复
- LDA+U到底能够对能带的计算有多大的改善呢
已经有22人回复
- 【求助成功】abinit 计算CaO能带结构图，间接带隙直接带隙问题
已经有9人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】用siesta计算六角晶格的能带结构（有新信息，欢迎讨论）
已经有4人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【求助】求助解决用mpich2.1 编译siesta，计算出现错误，
已经有12人回复
- 【讨论】利用castep计算ZnO能带带隙 算过的看看
已经有20人回复
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
2楼2011-10-14 15:03:06
lemonare
新虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1907.9
- 散金: 6
- 红花: 1
- 帖子: 52
- 在线: 65.8小时
- 虫号: 1407234
- 注册: 2011-09-19
- 性别: GG
- 专业: 半导体材料
3楼2011-10-17 13:31:27
