24小时热门版块排行榜    

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

spring5908

金虫 (小有名气)

[求助] 求助:带电缺陷HSE06总能计算结果匪夷所思,请赐教!已有1人参与

我最近在尝试用HSE06泛函做一些缺陷及掺杂结构的理论计算,结果得到了一些匪夷所思的数据,至今未找到原因,着急得不行,所以把计算结果和计算参数贴出来,请高手指点!

这个工作是关于带有不同电荷数的缺陷态的理论计算,由于超胞80个原子,用HSE06构型优化太耗时了,所以我用一般的DFT方法进行了全优化,然后用HSE06进行静态计算以得到HSE06方法的总能量。结果,用HSE06泛函得到的不同电荷态的缺陷构型的总能较一般DFT方法差别很大,主要体现在两个方面:一是带电缺陷总能比中性缺陷总能小100多ev,二是不同电荷态之间的总能差变的非常小,几乎可以忽略了!

在做超胞计算之前,我测试了用HSE06泛函计算了单胞的能带结构,结果是合理的,可见泛函本身的使用并没有问题。我将不同电荷缺陷的总能计算结果比较,以及HSE06泛函进行静态计算的INCAR文件贴在下面,请高手帮忙分析下到底是什么错误导致的,如能赐教,不胜感激啊!!!

q=0: -625.5eV (DFT) and -766.2eV (HSE06)
q=-1: -624.2eV (DFT) and -885.1eV (HSE06)
q=-2: -619.6eV (DFT) and -885.2eV (HSE06)
q=-3: -615.0eV (DFT) and -885.3eV (HSE06)

INCAR (HSE06, vasp5.2):
PREC=Accurate
NWRITE=1
ICHARG=2
ISTART=0
INIWAV=1
LREAL=.TRUE.
NELECT=489
EDIFF=1.E-7
NELM=1000
NELMIN=4
NELMDL=-5
ISMEAR=0
SIGMA=0.1
ISPIN=1
LHFCALC=.TRUE.
HFSCREEN=0.2
ALGO=Damped
TIME=0.4
PRECFOCK=Fast
SYMPREC=1E-6
NGXF=144
NGYF=144
NGZF=144
NSW=0
IBRION=-1
ADDGRID=.TRUE.
EDIFFG=-0.05
ISIF=0
POTIM=0.2
GGA=PE
LSCALAPACK=.TRUE.
LPLANE=.TRUE.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring5908

金虫 (小有名气)

引用回帖:
7楼: Originally posted by stractor at 2012-10-24 22:15:16
“一是带电缺陷总能比中性缺陷总能小100多ev”
因为HSE与标准GGA或者LDA中交换关联项不一样,所以总能不一样。再说能量绝对值没有任何意义,只有能量只差才有意义。

同样是HSE,带电缺陷比中性缺陷总能小100多ev,而在DFT中,二者的差值只有几个ev。
8楼2012-10-26 01:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
spring5908: 金币+1, 有帮助 2012-10-23 21:30:51
建议你这样试试:先DFT计算,HSE计算这步读DFT这步的波函数,祝你好运
2楼2012-10-23 21:19:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring5908

金虫 (小有名气)

引用回帖:
2楼: Originally posted by liliangfang at 2012-10-23 21:19:55
建议你这样试试:先DFT计算,HSE计算这步读DFT这步的波函数,祝你好运

谢谢,但是我已经测试过了,计算结果几乎是一样的。
3楼2012-10-23 21:30:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

引用回帖:
3楼: Originally posted by spring5908 at 2012-10-23 21:30:23
谢谢,但是我已经测试过了,计算结果几乎是一样的。...

那就麻烦了,调试其他参数试试吧,用pbe的势,HSE计算要在DFT的基础上。还有一点就是你试试用力收敛
4楼2012-10-23 21:40:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见