24小时热门版块排行榜    

查看: 1961  |  回复: 15

wup123

铁虫 (小有名气)

[求助] CI-neb 收敛问题求教

最近算几个NEB,搞的好头疼,不收敛。。。。
我取了4个image, 每个image中FORCES: max atom, RMS都在0.06左右,而且上下震荡,就是不收敛;一般我取POTIM=0.2,有的作业POTIM都调到了0.1,力还在震荡,我的收敛标准是EDIFFG=-0.03.
  请问, 还能去调些啥参数促进收敛呢?先谢谢大家了!!!
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

liuziran

金虫 (小有名气)

【答案】应助回帖


贺仪(金币+1): 多谢指教! 2011-12-02 21:45:42
uuv2010: 回帖置顶 2011-12-02 22:35:53
wup123(金币+5): 2011-12-04 22:50:34
1. 一般要做两次优化,初次收敛标准可以粗糙点(-0.5),第二次再苛刻些(-0.1);
2. 初末两个状态一定要搞准,然后产生的images也要一个一个的检查,不要相信脚本产生的东西。
3楼2011-12-02 21:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wup123

铁虫 (小有名气)

送鲜花一朵
谢谢机器人了,但是不大相关。。。。。
2楼2011-12-02 17:45:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

【答案】应助回帖


uuv2010(金币+1): 欢迎参加讨论 2011-12-03 10:18:02
wup123(金币+5): 2011-12-04 22:50:42
引用回帖:
1楼: Originally posted by wup123 at 2011-12-02 17:43:59:
最近算几个NEB,搞的好头疼,不收敛。。。。
我取了4个image, 每个image中FORCES: max atom, RMS都在0.06左右,而且上下震荡,就是不收敛;一般我取POTIM=0.2,有的作业POTIM都调到了0.1,力还在震荡,我的收敛 ...

震动的厉害可能可能是因为image取少了。


你POTIM = 0.2这个用的是vasp自带的优化器,改成vtstool的优化器行不?IBRION =3
POTIM = 0
OPTIM = 似乎是这样的。这个收敛比较快。只要鞍点收敛了就可以停下来。然后分两段分别做优化路径。

你的两个端点之间的距离多少?
4楼2011-12-02 23:04:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)

谢谢你了,我在仔细查查看。
5楼2011-12-04 22:49:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by liuziran at 2011-12-02 21:09:25:
1. 一般要做两次优化,初次收敛标准可以粗糙点(-0.5),第二次再苛刻些(-0.1);
2. 初末两个状态一定要搞准,然后产生的images也要一个一个的检查,不要相信脚本产生的东西。

你是说的IOPT这个参数吧,我是这么理解的,要是你刚开始做neb的时候,这时候力很大,可以用IOPT=2,这样是CG的方法,可以快速的收敛;当力收敛的很小的时候,再采用IOPT=3还是4呢?对于3和4的区别不大明白,这种方法我没试过,不过我可以尝试下。初始态间的距离有影响吗?我初始态间的距离是2.63.
6楼2011-12-04 23:03:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by cenwanglai at 2011-12-02 23:04:39:
震动的厉害可能可能是因为image取少了。


你POTIM = 0.2这个用的是vasp自带的优化器,改成vtstool的优化器行不?IBRION =3
POTIM = 0
OPTIM = 似乎是这样的。这个收敛比较快。只要鞍点收敛了就可以停下 ...

你是说的IOPT这个参数吧,我是这么理解的,要是你刚开始做neb的时候,这时候力很大,可以用IOPT=2,这样是CG的方法,可以快速的收敛;当力收敛的很小的时候,再采用IOPT=3还是4呢?对于3和4的区别不大明白,这种方法我没试过,不过我可以尝试下。初始态间的距离有影响吗?我初始态间的距离是2.63.
7楼2011-12-04 23:04:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by liuziran at 2011-12-02 21:09:25:
1. 一般要做两次优化,初次收敛标准可以粗糙点(-0.5),第二次再苛刻些(-0.1);
2. 初末两个状态一定要搞准,然后产生的images也要一个一个的检查,不要相信脚本产生的东西。

刚才我回复错了,本来回复下面这个人的。恩,我检查下。
8楼2011-12-04 23:05:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-12-05 19:27:39
引用回帖:
6楼: Originally posted by wup123 at 2011-12-04 23:03:21:
你是说的IOPT这个参数吧,我是这么理解的,要是你刚开始做neb的时候,这时候力很大,可以用IOPT=2,这样是CG的方法,可以快速的收敛;当力收敛的很小的时候,再采用IOPT=3还是4呢?对于3和4的区别不大明白,这种 ...

这个距离比较短,算起来应该不是很困难吧~~
9楼2011-12-04 23:12:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

引用回帖:
3楼: Originally posted by liuziran at 2011-12-02 21:09:25:
1. 一般要做两次优化,初次收敛标准可以粗糙点(-0.5),第二次再苛刻些(-0.1);
2. 初末两个状态一定要搞准,然后产生的images也要一个一个的检查,不要相信脚本产生的东西。

我想问一下您说的那个先粗优化,再细优化,能量的收敛标准等用不用变呢?
我之前做过粗优化优化到0.5,0.1这都收敛了,细优化想做到0.01,可是到了0.05后就怎么都不收敛了,而且力越来越大了。最后尝试做0.04还是不收敛,请问这个情况怎么处理啊?谢谢
10楼2011-12-05 09:27:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wup123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 290分有科研经历,论文在投 +6 腻腻gk 2026-03-14 6/300 2026-03-16 09:39 by 闲人终南山
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见