24小时热门版块排行榜    

查看: 3116  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wangyan6930

木虫 (著名写手)

[求助] HSE计算的一些问题,希望大家帮我解决谢谢 已有4人参与

经过最近的多次测试,兄弟想请教如下几个问题,我也是做了多次测试很出的结果,已经成功的计算了Si的能带,并和实验值符合的很好,但是测试新的例子的时候出下了以下几个问题,请大家帮忙:
1、关于布里渊区K点的选择,我现在做的测试是NaNbO3的能带计算,空间群是221PM-3M对称性,高对称点是X(0.5,0,0) R(0.5,0.5,0.5) M (0.5,0.5,0),G(0,0,0),我选择的路径是X-R-M-G-R(CASTEP默认路径)和R-X-G-M-R-G分别进行了计算,想请教一下KPOINTS的路径怎么样选取更加合理?
2、我HSE计算能带的文件里面我同时处理了态密度,发现从DOS上读取的带隙和能带图上读取的带隙有一些差距,这个是由于DOS计算的时候加展宽的原因吗?还是同一个计算结果里面同时选择的DOS和能带本来就不合适?也就是HSE的DOS和能带应该分开来算?3
HSE计算出来的带隙约为2.5左边,跟实验值的3.19有些差距,我选择了我的两次不同计算过程的文件给您发过去,请您帮我查看一下错误可以吗?积分路径是RXGMRG 计算出来的能带结果是2.5左右,能不能麻烦帮我看下问题出在什么地方,附件是我选的两次不同计算路径的参数,麻烦帮我看一下,谢谢!

我把我计算的参数贴出来请大家给我指出问题,积分路径R-X-G-M-R-G,主要问题是计算出来的能带和实验值有差距,反复测试了数十次和查资料实在找不到问题,请大家帮忙,谢谢!
PBE计算

INCAR

ENCUT = 520
EDIFF = 1e-5
EDIFFG = 2e-4
ISMEAR=0
SIGMA=0.1
IBRION=-1
NSW=0
ISIF=2
POTIM=0.1
LWAVE=T
LCHARG=T
#NPAR=1

KPOINTS
Automatic mesh
0
Monkhorst Pack
  5   5   5
  0.  0.  0.

POSCAR
Na Nb O3  
    1.00000000000000
     4.0096004013579449    0.0000000011715256    0.0000000011715255
     0.0000000723430511    4.0096004013579449    0.0000000011715256
     0.0000000723430511    0.0000000723430511    4.0096004013579449
   Na   Nb   O
     1     1     3
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.5000000000000000  0.5000000000000000  0.5000000000000000
  0.0000000000000000  0.5000000000000000  0.5000000000000000
  0.5000000000000000  0.5000000000000000 -0.0000000000000000
  0.5000000000000000  0.0000000000000000  0.5000000000000000

HSE自洽计算
INCAR

ENCUT = 520
EDIFF = 1e-5
EDIFFG = 2e-4
ISMEAR=0
SIGMA=0.1
IBRION=-1
NSW=0
ISIF=2
POTIM=0.1
LWAVE=T
LCHARG=T
#NPAR=1
LHFCALC=.TRUE.
HFSCREEN=0.2
ALGO = D
TIME = 0.4

KPOINTS
Automatic mesh
0
Monkhorst Pack
  5   5   5
  0.  0.  0.

HSE能带计算
INCAR
ENCUT = 520
EDIFF = 1e-5
EDIFFG = 2e-4
ISMEAR=0
SIGMA=0.1
IBRION=-1
NSW=0
ISIF=2
POTIM=0.1
LWAVE=T
LCHARG=T
#NPAR=1
LHFCALC=.TRUE.
HFSCREEN=0.2

KPOINTS

Automatically generated mesh
     60
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.20000000000000    0.00000000000000    0.00000000000000             6
    0.40000000000000    0.00000000000000    0.00000000000000             6
    0.20000000000000    0.20000000000000    0.00000000000000            12
    0.40000000000000    0.20000000000000    0.00000000000000            24
    0.40000000000000    0.40000000000000    0.00000000000000            12
    0.20000000000000    0.20000000000000    0.20000000000000             8
    0.40000000000000    0.20000000000000    0.20000000000000            24
    0.40000000000000    0.40000000000000    0.20000000000000            24
    0.40000000000000    0.40000000000000    0.40000000000000             8
   0.50000000  0.50000000  0.50000000       0
   0.50000000  0.44444444  0.44444444       0
   0.50000000  0.38888889  0.38888889       0
   0.50000000  0.33333333  0.33333333       0
   0.50000000  0.27777778  0.27777778       0
   0.50000000  0.22222222  0.22222222       0
   0.50000000  0.16666667  0.16666667       0
   0.50000000  0.11111111  0.11111111       0
   0.50000000  0.05555556  0.05555556       0
   0.50000000  0.00000000  0.00000000       0
   0.50000000  0.00000000  0.00000000       0
   0.44444444  0.00000000  0.00000000       0
   0.38888889  0.00000000  0.00000000       0
   0.33333333  0.00000000  0.00000000       0
   0.27777778  0.00000000  0.00000000       0
   0.22222222  0.00000000  0.00000000       0
   0.16666667  0.00000000  0.00000000       0
   0.11111111  0.00000000  0.00000000       0
   0.05555556  0.00000000  0.00000000       0
   0.00000000  0.00000000  0.00000000       0
   0.00000000  0.00000000  0.00000000       0
   0.05555556  0.05555556  0.00000000       0
   0.11111111  0.11111111  0.00000000       0
   0.16666667  0.16666667  0.00000000       0
   0.22222222  0.22222222  0.00000000       0
   0.27777778  0.27777778  0.00000000       0
   0.33333333  0.33333333  0.00000000       0
   0.38888889  0.38888889  0.00000000       0
   0.44444444  0.44444444  0.00000000       0
   0.50000000  0.50000000  0.00000000       0
   0.50000000  0.50000000  0.00000000       0
   0.50000000  0.50000000  0.05555556       0
   0.50000000  0.50000000  0.11111111       0
   0.50000000  0.50000000  0.16666667       0
   0.50000000  0.50000000  0.22222222       0
   0.50000000  0.50000000  0.27777778       0
   0.50000000  0.50000000  0.33333333       0
   0.50000000  0.50000000  0.38888889       0
   0.50000000  0.50000000  0.44444444       0
   0.50000000  0.50000000  0.50000000       0
   0.50000000  0.50000000  0.50000000       0
   0.44444444  0.44444444  0.44444444       0
   0.38888889  0.38888889  0.38888889       0
   0.33333333  0.33333333  0.33333333       0
   0.27777778  0.27777778  0.27777778       0
   0.22222222  0.22222222  0.22222222       0
   0.16666667  0.16666667  0.16666667       0
   0.11111111  0.11111111  0.11111111       0
   0.05555556  0.05555556  0.05555556       0
   0.00000000  0.00000000  0.00000000       0
回复此楼

» 收录本帖的淘帖专辑推荐

心理互助

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

wangyan6930

木虫 (著名写手)

引用回帖:
3楼: Originally posted by lyd924 at 2015-01-07 11:47:14
能隙应该主要从能带来判断吧,态密度展宽与smear有关,会造成能隙与能带图中的有点差异

楼主你好,我就是利用脚本程序处理出来的能带来看带隙的,不是从DOS上看的,看出来还是有点小,能不能从参数上分析一下我的有什问题。
4楼2015-01-07 14:39:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

chuanghua304

禁虫 (职业作家)

★ ★ ★ ★ ★
感谢参与,应助指数 +1
wangyan6930: 金币+5, 有帮助, 楼上有没有用VASP做HSE计算能带,我就是想讨论一下我碰到的问题,我的KPOINTS是4*4*4的,确实不是很;晶格常数是计算所得,跟实验值接近;ISMEAR=-5对能带计算有没有影响?HSE计算好像大家使用的ISMEAR=0 来做HSE计算的,我可以和楼主交流一下不? 2015-01-07 09:54:16
本帖内容被屏蔽

2楼2015-01-07 08:29:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lyd924

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2015-01-07 21:23:39
wangyan6930: 金币+5, 有帮助 2015-01-21 01:13:53
能隙应该主要从能带来判断吧,态密度展宽与smear有关,会造成能隙与能带图中的有点差异
3楼2015-01-07 11:47:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)

引用回帖:
4楼: Originally posted by wangyan6930 at 2015-01-07 14:39:36
楼主你好,我就是利用脚本程序处理出来的能带来看带隙的,不是从DOS上看的,看出来还是有点小,能不能从参数上分析一下我的有什问题。...

在朋友的指点下 我发现了自己一个参数设置的错 不知道对结果有没有影响EDIFFG=2e-4设置错了应该设置成EDIFFG=-2e-4,但是不知道对计算结果有没有影响,我会继续做测试,也请各位高手继续指出我的错误。
5楼2015-01-07 15:41:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见