| 查看: 1960 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
abinit 计算4H-SiC问题（六角晶系）
|
|||
|
我想对4H-SiC(六角晶系)的总能量问题进行计算，但是abinit始终不能启动，出现如下问题： symbrav.F90:229:BUG Could not succeed to determine the bravais lattice problem,iaxis,invariant= 1 1 1 bravais(1)= 3 iholohedry= 1 action : contact ABINIT group. 首先我个人查到了相关的文献支持，abinit官方没有说不支持六角晶系的计算，同样我还查到早在2009年沈阳工业大学的女生已经发表用abinit计算ZnO薄膜的文章和硕士论文。所以我可以肯定的说六角晶系肯定能计算，但就是我找不到我的输入文件到底出现了什么问题，我试过修改xred，但是还是出现同样的问题。不是以前遇到的问题，显示是error，这个只是contact ABINIT group。到底是什么问题，求高人指点。（附输入文件和2009沈阳工业大学硕士论文） # Crystalline silicon : computation of the total energy #Definition of the unit cell acell 2*5.8128 9.5393 # This is equivalent to 10.18 10.18 10.18 angdeg 90 90 120 rprim 0.866025403 0.5 0 -0.866025403 0.5 0 0.0 0.0 1.0 #Definition of the atom types ntypat 2 # There is only one type of atom znucl 6 14 # 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 natom 4 # There are two atoms typat 1 1 2 2 xred 0 0 0 0.3333 -0.3333 0.25 0 0 0.188 0.3333 -0.3333 0.438 #xred 同样可以写为 # 1/3 2/3 0 #2/3 1/3 0 #1/3 2/3 0.385 #2/3 1/3 0.885 会出现同样的错误。 #Definition of the planewave basis set ecut 20.0 # Maximal kinetic energy cut-off, in Hartree #Definition of the k-point grid kptopt 1 # Option for the automatic generation of k points, taking # into account the symmetry ngkpt 8 8 8 # This is a 2x2x2 grid based on the primitive vectors nshiftk 1 # of the reciprocal space (that form a BCC lattice !), # repeated four times, with different shifts : shiftk 0.0 0.0 0.5 # In cartesian coordinates, this grid is simple cubic, and # actually corresponds to the # so-called 4x4x4 Monkhorst-Pack grid #Definition of the SCF procedure nstep 10 # Maximal number of SCF cycles toldfe 1.0d-6 # 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 diemac 12.0 # Although this is not mandatory, it is worth to # precondition the SCF cycle. The model dielectric # function used as the standard preconditioner # is described in the "dielng" input variable section. # Here, we follow the prescription for bulk silicon. # add to conserve old < 6.7.2 behavior for calculating forces at each SCF step optforces 1 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ZnO薄膜光学性质及其第一原理计算研究.caa
2013-04-18 14:08:45, 376 bytes
» 猜你喜欢
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有9人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ABINIT和VASP相比，计算GW，哪个更快啊？
已经有13人回复
- Abinit计算Ni的磁性，太小了？
已经有14人回复
- abinit计算晶体生成焓的问题讨论
已经有7人回复
- Abinit计算声子问题
已经有10人回复
- 求助abinit计算电光系数
已经有7人回复
- 6H-SiC, 4H-SiC和3C-SiC中的6H、4H、3C代表什么意思
已经有7人回复
- 关于ABINIT计算四方结构弹性常数不收敛的问题
已经有5人回复
- Abinit 新旧版计算结果差别多大可以认为正常，没有安装错误之类的问题。
已经有4人回复
- 【求助】abinit声子计算
已经有7人回复
- 【求助】怎么确定abinit计算能带结构的费米能位置
已经有6人回复
- 【求助】请问各位高手ABINIT输出的Pressure 代表什么压力 和做压力计算有关吗
已经有9人回复
- 【求助】求助abinit和vasp计算的总能的差异
已经有6人回复
- 【求助】Abinit中计算PDOS时怎样确定投影原子半径？
已经有4人回复
- 【求助】Abinit能计算块体材料电导率吗？
已经有6人回复
|
感谢你的建议，我看到有MS建立abinit输入文件的文章，所以我用这个方法做的，所以能否给看下正确与否。 xred 0.00000 0.00000 0.18800 0.00000 0.00000 0.00000 0.66667 0.33333 0.43900 0.66667 0.33333 0.25080 参考文件是MS导出的cif文件 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : SiC(2).c
2013-04-19 10:51:34, 1.65 K
9楼2013-04-19 10:51:44
WDD880227
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 25 (小学生)
- 贵宾: 0.098
- 金币: 1348.5
- 散金: 194
- 红花: 16
- 帖子: 761
- 在线: 371.6小时
- 虫号: 1129580
- 注册: 2010-10-22
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
2楼2013-04-18 16:23:28
|
abinit 对rprim的定义很明确，是这样的。 Alternatively to rprim, directions of dimensionless primitive vectors can be specified by using the input variable angdeg. This is especially useful for hexagonal lattices (with 120 or 60 degrees angles). Indeed, in order for symmetries to be recognized, rprim must be symmetric up to tolsym (10 digits by default), inducing a specification such as rprim 0.86602540378 0.5 0.0 -0.86602540378 0.5 0.0 0.0 0.0 1.0 that can be avoided thanks to angdeg: angdeg 90 90 120 所以我觉得应该没问题，具体参考 http://www.abinit.org/documentat ... s/varbas.html#rprim |
3楼2013-04-18 16:42:14
4楼2013-04-18 21:09:34