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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于CI-NEB
221/3返回列表123下一页
查看: 5053  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
关闭
请填处理意见
关闭