24小时热门版块排行榜    

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

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

金虫 (小有名气)

引用回帖:
4楼: Originally posted by liliangfang at 2012-10-23 21:40:07
那就麻烦了,调试其他参数试试吧,用pbe的势,HSE计算要在DFT的基础上。还有一点就是你试试用力收敛...

你所说的用力收敛具体是怎么做,可否说的详细一点?
另外,我还有两个问题。之前的DFT必须用pbe么,还是用pw91也可以?另外,http://muchong.com/html/201205/4421729.html这个帖子中提到HSE不能进行静态计算,并且提到总能用DFT结果即可,你认为合理么?
谢谢!
5楼2012-10-23 21:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请科研项目纯靠实力能成功吗? 200+3 David1982 2026-10-13 6/300 2026-10-14 11:14 by 北京莱茵编辑
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见