【求助成功】计算能带结构时出现错误，请教 - 第一性原理 - Vasp&MedeA

11楼2010-08-10 22:42:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):多谢侯老师指教 2010-08-11 09:33:38

你的基矢：
liBH4-Hex
4.27600000000000
   0.8481987248411803 -0.4897079788164736 0.0000000000000000
   0.0000000000000000 0.9794159576329472 0.0000000000000000
   0.0000000000000000 0.0000000000000000 1.7339757780272302
即
a
0.5 * sqrt(3.0)          -0.5       0.000
0.00000000000       1.0       0.000
0.0                         0.0          c/a

这个也是合理的。

你将INCAR中的增加 SYMPREC 关键词，将它的值设置的比默认值大些。

12楼2010-08-11 09:20:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

1ST强帖: 1
应助: 5 (幼儿园)
金币: 401.2
散金: 15
红花: 3
帖子: 264
在线: 87.8小时
虫号: 1004137
注册: 2010-04-23
性别: GG
专业: 原子和分子物理

引用回帖:

Originally posted by valenhou001 at 2010-08-11 09:20:18:
你的基矢：
liBH4-Hex
4.27600000000000
   0.8481987248411803 -0.4897079788164736 0.0000000000000000
   0.0000000000000000 0.9794159576329472    ...

加入SYMPREC=1E-4之后，还是不行。我估计，是不是产生K点的方法不对？？？但是师兄说是原子位置的问题。。。我就是不知道是什么问题？？？

[ Last edited by qiqi2926 on 2010-8-11 at 15:14 ]

13楼2010-08-11 10:41:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

1ST强帖: 1
应助: 5 (幼儿园)
金币: 401.2
散金: 15
红花: 3
帖子: 264
在线: 87.8小时
虫号: 1004137
注册: 2010-04-23
性别: GG
专业: 原子和分子物理

引用回帖:

Originally posted by jghe at 2010-08-10 22:42:19:

是的，但是每一对高对称点间要空两行

我试过这种方法。每一对搞对称断电之间空一行行吗？我看侯老师的书上就是空一行。还是不行。。。这个问题的症结我还不是很清楚。您能讲讲吗？？？

14楼2010-08-11 10:43:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):呵呵，热心人。 2010-08-11 15:54:59

可以的话，你做自洽计算时的INCAR, POSCAR, KPOINTS, OUTCAR， POTCAR,还有做band structure计算时的 INCAR, KPOINTS, syml, 压缩之后发到我的email: valenhou@gmail.com。我来瞧瞧。

还有注明你用的是哪个vasp版本。

15楼2010-08-11 15:42:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢建议 2010-08-12 18:26:30

在计算能带时，不能设置ISMEAR=-5或-4，应该设置其他的，比如ISMEAR =0 或1。

16楼2010-08-11 20:08:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qiqi2926

铜虫 (小有名气)

1ST强帖: 1
应助: 5 (幼儿园)
金币: 401.2
散金: 15
红花: 3
帖子: 264
在线: 87.8小时
虫号: 1004137
注册: 2010-04-23
性别: GG
专业: 原子和分子物理

引用回帖:

Originally posted by valenhou001 at 2010-08-11 20:08:37:
在计算能带时，不能设置ISMEAR=-5或-4，应该设置其他的，比如ISMEAR =0 或1。

我把ISMEAR给去掉了，不设值。。。
其中INCAR文件为：
SYSTEM = LiBH4-Hex
ENCUT = 600
ISTART = 1; ICHARG = 2
#ISMEAR = -5
NBANDS = 20
PREC = Accurate

就可以运行了，但是不知道为什么？？？能不能解释一下？？？

17楼2010-08-11 23:45:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jghe

金虫 (著名写手)

应助: 10 (幼儿园)
贵宾: 0.225
金币: 112.2
散金: 5
红花: 3
帖子: 1481
在线: 123小时
虫号: 322338
注册: 2007-03-11
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 量子化学

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢指导 2010-08-12 18:26:18

就是ISMEAR得问题，去掉就按默认值了，ISMEAR=1 Sigma=0.2

18楼2010-08-12 09:59:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingfengbo

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 79
帖子: 70
在线: 38.4小时
虫号: 939022
注册: 2010-01-08
专业: 凝聚态物性 II ：电子结构

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by qiqi2926 at 2010-08-10 14:05:06:
计算六角晶系的能带结构，自洽得到电荷密度已经完成，再进行第二步，也就是用所得到的电荷密度再重新做非自洽计算时出现问题。
《参考侯老师的文章》：找出简约布里渊区的特殊k点，编辑了syml文件（利用gk.x产生 ...

您好！能不能把您的gk.x程序给我发一份？邮箱gy1357900000@163.com，我用gk.x和syml产生的输出文件KPOINTS和inp.kpt里面什么也没有，我想可能是我的gk.x程序有问题！谢谢！

赞一下(1人)

19楼2010-09-08 18:34:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nanost

木虫 (知名作家)

Dr.

应助: 35 (小学生)
金币: 4696.6
散金: 228
红花: 18
帖子: 5211
在线: 1195.5小时
虫号: 195830
注册: 2006-02-23
性别: GG
专业: 无机合成和制备化学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

19楼: Originally posted by dingfengbo at 2010-09-08 18:34:22
您好！能不能把您的gk.x程序给我发一份？邮箱gy1357900000@163.com，我用gk.x和syml产生的输出文件KPOINTS和inp.kpt里面什么也没有，我想可能是我的gk.x程序有问题！谢谢！

...

我也遇到了相同问题，您解决了吗？谢谢