24小时热门版块排行榜    

查看: 2384  |  回复: 5

liupeitao20

金虫 (小有名气)

[求助] CL_NEB方法搜索过渡态 怎么准确获得saddle-point 的几何及相应的振动频率?

我们知道 用CL_NEB计算能垒的或者搜索过渡态的时候,先优化好初始态和终了态,然后做插值获得几个Images,然后ISIF=2 等相应设置后计算,收敛后得到能量,即获得能垒,可是想准确得到鞍点处过渡态的几何结构,如何获得?我是这么做的, 对于能量最高的Image的CONTACR拷贝至POSCAR然后以力作为收敛标准,进行优化,优化收敛后进行频率计算,获得该处的频率只有一个虚频,即找到了该过渡态?不知我做的是否正确,还请牛人指点?因为我得到过渡态的结构的频率虚频不止一个,所以我很怀疑自己的方法!还有就是大家怎么去算一个体系的振动频率!是直接用VASP自带的还是自己编程构造dynamical matrix然后求解频率?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流 2012-09-25 09:52:44
liupeitao20: 金币+5, ★★★很有帮助, 说的很有道理 2012-09-26 20:34:52
CNEB收敛后,择能量最高的image,不要优化,算频率,如果有且仅有一个虚频,其振动模式沿断键或成键方向可是为TS,若有多个虚频,把不需要的虚频消掉即可,消虚频方法各人可能有不同做法.....
至于怎么算频率,查下IBRION......
几何结构吗,不就在OUTCAR里吗?
2楼2012-09-25 06:41:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liupeitao20

金虫 (小有名气)

引用回帖:
2楼: Originally posted by hakuna at 2012-09-25 06:41:46
CNEB收敛后,择能量最高的image,不要优化,算频率,如果有且仅有一个虚频,其振动模式沿断键或成键方向可是为TS,若有多个虚频,把不需要的虚频消掉即可,消虚频方法各人可能有不同做法.....
至于怎么算频率,查下 ...

谢谢回复!能量最高点处如有多个虚频,只能说是为完全达到过渡态,为什么要去人为的消除虚频!频率计算VASP自带的本人熟悉,想问一下有没有更好的code可以自己构造动力学矩阵,然后求解频率!
3楼2012-09-25 20:04:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

SeanMa2010

铜虫 (初入文坛)

★
WDD880227: 金币+1, 谢谢参与 2012-09-26 16:36:25
能量最高的Image不能优化啊,本来就是在saddle point附近,是非常不稳定的,稍一优化直接优化到某个能量极小点了,
4楼2012-09-25 20:22:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

★
WDD880227: 金币+1, 多谢耐心指教 2012-09-26 16:36:34
引用回帖:
3楼: Originally posted by liupeitao20 at 2012-09-25 20:04:55
谢谢回复!能量最高点处如有多个虚频,只能说是为完全达到过渡态,为什么要去人为的消除虚频!频率计算VASP自带的本人熟悉,想问一下有没有更好的code可以自己构造动力学矩阵,然后求解频率!...

多个虚频意味着高阶鞍点,并不是过渡态,过渡态应该只有一个虚频且其振动模式与成键或断键模式匹配....
构造力矩阵的code没用过....
5楼2012-09-26 09:15:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)

引用回帖:
2楼: Originally posted by hakuna at 2012-09-25 06:41:46
CNEB收敛后,择能量最高的image,不要优化,算频率,如果有且仅有一个虚频,其振动模式沿断键或成键方向可是为TS,若有多个虚频,把不需要的虚频消掉即可,消虚频方法各人可能有不同做法.....
至于怎么算频率,查下 ...

借此地问下,vasp算的过渡态可以消虚频吗?具体怎么处理?
相逢是缘!
6楼2012-09-29 08:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liupeitao20 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见