24小时热门版块排行榜    

查看: 1916  |  回复: 12

egogg

木虫 (小有名气)

[求助] NEB 结果分析

过渡态方面我是一个新手,求教一下NEB的结果如何分析:
CI-NEB的脚本给出了如下的反应路径,该如何分析?

反应过程的动画跟我想象的差不多,还需要做些什么工作?
非常感谢!

[ Last edited by egogg on 2012-5-15 at 15:15 ]
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egogg

木虫 (小有名气)

自己顶一个。
2楼2012-05-15 15:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongb123456

铜虫 (小有名气)

帮你顶一个 我也想知道!!
LZ知道了 告诉下我呗 谢谢
3楼2012-05-15 17:04:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egogg

木虫 (小有名气)

来个有经验的人啊,自己再顶顶。
4楼2012-05-16 12:58:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢指教 2012-05-20 08:41:07
egogg: 金币+10, ★★★★★最佳答案 2012-05-20 11:49:29
这个告诉你03下的结构可能是TS,可以对这个结构做一下频率计算,已确认是否真的是TS,如果是,你可以求出能垒(00与03能量差)以及这个过程的能量变化(00与07能量差)。另外,TS结构更像00,说是前过渡态的话,有点勉强。
如果你做的图是正确的话!
其它似乎没什么要说的了
5楼2012-05-19 21:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egogg

木虫 (小有名气)

引用回帖:
5楼: Originally posted by hakuna at 2012-05-19 21:28:16:
这个告诉你03下的结构可能是TS,可以对这个结构做一下频率计算,已确认是否真的是TS,如果是,你可以求出能垒(00与03能量差)以及这个过程的能量变化(00与07能量差)。另外,TS结构更像00,说是前过渡态的话,有 ...

明白了,非常感谢。

也就是说需要进行一次频率验证就行了。

但是能不能指点一下,为什么说03的结构是过渡态,这个图显示的不是04为最高峰么?还有就是如果00有些像前过渡态的话,我能不能从01或02开始重新作NEB计算?因为就像您说的那样,从反应的动画来看,00和后面的差很远(相关性很差,我取初始态的时候取得离我想研究的反应路径远有些远了)。

还有就是做频率验证用vasp自带的算算Gamma点就行,还是用vstool的官方网站上介绍的Dynamical Matrix方法比较好一点,非常感谢!
6楼2012-05-20 11:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)

★ ★ ★
贺仪: 金币+3, 谢谢指教! 2012-05-21 08:17:12
引用回帖:
6楼: Originally posted by egogg at 2012-05-20 11:48:50:
明白了,非常感谢。

也就是说需要进行一次频率验证就行了。

但是能不能指点一下,为什么说03的结构是过渡态,这个图显示的不是04为最高峰么?还有就是如果00有些像前过渡态的话,我能不能从01或02开始重新 ...

1.但是能不能指点一下,为什么说03的结构是过渡态,这个图显示的不是04为最高峰么?
我以为你的第一个点是00,如果不是,就按你的次序后推,总之,能量最高者可能是TS
2.还有就是如果00有些像前过渡态的话,我能不能从01或02开始重新作NEB计算?
00里是反应物,不会像过渡态吧。找TS是在反应物和产物间搜索的,一般用CNEB在中间插若干点,找TS命中率还是蛮高的。也有人采用逐步逼近的办法,譬如每次插入3个点,待NEB收敛后,根据情况,将两个端点向内收缩,减小TS所在的区间,反复多次,当有TS区间收缩到足够小时,认为得到了TS,不过,这样做显得比较麻烦,用的人很少。
3.将IBRION=5直接做频率计算,具体细节在参考一下手册吧
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
7楼2012-05-20 18:30:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egogg

木虫 (小有名气)

引用回帖:
7楼: Originally posted by cavediger at 2012-05-20 18:30:08:
1.但是能不能指点一下,为什么说03的结构是过渡态,这个图显示的不是04为最高峰么?
我以为你的第一个点是00,如果不是,就按你的次序后推,总之,能量最高者可能是TS
2.还有就是如果00有些像前过渡态的话,我 ...

谢谢,明白了。
8楼2012-05-20 19:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenyan0517

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by hakuna at 2012-05-19 21:28:16:
这个告诉你03下的结构可能是TS,可以对这个结构做一下频率计算,已确认是否真的是TS,如果是,你可以求出能垒(00与03能量差)以及这个过程的能量变化(00与07能量差)。另外,TS结构更像00,说是前过渡态的话,有 ...

想请教一下,NEB计算时,电子步收敛时一直在震荡,一正一负交替,有什么好办法解决吗?
9楼2012-05-21 09:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)

★
贺仪: 金币+1, 谢谢指教! 2012-05-21 10:57:26
引用回帖:
9楼: Originally posted by wenyan0517 at 2012-05-21 09:50:39:
想请教一下,NEB计算时,电子步收敛时一直在震荡,一正一负交替,有什么好办法解决吗?

这个还真没遇到过,有两点可以考虑一下:
1.更换电子步的优化方法
2.调整一下电荷密度的混合方法
仅供参考

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

我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
10楼2012-05-21 10:38:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 egogg 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 07:55 by GBJlw5zzHWpZ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-06 07:55 by GBJlw5zzHWpZ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 05:54 by qL1cU46YIdK0
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 05:33 by qL1cU46YIdK0
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 05:15 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-06 04:58 by qL1cU46YIdK0
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-06 02:37 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:26 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:40 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见