| 查看: 1090 | 回复: 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. |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有20人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
