24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lgxyz

铜虫 (正式写手)

应助: 6 (幼儿园)
金币: 54.6
散金: 66
红花: 8
帖子: 528
在线: 240.5小时
虫号: 1043996
注册: 2010-06-19
性别: GG
专业: 理论和计算化学

[交流] 【求助】关于CI-NEB 已有8人参与

各位，最近尝试着用CI-NEB找过渡态，先用IBRION=3 POTIM=0.1使各个images上的力小于1eV/ang，然后用IOPT=1优化，说是IOPT=1是更好的优化方法，在力小于一定程度的时候。但是结果是各个image上的力不但没有收敛的趋势，反而越来越大，不知大家有没有遇到这种情况？

回复此楼

» 本帖已获得的红花（最新10朵）

阿克飞顿

» 猜你喜欢

1楼 2010-12-07 16:54:14

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

阿克飞顿

金虫 (正式写手)

应助: 8 (幼儿园)
金币: 1426.6
散金: 37
红花: 2
帖子: 572
在线: 169.6小时
虫号: 2469736
注册: 2013-05-18
性别: GG
专业: 光学和光电子材料

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

送红花一朵

引用回帖:

19楼: Originally posted by hakuna at 2015-09-10 11:24:44
1.TS搜索过程中，力或能量出现跳跃应该都是正常现象，搜索TS本质上也是优化，而目前的优化都是与梯度相关的局域搜索方法，所以出现跳跃很正常.......
2.建议还是等NEB收敛后再决定...

那这种要经过两三百步才能达到力判据收敛的状况正常吗？为什么有的人说超过五六十就说明没必要继续算了呢？它们说的是电子步吧！
这是我的INCAR:
SYSTEM = TS1
ENCUT = 400
EDIFF = 1E-4
EDIFFG = -0.05
GGA=PE
ISTART = 0 ; ICHARG = 2
ISMEAR = 0;SIGMA=0.2
ISIF=2
LREAL = AUTO
NSW=1000
IBRION=3
POTIM=0.0
NELM = 500
IALGO=48
LDIAG=.TRUE.
ALGO = Fast
LWAVE = .TRUE.
LCHARG = .TRUE.
NPAR = 5;NSIM = 4
LPLANE = .TRUE.
LSCALU = .TRUE.
IOPT=3
ICHAIN=0
MAXMOVE=0.2
SPRING=-5
IMAGES=3
LCLIMB = .TRUE.

KPOINT
我用的2×2×1
体系有80多个原子共七层，下面四层固定

发自小木虫Android客户端