24小时热门版块排行榜    

查看: 2280  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liupeitao20

金虫 (小有名气)

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

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

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
信息提示
请填处理意见