计算化学公社

标题: 过渡态有多个虚频，尝试了很多办法都解决不了 [打印本页]

作者
Author: 我是大帅逼 时间: 2020-10-15 14:34
标题: 过渡态有多个虚频，尝试了很多办法都解决不了
{VASP!!!!!!!!!!!!!!!!!表面催化！！！！！！}要疯了，为啥我无论怎么优化,最后频率验证总是有多个虚频啊，而且是没法忽略的那种。试了以下方法：
1，IBRION=1,有的优化了无数步，力很小了。有的收敛了，一算频率。。。2个

然后提高收敛精度，还是不行
2，dimer，最后的确优化到只剩一个虚频，但是结构变了，并不是我想要的结构了（大家用dimer的时候没有遇到这种情况吗！！！最后的结构显然不合理）

我要放弃了，是不是有多个虚频就不是过渡态啊，频率这种东西也太难搞了，我直接吐了

另外，为解决这个问题，主要参考了刘大佬的贴子：http://bbs.keinsci.com/thread-13065-1-1.html
针对其提出的过渡态优化不收敛依次尝试以下方法做了以下改进
1，避免不合理的初始结构：我虽然用的nebmake，但所有结构我都检查并将不合理结构手动调整了。
2，CI-NEB插点的数目不合理，emmm，我觉得还行，因为我的距离都差不多，基本每个基元步骤都插4个点，一次CINEB就差不多了。
3，先用NEB或者CINEB进行粗收敛，再用DIMER/CINEB精收敛,NEB+DIMER/CINEB+DIMER都尝试了，最大的问题是DIMER之后结构变得不合理，举个例子CH3的分解，原本CINEB结束H原子（最终落在红圈的地方）在催化剂原子的顶位(侧顶位)，dimer完，好家伙，直接给我优化到穴位了，这沃妮马

最后附上这狗日的频率，我本来已经很宽容了，好多基元反应100以下的虚频我都忽略了，可是还是有很多小瘪犊子不行

作者
Author: qmlearner 时间: 2020-10-15 14:47
正在跑dimer，第三步curvature还是正值的我正在瑟瑟发抖(T_T)

作者
Author: lyzhang 时间: 2023-3-13 16:56
请问没有虚频怎么办呢？提高了电子和力精度也不行，但是能垒数值和路径是和文献一致的

作者
Author: sjlgaga 时间: 2023-6-27 15:56
请问楼主是怎么解决的

作者
Author: xmuqyj 时间: 2023-8-2 21:20

qmlearner 发表于 2020-10-15 14:47
正在跑dimer，第三步curvature还是正值的我正在瑟瑟发抖(T_T)

很正常

作者
Author: xmuqyj 时间: 2023-8-2 21:21
忽略100以下的虚频，这样会导致零点能不准确

作者
Author: zhang1999501 时间: 2023-8-2 21:35
可以尝试一下qvasp程序，里面有一个命令qvasp -zc 可以消除过渡态虚频。

作者
Author: xmuqyj 时间: 2023-8-5 10:39

zhang1999501 发表于 2023-8-2 21:35
可以尝试一下qvasp程序，里面有一个命令qvasp -zc 可以消除过渡态虚频。

结构优化虚频可以消除吗

作者
Author: zhang1999501 时间: 2023-8-5 17:43

xmuqyj 发表于 2023-8-5 10:39
结构优化虚频可以消除吗

手册中写的是消除过渡态的小虚频，确实有效，但是否可以消除结构优化的虚频，你可以咨询qvasp作者易文才。

作者
Author: 火星哥 时间: 2023-8-6 13:56

xmuqyj 发表于 2023-8-2 21:21
忽略100以下的虚频，这样会导致零点能不准确

您好，如果优化出【同种反应途径，成键角度大的】100以下的唯一虚频，当提升泛函基组进行溶剂化时发现零点能不准了，但是我们也定位到【同种反应途径不同成键角度的过渡态，角度小的其虚频在200-500之间】，同类反应途径随泛函基组的改变没有影响随角度变化变化的趋势。但不同的反应性途径，随着泛函基组的提升，反应的优势路径发生了改变。这种情况需要将我的【大成键角度】过渡态虚频提高才可以吗？

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