24小时热门版块排行榜    

查看: 2006  |  回复: 7

wd12309812

铜虫 (正式写手)

[求助] 算NEB不收敛,大家帮忙看一下这个INCAR有什么问题 已有2人参与

NELM = 200
结果OSZICAR都是到200自动停止,d E只达到-2或-3次方(我设置-5),求帮忙看一下INCAR有什么问题!

非常感谢!!!



SYSTEM=LTO7

Startparameter for this run
#NWRITE    =    2     LPETIM=.TURE.        # write-flag&timer
#INIWAV   =    1                          # INIWAV, default is 1,only useful when Istart=0
ISTART    =    0        #job : 0-new  1-constant cut-off  2-constant cell size/shape.and cutoff read from WAVECAR
ICHARG    =    2        #job : 0-from initial wf  1-from CHGCAR   2-atomic charge density , relax/DOS=2  EK=11
LCHARG    =  .F.
LWAVE     =   .F.
ISPIN     =    1        
#ISYM      =    0        #default 1-US-pp  2-PAW  symmetry stuff ON(1 or 2) or OFF(0)
PREC      =   Normal

#GGA       =    PE
#VOSKOWN  =    1        #when GGA91 used for magnetic calculation
#LORBIT    =    11      # for pdos
ENCUT     =   300 eV

Electronic Relaxation
NELMIN    =    4
NELM      =   200       #   the max number of electronic SC
EDIFF     =    1E-05    #   stopping-criterion for ELM
LREAL     =    Auto     #   F-reciprocal  T-real  A-auto

#L(S)DA+U
#LDAU = .TRUE.
#LDAUTYPE = 2
#LDAUL=2 -1
#LDAUU = 7.5 0
#LDAUJ = 0 0
#LDAUPRINT = 0

Ionic Relaxation
IBRION    =    3        #   ionic relaxation:0-MD 1-quasi-Newton 2-CG   DOS/EK =-1
ISIF      =    2
NSW       =    200    #   DOS/EK=0  relax = 400
EDIFFG    =   2E-05  
POTIM     =    0.01

LVDW =.TRUE.

NEB related tags
SPRING = -5
IMAGES = 6
ICHAIN = 0
LCLIMB = T

NPAR = 3


DOS related values
ISMEAR    =    0
SIGMA     =    0.05
#NEDOS     =    1200   # for more dos

ALgorithm
ALGO      =    FAST

#spin-orbital interaction
#LSORBIT   = .TRUE.
#SAXIS     = 0 0 1
#MAGMOM    =
#GGA_COMPAT = .FALSE.
#LMAXMIX    =  4
#ADDGRID    = .TRUE.
回复此楼

» 猜你喜欢

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

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

hakuna

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
wd12309812: 金币+5, ★★★很有帮助 2015-04-21 10:13:53

晕,真服了你了,你指定就跑200步还抱怨200终止了
2楼2015-03-19 18:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
wd12309812: 金币+5, ★★★很有帮助 2015-04-21 10:14:01
正如楼上所说。。。。NELM=200,还有你反应物和产物在结构优化的时候电子能200步以内能收敛不?
3楼2015-04-21 09:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

寂寞烟雨

铜虫 (小有名气)


一般电子步超过30步左右还不能收敛就很难收敛啦
4楼2015-04-21 10:03:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wd12309812

铜虫 (正式写手)

引用回帖:
4楼: Originally posted by 寂寞烟雨 at 2015-04-21 10:03:12
一般电子步超过30步左右还不能收敛就很难收敛啦

就是说  NELM不要太大,计算时候最好不要超过30是吗?
5楼2015-04-21 10:14:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wd12309812

铜虫 (正式写手)

问题已解决,可能跟POTIM有关系,
通过调EDIFF EDIFFG逐步收敛


想问一下,电子步一般设置多少呢?
6楼2015-04-21 10:16:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

寂寞烟雨

铜虫 (小有名气)


引用回帖:
5楼: Originally posted by wd12309812 at 2015-04-21 10:14:47
就是说  NELM不要太大,计算时候最好不要超过30是吗?...

一般采用默认值60就可以    如果超过了40步还不收敛   你看看是不是结构写的太差或者可以对IALGO、LDIAG和混合参数进行一些合理的设置    具体的你可以看手册
7楼2015-04-21 10:58:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

想请教一下楼主
我看你有开 LVDW
你开VDW时 是挑哪种VDW的?
你是用哪种POTCAR 去配合该种VDW?
8楼2015-04-21 11:35:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wd12309812 的主题更新
信息提示
请填处理意见