| 查看: 2254 | 回复: 4 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
忽悠小皮铜虫 (小有名气)
|
[交流]
【求助成功】石墨烯能带计算的问题 已有4人参与
|
|||
|
大家好。我是siesta新手，最近学习计算石墨烯的能带和态密度。我的计算结果很是怪异。请大家帮忙看看我的输入文件有什么问题。谢谢！ SystemName graphene band structure SystemLabel g-b NumberOfAtoms 2 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 6 C %endblock ChemicalSpeciesLabel MeshCutoff 300. Ry SolutionMethod diagon DM.UseSaveDM .true. UseSaveData .true. SpinPolarized false xc.functional LDA xc.authors CA SaveDeltaRho .true. SaveRho .true. SaveTotalPotential .true. WriteCoorStep .true. %block LocalDensityOfStates -8.5 1.50 eV %endblock LocalDensityOfStates %block ProjectedDensityOfStates -8.5 1.5 0.1 800 eV %endblock ProjectedDensityOfStates BandLinesScale ReciprocalLatticeVectors WriteBands .true. WriteKbands .true. %block BandLines 1 0.333333 0.333333 0.0000 K 100 0.000000 0.000000 0.0000 \Gamma 100 0.000000 0.500000 0.000 M 100 0.333333 0.333333 0.0000 K %endblock BandLines LatticeConstant 2.46 Ang %block LatticeParameters 1.0 1.0 10.0 90. 90. 60. %endblock LatticeParameters AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.333333 0.333333 0.0000 1 0.666666 0.666666 0.0000 1 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 30 0 0 0.0 0 30 0 0.0 0 0 1 0.0 %endblock kgrid_Monkhorst_Pack [ Last edited by 忽悠小皮 on 2010-5-20 at 17:29 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 SIESTA |
» 猜你喜欢
收到复试调剂但是去不了
已经有8人回复
-
260求调剂
已经有5人回复
-
急需调剂
已经有10人回复
-
申博/考博
已经有4人回复
-
化工学硕294分，求导师收留
已经有37人回复
-
求调剂
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有19人回复
-
304求调剂
已经有7人回复
-
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关石墨烯能带结构的计算
已经有11人回复
- 诚求CASTEP如何计算石墨稀的能带图
已经有16人回复
- 双层石墨烯的能带、pdos、以及差分电荷密度
已经有5人回复
- dmol相关问题交流学习
已经有14人回复
- 速求 用matlab写 armchair型石墨烯纳米条带 的能带图
已经有8人回复
- 关于MS中通过 Brillouin Zone Path 设置和修改的问题
已经有31人回复
- Raman光谱可以用来判断石墨烯的层数吗？用化学还原法合成的石墨烯。
已经有20人回复
- 石墨烯中为什么有两个不等价的原子
已经有32人回复
- 【求助】石墨烯的MS建模
已经有13人回复
- 【讨论】从石墨烯得奖想出去
已经有15人回复
- [求助]什么结构可以在VNL 中计算能带和DOS
已经有16人回复
- 【重点讨论】计算能带时如何确定高对称点和k-points path?
已经有40人回复
- 【交流】关于求解石墨烯的能带矩阵
已经有8人回复
- 【求助】初学者求助： 紧束缚方法计算能带结构
已经有5人回复
xhzha
铁杆木虫 (正式写手)
- 应助: 37 (小学生)
- 金币: 9454.6
- 散金: 30
- 红花: 7
- 帖子: 409
- 在线: 577.7小时
- 虫号: 1101711
- 注册: 2010-09-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2011-12-16 13:06:04
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指导 2010-05-20 08:25:46
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指导 2010-05-20 08:25:46
|
可能基矢弄错了。 LatticeConstant 2.46 Ang %block LatticeParameters 1.0 1.0 10.0 90. 90. 60. %endblock LatticeParameters 需改为： %block LatticeParameters 2.46 2.46 10.0 90. 90. 120.0 %endblock LatticeParameters 或者： LatticeConstant 2.46 Ang %block LatticeVectors 1.0 0.0 0.0 -0.5 0.866025 0.0 0.0 0.0 4.06504 %endblock LatticeVectors 你也需要检查原子坐标。 %block AtomicCoordinatesAndAtomicSpecies 0.0000000000000000 0.0000000000000000 0.0000000000000000 1 0.3333333333333333 0.6666666666666667 0.0000000000000000 1 %endblock AtomicCoordinatesAndAtomicSpecies [ Last edited by valenhou001 on 2010-5-19 at 21:53 ] |
2楼2010-05-19 21:51:52
忽悠小皮
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 134.7
- 散金: 168
- 帖子: 79
- 在线: 31.9小时
- 虫号: 688416
- 注册: 2009-01-06
- 专业: 凝聚态物性 II ：电子结构
3楼2010-05-20 15:03:14
jixiao0526
金虫 (小有名气)
- 应助: 12 (小学生)
- 金币: 808.3
- 散金: 158
- 红花: 2
- 帖子: 151
- 在线: 235.6小时
- 虫号: 1407476
- 注册: 2011-09-19
- 性别: GG
- 专业: 电磁场与波
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+0.5):给个红包，谢谢回帖
|
初学计算，计算的结果中没有.band文件，哪位高手帮我看看问题出在哪 NumberOfAtoms 9 NumberOfSpecies 2 %block ChemicalSpeciesLabel 1 6 C 2 3 Li %endblock ChemicalSpeciesLabel XC.functional LDA XC.authors CA # ===================================================================================== # Cell & K-MP-grid LatticeConstant 4.9200 Ang # Lattice constant alat %block LatticeVectors # in units of alat 1.0 0.0000 0.0 -0.5 0.866 0.00 0.0 0.0 3.048780488 # c=15.08000A %endblock LatticeVectors #Basis set generation PAO.BasisSize DZP PAO.BasisType split PAO.EnergyShift 100 meV #Molecular Dynamic Calculation #MD #MD.TypeOfRun Broyden #MD.NumCGsteps 100 #MD.MaxForceTol 0.05 eV/Ang Functional and Solution Method(Order-N/diagoalization) SpinPolarized false SolutionMethod diagon MeshCutoff 150.0 Ry %block kgrid_Monkhorst_Pack 8 0 0 0.0 0 8 0 0.0 0 0 1 0.5 %endblock kgrid_Monkhorst_Pack # ===================================================================================== SCF options MaxSCFIterations 100 # Maximum number of SCF iter DM.NumberPulay 3 # One Pulay every 5 iter DM.MixingWeight 0.25 # New DM amount for next SCF cycle DM.Tolerance 1.d-4 # Tolerance in maximum difference DM.UseSaveDM true # to use continuae..ion files SolutionMethod diagon # OrderN or Diagon DM.UseSaveDM true # Restart DM.AllowReuse T DM.AllowExtrapolation F MeshCutoff 150. Ry # OK for carbon (see Emilio) OccupationFunction FD ElectronicTemperature 250 meV # Temp. for Fermi smearing # ===================================================================================== %block BandLine 1 0.0000000000 0.0000000000 0.0000000000 \GAMMA 30 -0.333 0.667 0 F 20 0.0 0.5 0 Q 20 0 0 0.0 Z %endblock BandLines # ============================================================================ %block ProjectedDensityOfStates -10.0 10.0 0.05 100 eV %endblock ProjectedDensityOfStates # ===================================================================================== # Output Options LongOutput true BandLinesScale ReciprocalLatticeVectors # ===================================================================================== AtomicCoordinatesFormat Ang # Format for inates %block AtomicCoordinatesAndAtomicSpecies < aa.xyz |
4楼2011-12-10 13:29:10
