IRC分析后，反应产物结构优化不收敛

OvoTin · 发表于 Post on 2025-3-10 21:10:11

各位老师好，我需要计算的是羟基自由基的RAF反应能垒。在对IRC分析后的反应物进行结构优化时没有出现明显震荡，但一直不收敛，也没有报错。请问需要加收敛条件重新计算吗？
前期过渡态的结构优化只有一个虚频且震动方向符合预期。
软件版本号为 Gaussian 16, Revision A.03

↑这是反应物的能量曲线

↑这是IRC的结果

输入及输出文件.rar (1.41 MB, 下载次数 Times of downloads: 5)
↑任务的输入及输出文件

Solitude198 · 发表于 Post on 2025-3-10 22:01:00

你放大看看确实有点震荡，但是整体能量是下降的，m062x是这样的，继续跑还是不行按照sob的指示来或者试试其他泛函

Stardust0831 · 发表于 Post on 2025-3-11 00:20:32

M06-2X对格点要求极高，我个人认为一个有效的方法是int=superfine来使用更好的积分格点。我实测发现g16默认的int=ultrafine对M06-2X经常是不够用的。

sobereva · 发表于 Post on 2025-3-11 01:10:07

老生常谈的问题
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

M06-2X结合SMD是最经典的几何优化难收敛的组合（上文和下文都说了），完全是自找难收敛
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

用B3LYP-D3(BJ)或wB97XD结合IEFPCM干这个是好得多的选择，耗时更低、结果不会更差，且出现震荡难收敛的概率连当前1/10都没有。

OvoTin · 发表于 Post on 2025-3-11 12:10:34

Solitude198 发表于 2025-3-10 22:01
你放大看看确实有点震荡，但是整体能量是下降的，m062x是这样的，继续跑还是不行按照sob的指示来或者试试 ...

感谢

OvoTin · 发表于 Post on 2025-3-11 12:11:08

Stardust0831 发表于 2025-3-11 00:20
M06-2X对格点要求极高，我个人认为一个有效的方法是int=superfine来使用更好的积分格点。我实测发现g16默认 ...

多谢帮助，我会试试看的

OvoTin · 发表于 Post on 2025-3-11 12:12:01

sobereva 发表于 2025-3-11 01:10
老生常谈的问题
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

感谢老师回复，我再认真学习一下您的博文

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] IRC分析后，反应产物结构优化不收敛

浏览过的版块