24小时热门版块排行榜    

查看: 1019  |  回复: 1

liufeila

新虫 (小有名气)

[求助] CINEB过渡态不收敛

NEW=2000,现在算了5天了,算到1350多步了,还没有收敛。做的是甲醇分子在锐钛矿晶面上脱去一个H原子。始态和终态都是优化好的结构。先把INCAR贴出来,请大侠们帮忙看看。万分感谢

Start parameter for this run
LPLANE = .TRUE.
ISTART = 0
ICHARG= 2

NEB tag
LCLIMB = .TRUE.
ICHAIN = 0
IMAGES = 4
SPRING = -5
EFIRST = -679.747953 ev
ELAST = -679.966761 ev
MAXMOVE = 0.2
TIMESTEP = 0.1
NSIM = 4

Electronic Relaxation
ENCUT = 400 ev
PREC = Normal
ISIF = 2
EDIFF = 1E-4
EDIFFG = -0.01
NELMIN= 6
LREAL= Auto
ALGO = Fast

DOS related values
ISMEAR = 0
SIGMA = 0.1
NEDOS = 300
ISPIN = 2

Ionic Relaxtion
ISMEAR = 0
NSW = 2000
IBRION = 3
POTIM = 0.0
IOPT = 7
NBLOCK = 1

Write flag
LWAVE = .FALSE.
LCHARG = .FALSE.
LVTOT = .FALSE.
LORBIT = .FALSE.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。 2012-06-25 17:36:50
liufeila: 金币+2, 有帮助, 好的, 我先试试 2012-06-25 18:07:46
楼主可以试试IOPT=2,用CG方法优化给出势能面的曲率应该会好些,IOPT=7是基于quick-min的,这个方法对于不精确的力还好些,就是说粗略找到过渡态还是可以的,但是对于楼主这个-0.01的力这个方法就不那么适用了。For high forces (far from the minimum) or inaccurate forces (close to the minimum) the quick-min or FIRE methods are recommended. These two methods do not rely on curvatures, and tend to be less aggressive, better behaved, but also less efficient than CG/LBFGS.
这是网站上那个说明,希望对你有帮助。
2楼2012-06-25 16:33:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liufeila 的主题更新
信息提示
请填处理意见