24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1797  |  回复: 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 的主题更新
信息提示
请填处理意见