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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于CI-NEB
221/3返回列表123下一页
查看: 5109  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 我不理解！ +9 Edward_pc 2026-08-26 14/700 2026-08-26 18:03 by hsding2026
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭