计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1336|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 用dimer找出的过渡态怎么确定是正确连接产物和反应物?

[复制链接 Copy URL]

102

帖子

0

威望

2672

eV
积分
2774

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
一个反应,用CI-NEB死活不收敛,总是算一阵子就跳动,然后能量变高,力变大。用neb2dim搞了个dimer计算,找出来的用频率计算验证只有一个虚频。但是这种情况怎么确定这个过渡态就是给定的产物和反应物的?非周期性一般用IRC来确定,周期性体系真不知道了。。。看虚频的振动方向?感觉也只适合特别简单的1-2个原子移动的情况

131

帖子

0

威望

2523

eV
积分
2654

Level 5 (御坂)

2#
发表于 Post on 2023-1-30 20:54:03 | 只看该作者 Only view this author
目前的VASP wiki上说用IBRION=40可以进行IRC计算,见https://www.vasp.at/wiki/index.php/IRC_calculations,具体从哪个版本开始支持我在官网上没查到。但是有三个问题:
(1) 提示应该用最低EDIFFG=-0.005的收敛限,比通常用的-0.02甚至-0.03严不少;
(2) 官网上没给出任何例子;
(3) 需要提供想验证的虚频的振动模式,但使用有限差分法计算数值频率时,常常仅放开部分原子,甚至为了节约计算资源,一般放开的原子比找过渡态时放开的原子还少,这样通过有限差分计算得到的虚频振动模式必然和将原子固定设置与找过渡态时保持一致有一定差别,对IRC路径可能有影响。当然,仅用于验证过渡态是否计算正确的话,在IRC计算中仅放开在数值频率计算中放开的原子,或许是一种可行的做法。

31

帖子

0

威望

264

eV
积分
295

Level 3 能力者

3#
发表于 Post on 2023-2-10 16:35:01 | 只看该作者 Only view this author
一般单虚频对应反应方向就行了 体系很复杂的话我也不知道怎么确定了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-25 19:16 , Processed in 0.186284 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list