应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于CI-NEB
221/3返回列表123下一页
查看: 4705  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lgxyz

铜虫 (正式写手)

[交流] 【求助】关于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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hakuna

木虫 (知名作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
20楼: Originally posted by 阿克飞顿 at 2015-09-10 11:49:21
那这种要经过两三百步才能达到力判据收敛的状况正常吗？为什么有的人说超过五六十就说明没必要继续算了呢？它们说的是电子步吧！
这是我的INCAR:
SYSTEM = TS1
ENCUT = 400
EDIFF  = 1E-4
EDIFFG = -0.05
GGA=P ...

1.NEB优化与一般结构优化还是有差别的，几百步甚至上千步才收敛都是正常的，这和你中间插入的image有关，插入的image比较准确，优化就快，否则就慢，遗憾的是现在还没有特别好的办法使插入的image比较准确，一些脚本或小程序都是采用线性插值方法，插出来的image是否合理，怎么适当调整，都是个人经验了，所以，同一个问题，不同人去做时，收敛结果还是可能存在差异的
2.电子步默认好像是60吧，定义太多电子步一般情况下意义不大
3.有建议初搜TS时用少量K点，待收敛后，加密K点接着精搜的做法
4.原则上讲，搜TS时，表面层要放开，但这样会极大地增加TS搜索成本；为节省时间，有将表面层固定搜TS的，完了以后在局部优化一下表面，扣除表面弛豫的影响；也有忽略表面影响的.....一般情况下，差别都很小
赞一下(2人) 回复此楼
21楼2015-09-11 10:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
18楼: Originally posted by 阿克飞顿 at 2015-09-10 10:54:01
您好，请问您一个问题。
我用CI-NEB找过渡态，分成三步计算的，第一步EDIFFG设为-0.5，经历二三十步达到力的收敛判据后，第二部设为-0.1.一般要经历100多步才能达到力收敛，第三部设为-0.05一般要两三百才能收敛。 ...

1.TS搜索过程中，力或能量出现跳跃应该都是正常现象，搜索TS本质上也是优化，而目前的优化都是与梯度相关的局域搜索方法，所以出现跳跃很正常.......
2.建议还是等NEB收敛后再决定
赞一下(1人) 回复此楼

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

阿克飞顿
19楼2015-09-10 11:24:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

bingmou

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
嗯，我遇到过，看着力越来越小，然后突然变大，然后越来越大，很郁闷
赞一下(1人) 回复此楼
2楼2010-12-07 22:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by bingmou at 2010-12-07 22:02:53:
嗯，我遇到过，看着力越来越小，然后突然变大，然后越来越大，很郁闷

谢谢 请问后来解决了吗？
赞一下 回复此楼
3楼2010-12-07 22:11:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-12-08 08:58:42
引用回帖:
Originally posted by bingmou at 2010-12-07 22:02:53:
嗯，我遇到过，看着力越来越小，然后突然变大，然后越来越大，很郁闷

你的收敛标准是不是力，改成能量试试
我之前算的收敛标准是力的都出现了这样的问题，结果出来很离谱
赞一下(2人) 回复此楼
心之所向无惧无悔
4楼2010-12-08 08:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by sunyang1988 at 2010-12-08 08:40:05:

你的收敛标准是不是力，改成能量试试
我之前算的收敛标准是力的都出现了这样的问题，结果出来很离谱

这个CI-NEB只能用力作为收敛标准。
赞一下 回复此楼
5楼2010-12-08 09:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢提示 2010-12-08 13:34:06
而且力作为收敛标准应该是比较合理的才对，能量收敛了有可能力还没有收敛，这样的话体系就存在很大应力
赞一下(2人) 回复此楼
6楼2010-12-08 13:32:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):鼓励讨论交流 2010-12-08 14:29:23
引用回帖:
Originally posted by lgxyz at 2010-12-08 09:20:28:


这个CI-NEB只能用力作为收敛标准。

不对吧，可以用能量标准的。CI-NEB收敛后一般情况下能量高的image就是过渡态了，个别情况会是二阶鞍点，消掉多出来的虚频就可以了，不理解你为什么那样做....
赞一下(2人) 回复此楼

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

阿克飞顿
7楼2010-12-08 13:39:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):鼓励讨论交流 2010-12-08 14:29:29
引用回帖:
Originally posted by bingmou at 2010-12-08 13:32:32:
而且力作为收敛标准应该是比较合理的才对，能量收敛了有可能力还没有收敛，这样的话体系就存在很大应力

普通的NEB是不能用力作为收敛标准的，因为存在的弹簧力，用力作为标准最后会得到非常荒谬的结果
赞一下(2人) 回复此楼
心之所向无惧无悔
8楼2010-12-08 14:00:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by hakuna at 2010-12-08 13:39:49:


不对吧，可以用能量标准的。CI-NEB收敛后一般情况下能量高的image就是过渡态了，个别情况会是二阶鞍点，消掉多出来的虚频就可以了，不理解你为什么那样做....

对呀 当每个image上的力小于收敛标准了，对应的能量最高的那个image就是过渡态。因为CI-NEB是基于力的算法。
赞一下 回复此楼
9楼2010-12-08 14:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
★
bingmou(金币+1):感谢分享 2010-12-08 19:04:47
引用回帖:
Originally posted by sunyang1988 at 2010-12-08 14:00:45:

普通的NEB是不能用力作为收敛标准的，因为存在的弹簧力，用力作为标准最后会得到非常荒谬的结果

弹簧力是沿着两个image连线的切线方向，是为了保证每个image之间的联系；而这里讲到的力是构型在势能面上，是为了让这个构型优化为局域极小；这二者同时作用才能找到一级鞍点。我是这样理解的，不对之处，请多多指点！
赞一下(1人) 回复此楼
10楼2010-12-08 14:50:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgxyz 的主题更新
221/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +5 Ncdx123456 2026-03-13 6/300 2026-03-15 23:39 by lovewei0727
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 求调剂，药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 083000环境科学与工程调剂 +8 mingmingry 2026-03-09 9/450 2026-03-11 10:23 by 沙漠之狐994
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭