计算化学公社

标题: IRC分析后，反应产物结构优化不收敛 [打印本页]

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

(, 下载次数 Times of downloads: 20)
↑这是反应物的能量曲线
(, 下载次数 Times of downloads: 19)
↑这是IRC的结果

(, 下载次数 Times of downloads: 8)
↑任务的输入及输出文件

作者
Author: Solitude198 时间: 2025-3-10 22:01
你放大看看确实有点震荡，但是整体能量是下降的，m062x是这样的，继续跑还是不行按照sob的指示来或者试试其他泛函
(, 下载次数 Times of downloads: 22)
(, 下载次数 Times of downloads: 19)

作者
Author: Stardust0831 时间: 2025-3-11 00:20
M06-2X对格点要求极高，我个人认为一个有效的方法是int=superfine来使用更好的积分格点。我实测发现g16默认的int=ultrafine对M06-2X经常是不够用的。

作者
Author: sobereva 时间: 2025-3-11 01:10
老生常谈的问题
量子化学计算中帮助几何优化收敛的常用方法
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都没有。

作者
Author: OvoTin 时间: 2025-3-11 12:10

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

感谢

作者
Author: OvoTin 时间: 2025-3-11 12:11

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

多谢帮助，我会试试看的

作者
Author: OvoTin 时间: 2025-3-11 12:12

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

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

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