低级基组的过渡态IRC结果正常，高级基组过渡态的IRC出现问题

gjf · 发表于 Post on 2022-12-2 10:29:33

本帖最后由 gjf 于 2022-12-2 10:32 编辑

给位老师，我在计算一个反应的过渡态结构，在计算的时候，使用TS方法进行计算。刚开始使用低级基组6-31g对初猜结构进行过渡态计算，并且计算出了结果，跑出来的IRC也很平滑，产物和反应物符合预期。低级基组过渡态计算参数

IRC计算参数

为了得到更高精度的过渡态结构，对低级基组计算的过渡态结构，采用了高级基组进行重新计算，也计算出了有且仅有一个虚频，虚频振动方向也符合成键规律，对这个过渡态跑IRC时出现了问题。
高级基组参数

IRC 计算参数

这个图就是高级基组计算的过渡态跑出来的IRC，红色圆圈是过渡态结构的能量，不是整个IRC能量的最高点，并且发现IRC曲线上还有很多小凸起。所以想问各位老师，为什么低级基组计算的过渡态跑出来的IRC没问题，以低级基组计算的过渡态为初猜结构，采用高级基组进行计算后的过渡态结构，IRC却出现了问题。

sai77 · 发表于 Post on 2022-12-2 12:40:18

6-31g不加极化函数，结果完全没法用
第二次算的说明过渡态定位不准确，按照社长的博文解决

在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400

gjf · 发表于 Post on 2022-12-2 14:26:49

sai77 发表于 2022-12-2 12:40
6-31g不加极化函数，结果完全没法用
第二次算的说明过渡态定位不准确，按照社长的博文解决

谢谢老师，那我去试试加上极化函数进行优化。还有个问题想问一下，我以6-31g计算出的过渡态结构为初猜结构，采用高级基组重新进行计算过渡态，这样得出的结果还是不可以采用的吗？

chands · 发表于 Post on 2022-12-2 14:55:37

gjf 发表于 2022-12-2 14:26
谢谢老师，那我去试试加上极化函数进行优化。还有个问题想问一下，我以6-31g计算出的过渡态结构为初猜结 ...

如果你的高级别能确认过渡态，结果可以用。
把IRC stepsize改小一点看看。

gjf · 发表于 Post on 2022-12-2 15:24:29

chands 发表于 2022-12-2 14:55
如果你的高级别能确认过渡态，结果可以用。
把IRC stepsize改小一点看看。

高级别基组计算的过渡态，有且仅有一个虚频，且虚频的振动方向符合成键规律，只是跑出来IRC如上图所示，出现了问题。那老师我降低步长，在跑一次IRC

gjf · 发表于 Post on 2022-12-5 15:55:26

chands 发表于 2022-12-2 14:55
如果你的高级别能确认过渡态，结果可以用。
把IRC stepsize改小一点看看。

老师，我把maxpoints从100改为默认值后，跑了一次IRC，发现整个曲线就很正常，把我搞糊涂了，这个过渡态结构可以用吗？

chands · 发表于 Post on 2022-12-5 16:52:22

gjf 发表于 2022-12-5 15:55
老师，我把maxpoints从100改为默认值后，跑了一次IRC，发现整个曲线就很正常，把我搞糊涂了，这个过渡态 ...

应该可以。

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

[新手求助] 低级基组的过渡态IRC结果正常，高级基组过渡态的IRC出现问题