计算化学公社

标题: 过渡态优化中，如何让它保留已有虚频继续优化 [打印本页]

作者
Author: Dave666 时间: 2025-8-23 14:35
标题: 过渡态优化中，如何让它保留已有虚频继续优化
本帖最后由 Dave666 于 2025-8-23 14:51 编辑

在找一个过渡态的过程中，在监测过程中，发现其中一个优化比较接近收敛的点，四个指标如下：
Maximum Force             0.000025    0.000450    YES
RMS    Force                0.000003    0.000300    YES
Maximum Displacement 0.038291    0.001800    NO
RMS    Displacement    0.009457    0.001200    NO

对它进行振动分析，有一个虚频，振动频率数值只有-13（同类反应虚频在-70左右），但符合预期振动方向。
以这个点为起始，继续优化，期望找到它的过渡态，但最后优化得出的结果，虚频振动已经和相关的原子完全无关了，而是配体上的甲基的旋转。

请问老师们，有没有办法保留这个虚频，优化得到相应的过渡态？

作者
Author: wzkchem5 时间: 2025-8-23 14:42
可以先跑个IRC看看。有的反应，过渡态虚频振动方向与预期不符，但IRC能走到正确的反应物和产物，此时以IRC结果为准

作者
Author: Dave666 时间: 2025-8-23 14:50
本帖最后由 Dave666 于 2025-8-23 15:02 编辑

wzkchem5 发表于 2025-8-23 14:42
可以先跑个IRC看看。有的反应，过渡态虚频振动方向与预期不符，但IRC能走到正确的反应物和产物，此时以IRC ...

感谢老师的建议。
刚我的问题没叙述完整：最后优化得到的结果，虚频振动方向已经完全和核心反应原子无关了，而是配体上甲基的旋转。
irc跑的话应该是不太可能走到预期反应物和产物了。

作者
Author: 北大-陶豫 时间: 2025-8-23 14:59
冻结反应中心的键长键角等，进行 opt=ModRed freq 计算，看看振动模式是否还正常，正常的话再继续做 opt=(TS,ReadFC,NoFreeze,NoEigen) 试试

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