在Gaussian中计算IRC的方法和常见问题

chands · 发表于 Post on 2022-9-16 20:00:31

Lym 发表于 2022-9-16 19:37
谢谢老师的回答，能垒是很小，但是最奇怪的是，最高点处不是过渡态。并且从图上来看，生成物(右边的点)的 ...

最高点处不是过渡态，可能跟IRC的累积误差有关，毕竟你的能垒太小了，步长一大，累积误差就可能导致最高点不在过渡态处。

你是说用低级别算生成物的能量（单点能）低于反应物，用高级别算生成物的能量高于反应物？有这个可能，看你的图生成物就比反应物低10 kJ/mol左右，一般情况下高级别算得准一些，可能IRC的级别太低，考虑更高的级别找过渡态和IRC。

你的反应势能面比较平缓，可能是无势垒反应，就不能用过渡态理论了。

Lym · 发表于 Post on 2022-9-16 20:13:38

chands 发表于 2022-9-16 20:00
最高点处不是过渡态，可能跟IRC的累积误差有关，毕竟你的能垒太小了，步长一大，累积误差就可能导致最高 ...

谢谢老师的回答
1.我从IRC跑出来的路径中，提取反应物，生成物以及TS的结构采用PBE0-D3的方式进行优化，然后再在M062X-D3/DEF2TZVP隐形溶剂模型下计算单点能。然后进行自由能的比较的时候发现反应物的能量为:-1980.732 a.u.; TS:-1980.730 a.u.; 生成物的能量为:-1980.726 a.u.

那我还需要在更高的计算方法下找到TS的结构->相同的方向扫描IRC->提取三个点进行优化嘛？

jaywang7 · 发表于 Post on 2023-11-6 12:28:07

看博文中说找过渡态与跑irc的计算级别必须严格一致，想问下我找过渡态结构优化时用了tight recalc=3 那么是不是跑irc的时候也要用tight recalc=3这样？

Freeman · 发表于 Post on 2023-11-6 15:58:42

jaywang7 发表于 2023-11-6 12:28
看博文中说找过渡态与跑irc的计算级别必须严格一致，想问下我找过渡态结构优化时用了tight recalc=3 那么是 ...

那倒不用。计算级别指泛函/方法、基组、格点、溶剂模型，其他细节（如收敛限）只要合理就行。

jaywang7 · 发表于 Post on 2023-11-6 20:24:25

Freeman 发表于 2023-11-6 15:58
那倒不用。计算级别指泛函/方法、基组、格点、溶剂模型，其他细节（如收敛限）只要合理就行。

okok

sobereva · 发表于 Post on 2023-11-7 06:23:24

jaywang7 发表于 2023-11-6 12:28
看博文中说找过渡态与跑irc的计算级别必须严格一致，想问下我找过渡态结构优化时用了tight recalc=3 那么是 ...

这里说的“级别”是指对能量有系统性影响的因素，跟tight、recalc没关系

IRC里面的tight跟几何优化的tight更是根本没直接关系

jaywang7 · 发表于 Post on 2023-11-7 13:01:35

sobereva 发表于 2023-11-7 06:23
这里说的“级别”是指对能量有系统性影响的因素，跟tight、recalc没关系

IRC里面的tight跟几何优化的t ...

好的谢谢sob老师

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

[Gaussian/gview] 在Gaussian中计算IRC的方法和常见问题