IRC 过程中在 TS 附近出现能量尖峰与梯度异常，这种 IRC 是否可信？

HuiYr · 发表于 Post on 2025-12-24 17:45:28

本帖最后由 HuiYr 于 2025-12-24 17:49 编辑

各位老师好，目前我利用Gaussian计算一个金属催化的反应，

图一是最近搜索的一个过渡态的结果，有且仅有一个虚频，且是我想要的；图二是根据该结果跑IRC，利用#p IRC(calcfc,lqa,maxpoints=60) b3lyp/genecp nosymm empiricaldispersion=gd3bj关键词得到的结果。左右两边分别是TS前后的两个中间体，这样的IRC结果是我没有见过的，想要请教各位老师，① 这样的结果可信吗？② 得到这样的结果会是什么原因呢？③ 会是由于我没有进行guess=read的原因吗？

附：计算的是Fe(Ⅲ)的二重态，发现S**2=0.8787，偏离0.75>10%。④ 请问会是自旋污染导致这样的结果吗？

pikachuupup · 发表于 Post on 2025-12-24 18:32:51

检验一下TS的波函数稳定性，然后用guess=read读取稳定波函数，并用forward和reverse两个方向分开跑IRC

HuiYr · 发表于 Post on 2025-12-24 19:14:28

pikachuupup 发表于 2025-12-24 18:32
检验一下TS的波函数稳定性，然后用guess=read读取稳定波函数，并用forward和reverse两个方向分开跑IRC

#p stable=opt b3lyp/genecp nosymm guess=check empiricaldispersion=gd3bj geom=check 利用该关键词检查了波函数，发现波函数是稳定的。（The wavefunction is stable under the perturbations considered.
The wavefunction is already stable. ）那我重新读check，然后分开两边跑一下irc。谢谢你~

XM123 · 发表于 Post on 2026-7-25 12:26:39

pikachuupup 发表于 2025-12-24 18:32
检验一下TS的波函数稳定性，然后用guess=read读取稳定波函数，并用forward和reverse两个方向分开跑IRC

你好我是新手小白目前也遇到了相同的问题请问我要怎样检查波函数的稳定性呢？用了wavefunction关键词没有找到如果不问稳定加什么关键词修改呢？

李唯一 · 发表于 Post on 2026-7-25 15:43:34

XM123 发表于 2026-7-25 12:26
你好我是新手小白目前也遇到了相同的问题请问我要怎样检查波函数的稳定性呢？用了wavefunction关键词 ...

stable=opt

XM123 · 发表于 Post on 2026-7-27 15:47:45

李唯一发表于 2026-7-25 15:43
stable=opt

好的谢谢~

江月待何人 · 发表于 Post on 2026-8-9 20:01:45

确保IRC和优化在完全一样的水平下计算
3不是
4DFT理论上不考虑有自选污染

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

[Gaussian/gview] IRC 过程中在 TS 附近出现能量尖峰与梯度异常，这种 IRC 是否可信？

浏览过的版块