IRC扫描的右端反应物能量，与opt中的能量不一致

yukinolike · 发表于 Post on 2023-3-3 15:21:29

本帖最后由 yukinolike 于 2023-3-3 15:25 编辑

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

chands · 发表于 Post on 2023-3-3 17:34:09

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

fatpig · 发表于 Post on 2023-3-3 18:07:10

用guess=read读取irc的波函数，或者用最后一点的结构先做stable=opt再做优化...

yukinolike · 发表于 Post on 2023-3-3 21:41:19

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

不错的建议，我再试试。

yukinolike · 发表于 Post on 2023-3-3 21:43:55

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

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

chands · 发表于 Post on 2023-3-3 22:00:14

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

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

我爱糖葫芦树 · 发表于 Post on 2023-7-14 10:01:32

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

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

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

[Gaussian/gview] IRC扫描的右端反应物能量，与opt中的能量不一致