| 查看: 1308 | 回复: 4 | |||
hisap-nano新虫 (正式写手)
|
[交流]
为什么算能带时老是出错！ 已有2人参与
|
|
在算简单的金属的能带时老是提示出错，我看了半天也没整明白错在哪里了，大家帮我看看好吗，先谢了！ ndtset 2 #Optimization of the lattice parameters optcell1 1 ionmov1 3 ntime1 10 dilatmx1 1.05 ecutsm1 0.5 #Definition of the k-point grids kptopt1 1 # Option for the automatic generation of k points, taking # into account the symmetry nshiftk1 1 # These shifts will be the same for all hcp grids shiftk1 0.0 0.0 0.5 ngkpt1 6 6 6 prtden1 1 #Dataset 2 : the band structure iscf2 -2 getden2 -1 kptopt2 -3 nband2 24 ndivk2 15 15 15 15 15 kptbounds2 0.0 0.0 0.5 # A point 0.0 0.0 0.0 # Gamma point 0.5 0.0 0.0 # M point 2/3 1/3 0.0 # K point 2/3 1/3 0.5 # H point 1.0 1.0 1.0 # Gamma point in another cell. tolwfr2 1.0d-12 enunit2 1 # Will output the eigenenergies in eV prtdos2 2 prtdos2 3 #Definition of the unit cell acell 8.4776 8.4776 13.2906 #in bohr atomic unit angdeg 90.0000 90.0000 120.0000 #Definition of the atom types ntypat 1 # There is only one type of atom znucl 4 # The keyword "znucl" refers to the atomic number of the # possible type(s) of atom. The pseudopotential(s) # mentioned in the "files" file must correspond # to the type(s) of atom. Here, the only type is Silicon. #Definition of the atoms chkprim 0 natom 16 # There are 16 atoms typat 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 xred 0.16667 0.33333 0.12500 0.33333 0.16667 0.37500 0.66667 0.33333 0.12500 0.83333 0.16667 0.37500 0.16667 0.83333 0.12500 0.33333 0.66667 0.37500 0.66667 0.83333 0.12500 0.83333 0.66667 0.37500 0.16667 0.33333 0.62500 0.33333 0.16667 0.87500 0.66667 0.33333 0.62500 0.83333 0.16667 0.87500 0.16667 0.83333 0.62500 0.33333 0.66667 0.87500 0.66667 0.83333 0.62500 0.83333 0.66667 0.87500 #Definition of the planewave basis set ecut 10.0 # Maximal kinetic energy cut-off, in Hartree #Definition of the SCF procedure nstep 15 # Maximal number of SCF cycles toldfe 1.0d-8 # Will stop when, twice in a row, the difference # between two consecutive evaluations of total energy # differ by less than toldfe (in Hartree) # This value is way too large for most realistic studies of materials |
回复此楼» 猜你喜欢
有爆料，一个青年教师卖房得400万，然后换了一个四青帽子
已经有10人回复
-
一志愿西交机械专硕求调剂
已经有9人回复
-
245求调剂
已经有6人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有15人回复
-
材料复试求调剂
已经有24人回复
-
327求调剂
已经有11人回复
-
288求调剂，一志愿华南理工大学071005
已经有20人回复
-
有没有学校材料专业收跨调(一志愿085410)
已经有16人回复
-
各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。
已经有5人回复
-
化工学硕294分，求导师收留
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于使用MS中的CASTEP如何计算能带的问题
已经有9人回复
- 【求助】理论计算新手，求帮忙解决带隙及DOS的问题
已经有21人回复
- 【求助】小弟ATK计算的文章，求助审稿人关于能带的一个问题！江湖告急！ 望大家多帮忙
已经有9人回复
- 【求助】关于MS软件中能带和态密度的计算问题
已经有8人回复
- 【求助】VASP计算能带时KPOINTS设置出错
已经有14人回复
- 【求助成功】VASP算出来的能带和DOS图对应不上，哪里出了问题？
已经有20人回复
- 【求助】castep计算能带出错？
已经有3人回复
- 【求助成功】计算能带结构时出现错误，请教
已经有19人回复
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
2楼2011-12-17 09:27:29
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
3楼2011-12-17 11:29:54
hisap-nano
新虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 924.1
- 红花: 2
- 帖子: 412
- 在线: 226.3小时
- 虫号: 297318
- 注册: 2006-11-18
- 专业: 金属材料的力学行为
4楼2011-12-17 16:11:58
hisap-nano
新虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 924.1
- 红花: 2
- 帖子: 412
- 在线: 226.3小时
- 虫号: 297318
- 注册: 2006-11-18
- 专业: 金属材料的力学行为
hisap-nano: 回帖置顶 2011-12-18 12:21:05
hisap-nano: 取消置顶 2011-12-19 04:39:17
hisap-nano: 取消置顶 2011-12-19 04:39:17
|
斑竹，真对不起，还要向您请教，我修改了input文件后，程序能出结果，我用abinitbandstructuremaker 去画能带，结果提示算能带的数组不完整，您帮我看看是怎么回事好吗，先谢了！ input file ndtset 2 #Optimization of the lattice parameters optcell1 1 ionmov1 3 ntime1 10 dilatmx1 1.05 ecutsm1 0.5 iscf1 5 nstep1 10 toldfe1 1.0d-8 #Definition of the k-point grids kptopt1 1 nshiftk1 1 shiftk1 0.0 0.0 0.5 ngkpt1 6 6 6 prtden1 1 #Dataset 2 : the band structure iscf2 -2 nstep2 50 getden2 -1 kptopt2 -5 # the absolut value represent the segments ndivk2 15 15 15 15 15 # 15, 15, 15, 15 and 15 divisions of the 5 segments, delimited by 6 points. kptbounds2 0.0 0.0 0.5 0.0 0.0 0.0 0.5 0.0 0.0 1/3 1/3 0.0 1/3 1/3 0.5 1.0 1.0 1.0 tolwfr2 1.0d-12 enunit2 1 #Definition of the unit cell acell 8.4776 8.4776 13.2906 rprim 1.0 0.0 0.0 -0.5 sqrt(0.75) 0.0 0.0 0.0 1.0 #Definition of the atom types ntypat 1 znucl 4 #Definition of the atoms chkprim 0 natom 16 # There are 16 atoms typat 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 xred 0.16667 0.33333 0.12500 0.33333 0.16667 0.37500 0.66667 0.33333 0.12500 0.83333 0.16667 0.37500 0.16667 0.83333 0.12500 0.33333 0.66667 0.37500 0.66667 0.83333 0.12500 0.83333 0.66667 0.37500 0.16667 0.33333 0.62500 0.33333 0.16667 0.87500 0.66667 0.33333 0.62500 0.83333 0.16667 0.87500 0.16667 0.83333 0.62500 0.33333 0.66667 0.87500 0.66667 0.83333 0.62500 0.83333 0.66667 0.87500 #Definition of the planewave basis set ecut 10.0 错误提示 2 dataset(s) detected the last dataset is not completed and will be ignored Traceback (most recent call last): File "AbinitBandStructureMaker.py", line 1012, in ? dico[k] = '{l%s}' %ctrl.alphabet[L] KeyError: 27 |
5楼2011-12-18 12:20:28
