计算化学公社

标题: IRC扫描的右端反应物能量，与opt中的能量不一致 [打印本页]

作者
Author: yukinolike 时间: 2023-3-3 15:21
标题: IRC扫描的右端反应物能量，与opt中的能量不一致
本帖最后由 yukinolike 于 2023-3-3 15:25 编辑

如图，我做完的一个单方向的IRC扫描中，右端为反应物能量是-3552.2050 hartree左右，但是拿irc右端的结构去做opt，优化到的反应物能量却是-3552.1700左右。能量差了有0.03 hartree, 明明用到的基组和计算方法都是一样，为什么会出现这种差异呢？另外这个反应物优化到的能量甚至比过渡态能量还高...

作者
Author: chands 时间: 2023-3-3 17:34
IRC不是优化。IRC每一步的步长是固定的，而优化到最后步长必定变小直至收敛，因此IRC两端的结构不一定是极小点，必须将IRC两端的结构提取出来进一步优化已获得反应物（产物、中间体、反应前/后结合物）的结构。

作者
Author: fatpig 时间: 2023-3-3 18:07
用guess=read读取irc的波函数，或者用最后一点的结构先做stable=opt再做优化...

作者
Author: yukinolike 时间: 2023-3-3 21:41

fatpig 发表于 2023-3-3 18:07
用guess=read读取irc的波函数，或者用最后一点的结构先做stable=opt再做优化...

不错的建议，我再试试。

作者
Author: yukinolike 时间: 2023-3-3 21:43

chands 发表于 2023-3-3 17:34
IRC不是优化。IRC每一步的步长是固定的，而优化到最后步长必定变小直至收敛，因此IRC两端的结构不一定是极 ...

我知道这两个计算本身是不一样的，但是做完的irc，右端结构肯定是在极小点附近的。我最近卡在这里一两周了，拿右端结构去优化，不管怎么opt优化得到的结构，都是能量比过渡态能量高。

作者
Author: chands 时间: 2023-3-3 22:00

yukinolike 发表于 2023-3-3 21:43
我知道这两个计算本身是不一样的，但是做完的irc，右端结构肯定是在极小点附近的。我最近卡在这里一两周 ...

3楼的波函数稳定性可以试一下。我不知道具体体系不好说，理论级别是否严格一致？格点是否一致，有没有加溶剂模型，按理不至于差这么大。优化出来的结构跟原来IRC右端的结构差不多吗？还有一种可能是反应前复合物，不过差了80 kJ/mol有点大。

作者
Author: 我爱糖葫芦树 时间: 2023-7-14 10:01

yukinolike 发表于 2023-3-3 21:43
我知道这两个计算本身是不一样的，但是做完的irc，右端结构肯定是在极小点附近的。我最近卡在这里一两周 ...

请问贴主这个问题解决了吗，我的也是产物比ts能量高，irc拽出来那个点优化后也是这个样子

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