计算化学公社

标题: 用dimer找出的过渡态怎么确定是正确连接产物和反应物？ [打印本页]

作者
Author: ikea1984 时间: 2023-1-28 14:13
标题: 用dimer找出的过渡态怎么确定是正确连接产物和反应物？
一个反应，用CI-NEB死活不收敛，总是算一阵子就跳动，然后能量变高，力变大。用neb2dim搞了个dimer计算，找出来的用频率计算验证只有一个虚频。但是这种情况怎么确定这个过渡态就是给定的产物和反应物的？非周期性一般用IRC来确定，周期性体系真不知道了。。。看虚频的振动方向？感觉也只适合特别简单的1-2个原子移动的情况

作者
Author: ahxb 时间: 2023-1-30 20:54
目前的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计算中仅放开在数值频率计算中放开的原子，或许是一种可行的做法。

作者
Author: yuexiatx 时间: 2023-2-10 16:35
一般单虚频对应反应方向就行了体系很复杂的话我也不知道怎么确定了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)