计算化学公社

标题: IRC每次跑到0.8A时就报错不收敛 [打印本页]

作者
Author: goodgood 时间: 2025-4-7 20:51
标题: IRC每次跑到0.8A时就报错不收敛
本帖最后由 goodgood 于 2025-5-13 09:52 编辑

各位老师好，我最近在计算反应的过渡态，目前过渡态已经找到，有且仅有一个虚频，振动方向也是没问题的。
我之前在M06-2X/6-311+G**级别上找到了过渡态，并且IRC都没有问题，过渡态虚频为1010。
后来准备用更高精度revDSD-PBEP86-D3/jun-cc-pV(T+d)Z的试一试，过渡态成功找到，但是虚频竟然高达2370。于是我打算跑一个IRC试试，正向跑没问题，但是每次反向跑的时候，计算把步长调小，但是只要到0.8A这里能量陡增，不收敛报错。尝试用scf(xqc)就能成功结束。

我想请问：
（1）revDSD-PBEP86-D3/jun-cc-pV(T+d)Z比M06-2X/6-311+G**优化得到的过渡态虚频相差这么多是合理的吗？选择的泛函和基组是否有不合理的地方
（2）为什么每次在那一步的时候就能量激增，是否是我找的过渡态不合理

作者
Author: Stardust0831 时间: 2025-4-7 22:16
本帖最后由 Stardust0831 于 2025-4-7 22:17 编辑

既然想用高精度的理论方法，不如试试让高斯联用orca《将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分析》
我最推荐跑RI-wB97M-D4/ma-def2-QZVP级别的TS和IRC，这个级别有解析梯度和解析hessian，精度也是杂化泛函的天花板（orca关键词记得加defgrid3和verytightscf）。
想拿到更好的精度可以用RI-revDSD-PBEP86-D4/ma-def2-QZVPP甚至CCSD(T)-F12/cc-pVTZ-F12 with RI，在orca里都有解析梯度，可以得到极高精度的计算结果。
上面提到的三个级别可靠性更高，同时得益于RI加速，耗时比不带RI大大降低。

能量陡增可能是因为SCF还没有收敛，或者收敛到了不稳定的波函数。我总建议在跑完TS以后检查一下波函数稳定性。

作者
Author: sobereva 时间: 2025-4-8 06:39
一些MP2成份较高的双杂化泛函算一些开壳层体系未必结果理想，频率可能糟糕，见J. Phys. Chem. A, 124, 6899 (2020)。所以不要盲目认为用revDSD-PBEP86-D3结果就会更好。可以结合wB97XD/def2-TZVPP进行验证，以及尝试MP2成份低的双杂化泛函如PBE0-DH（MP2成份12.5%）

作者
Author: goodgood 时间: 2025-4-8 19:17

Stardust0831 发表于 2025-4-7 22:16
既然想用高精度的理论方法，不如试试让高斯联用orca《将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分 ...

谢谢您的回答。我检查了ts优化后的波函数是稳定的。后来我加了scf=xqc进行二次收敛后，IRC正常结束了，并且也检查了二次收敛后的波函数也是稳定的。请问这样找到的过渡态可用吗，会不会被审稿人质疑

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