计算化学公社

标题: IRC结果与单独优化结果不一致 [打印本页]

作者
Author: S321037017 时间: 2023-6-5 17:49
标题: IRC结果与单独优化结果不一致
1、找到了过渡态结构，用过渡态结果设置了IRC，得到的IRC结果把反应物和产物的结构复制出来设置相同的计算方法进行优化，结果发现IRC上显示的过渡态和反应物的自由能和和单独优化的过渡态和反应物的结果相差很大，不知道什么原因。（反应物优化没有虚频，过渡态只有一个虚频）
IRC：irc=(calcfc,maxpoints=30,lqa) b3lyp scrf=(smd,solvent=water) nosymm def2tzvp em=gd3
TS:opt=(calcfc,ts,noeigen,maxstep=100) freq b3lyp scrf=(smd,solvent=water) nosymm def2tzvp em=gd3
反应物优化：opt freq b3lyp scrf=(smd,solvent=water) nosymm def2tzvp em=gd3
IRC反应物和过渡态结构以及自由能：
(, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 4)
单独优化的反应物和过渡态：
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 4)
2、还有个问题就是IRC反应物的结构复制出来进行单独优化，得到的结果显是分子偏移了，想问下单独优化之后的结果能不能用啊，影响大不大？
IRC结构：
(, 下载次数 Times of downloads: 3)
单独优化的结构：

(, 下载次数 Times of downloads: 4)
希望有大佬可以解答，谢谢大家！

作者
Author: hebrewsnabla 时间: 2023-6-5 19:13
本帖最后由 hebrewsnabla 于 2023-6-5 19:16 编辑

1. 你的前两张图是电子能量，只能和单独优化结果的电子能量比较。

2. 反应物、生成物的单独优化结果和IRC末端有少许区别是正常的，区别太大不行。

作者
Author: S321037017 时间: 2023-6-5 20:40

hebrewsnabla 发表于 2023-6-5 19:13
1. 你的前两张图是电子能量，只能和单独优化结果的电子能量比较。

2. 反应物、生成物的单独优化结果和IR ...

好的，谢谢

作者
Author: sobereva 时间: 2023-6-6 10:23
没有可比性。IRC上根本就不可能有准确的反应物和产物结构，因为IRC步长是固定的，这点在此文里都已经充分强调了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

IRC路径两边端点的结构没有任何实际化学意义，用处仅限于判断过渡态找没找对、作为进一步优化极小点的初猜

而且哪怕maxpoints设得非常大尽量让IRC跑得完整，也有可能过早就满足Gaussian的IRC任务的结束判据而自动停掉，而后面距离极小点实际上可能还很远。

作者
Author: S321037017 时间: 2023-6-6 14:28

sobereva 发表于 2023-6-6 10:23
没有可比性。IRC上根本就不可能有准确的反应物和产物结构，因为IRC步长是固定的，这点在此文里都已经充分强 ...

好的，非常感谢卢老师

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