能否实现每次优化都做频率分析

mfdsrax2

Gaussian或者其他量化软件能否实现在calcall的基础上，每优化一次都自动做频率分析，方便判断优化方向

Uus/pMeC6H4-/キ

如果是寻找过渡态时要追踪优化方向，可以在输出文件里Initial Parameters列出的内坐标中找到自己关心的键长/键角/二面角项，然后观察每一步Eigenvectors required to have negative eigenvalues列出的本征值，看关心的内坐标项是否存在、贡献是否为负且绝对值较大（参考http://bbs.keinsci.com/thread-14081-1-1.html的介绍）。这并不需要什么每步计算Hessian矩阵（即calcall关键词的作用）甚至每步做振动分析之类的麻烦事，只要有#p就自然会输出了（我不清楚不开#p的情况如何，不过也没理由不用#p）。

编辑：如果是寻找极小点，追踪优化方向似乎意义不大，毕竟单纯的几何优化本就不足以找全所有极小点。真要看从一个初始结构出发怎么走到不同的极小点附近的话，或许可以尝试http://bbs.keinsci.com/thread-19372-1-1.html说的downhill算法（注意14和15楼的问答），不过downhill还是不能代替几何优化。

wal

优化方向是啥意思 是想看每一步的虚频分别长啥样吗

试了下 #p也不输出 可能要自己写脚本算每一步的freq吧

Uus/pMeC6H4-/キ

wal 发表于 2025-5-6 12:46
优化方向是啥意思 是想看每一步的虚频分别长啥样吗

试了下 #p也不输出 可能要自己写脚本算每一步的freq ...

Eigenvectors required to have negative eigenvalues应该是opt(ts)任务寻找过渡态才在L103输出的，opt任务寻找极小点确实没有。也可能和优化所用坐标有关，我上面的结论是对用冗余内坐标的默认情形说的。

说实话，我觉得优化路径上各点做振动分析这种操作比较可疑，毕竟这些点可能离驻点很远，附近的势能面形状可能比较奇怪，此时谐振近似并不足以很好描述局域性质。这个问题论坛上也有讨论过，我记得是牵涉无能垒反应的自由能变和速率常数计算之类的语境，不过暂时找不到帖子（）