| 查看: 1030 | 回复: 1 | |||
[求助]
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. |
回复此楼
likai0106
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 38 (小学生)
- 金币: 1352.4
- 帖子: 233
- 在线: 68.5小时
- 虫号: 634278
- 注册: 2008-10-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流。 2012-06-25 17:36:50
liufeila: 金币+2, ★有帮助, 好的， 我先试试 2012-06-25 18:07:46
感谢参与，应助指数 +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
