| 查看: 2467 | 回复: 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
» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流，， 2013-04-22 21:06:57
天气就是风景: 金币+1, ★★★很有帮助, 谢谢 2013-05-02 21:53:48
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流，， 2013-04-22 21:06:57
天气就是风景: 金币+1, ★★★很有帮助, 谢谢 2013-05-02 21:53:48
|
这段话的意思是rprim默认需要写10位小数。LZ抄的时候只抄了8位，就出错了。要用rpim的话就要严格按照手册的来。 同时楼上已经指出了，rprim和angdeg是重复的，设一个就可以了。 另外我不知道这个晶体应该长什么样，但感觉xred也是错的，很明显这样写的xred有z方向总是有0.5c的范围内有4个原子，另外0.5c没有原子，SiC不应该具有这种分离的层状结构，可能是另一半对称的在来源中没有标明。 |
5楼2013-04-18 22:03:00
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
天气就是风景
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2188.2
- 散金: 200
- 红花: 3
- 帖子: 124
- 在线: 41.8小时
- 虫号: 1428904
- 注册: 2011-10-06
- 专业: 光学
|
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
200