计算化学公社

标题: 能否实现每次优化都做频率分析 [打印本页]

作者
Author:
mfdsrax2    时间: 2025-5-6 11:16
标题: 能否实现每次优化都做频率分析
Gaussian或者其他量化软件能否实现在calcall的基础上,每优化一次都自动做频率分析,方便判断优化方向

作者
Author:
Uus/pMeC6H4-/キ    时间: 2025-5-6 11:46
本帖最后由 Uus/pMeC6H4-/キ 于 2025-5-6 12:00 编辑

如果是寻找过渡态时要追踪优化方向,可以在输出文件里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还是不能代替几何优化。

作者
Author:
wal    时间: 2025-5-6 12:46
本帖最后由 wal 于 2025-5-6 13:07 编辑

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

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

作者
Author:
Uus/pMeC6H4-/キ    时间: 2025-5-6 13:47
本帖最后由 Uus/pMeC6H4-/キ 于 2025-5-6 13:52 编辑
wal 发表于 2025-5-6 12:46
优化方向是啥意思 是想看每一步的虚频分别长啥样吗

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

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

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3