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

24小时热门版块排行榜    

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

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bingmou(金币+1):感谢分享 2010-12-08 19:05:09
引用回帖:
Originally posted by lgxyz at 2010-12-08 14:50:44:


弹簧力是沿着两个image连线的切线方向，是为了保证每个image之间的联系；而这里讲到的力是构型在势能面上，是为了让这个构型优化为局域极小；这二者同时作用才能找到一级鞍点。我是这样理解的，不对之处，请多 ...

按道理说应该是这样的，但VASP自带的EB计算如果取力的收敛标准不知道它指的力是什么力，能量标准时几个镜像点能量之和，如果受力也是几个镜像点之和弹簧力应该能够消掉。但我算过几次无一例外都是在采用力的收敛标准时结果出现了严重错误。用CLNEB是可以的
具体为什么会这样还没搞懂
赞一下(2人) 回复此楼
心之所向无惧无悔
11楼2010-12-08 15:44:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
嗯，我说的就是CI-NEB
赞一下 回复此楼
12楼2010-12-08 19:05:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wqg1806

木虫 (著名写手)

科研农民工
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-12-08 21:22:55
我以前也遇到过，我觉得是路径设计的不合理。中间可能还要经过一些过渡态
赞一下(2人) 回复此楼
一路走到黑
13楼2010-12-08 21:07:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by wqg1806 at 2010-12-08 21:07:27:
我以前也遇到过，我觉得是路径设计的不合理。中间可能还要经过一些过渡态

多谢！我隐约也觉得是这样的！
赞一下 回复此楼
14楼2010-12-08 21:20:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

04nylxb

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by hakuna at 2010-12-08 13:39:49:


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

请教，如何消掉多出来的虚频？非常感谢
赞一下(1人) 回复此楼
集中精力发文章
15楼2010-12-09 08:57:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢指点 2010-12-09 11:21:46
怎么调我就不知道了，但NEB确实是以image受的力为收敛标准的。NEB的力是包含两部分：切线方向的弹性和垂直于切线方向的真实力。
赞一下(2人) 回复此楼
16楼2010-12-09 11:13:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-12-10 00:10:27
引用回帖:
Originally posted by sunyang1988 at 2010-12-08 14:00:45:

请教，如何消掉多出来的虚频？非常感谢

这没什么标准方法，经验，直接调坐标、部分优化......

[ Last edited by hakuna on 2010-12-10 at 08:16 ]
赞一下(2人) 回复此楼
17楼2010-12-09 23:47:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿克飞顿

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
引用回帖:
7楼: Originally posted by hakuna at 2010-12-08 13:39:49
不对吧，可以用能量标准的。CI-NEB收敛后一般情况下能量高的image就是过渡态了，个别情况会是二阶鞍点，消掉多出来的虚频就可以了，不理解你为什么那样做.......

您好，请问您一个问题。
我用CI-NEB找过渡态，分成三步计算的，第一步EDIFFG设为-0.5，经历二三十步达到力的收敛判据后，第二部设为-0.1.一般要经历100多步才能达到力收敛，第三部设为-0.05一般要两三百才能收敛。期间我用nebconverge.pl查看力的收敛趋势，发现中间多次出现跳跃，即本来已经收敛到-0.07了，结果突然跳跃到-0.3然后又慢慢收敛下来。
请问我这样得到的过渡态可信吗？？？
另外加入我插入了四个构型，1 2 4都达到力判据了，3还差-0.01左右就能达到连判据，而2才是想要的过渡态，这是我能不能就停止计算，用2作为我的过渡态或者把2的构型再拿出去结构优化完作为过渡态？？
赞一下 回复此楼
18楼2015-09-10 10:54:01
 已阅   回复此楼   关注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的回帖

阿克飞顿

金虫 (正式写手)
★
小木虫: 金币+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客户端
赞一下 回复此楼
20楼2015-09-10 11:49:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgxyz 的主题更新
222/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭