计算化学公社

标题: 低级基组的过渡态IRC结果正常，高级基组过渡态的IRC出现问题 [打印本页]

作者
Author: gjf 时间: 2022-12-2 10:29
标题: 低级基组的过渡态IRC结果正常，高级基组过渡态的IRC出现问题
本帖最后由 gjf 于 2022-12-2 10:32 编辑

给位老师，我在计算一个反应的过渡态结构，在计算的时候，使用TS方法进行计算。刚开始使用低级基组6-31g对初猜结构进行过渡态计算，并且计算出了结果，跑出来的IRC也很平滑，产物和反应物符合预期。低级基组过渡态计算参数
(, 下载次数 Times of downloads: 13)
IRC计算参数
(, 下载次数 Times of downloads: 17)
为了得到更高精度的过渡态结构，对低级基组计算的过渡态结构，采用了高级基组进行重新计算，也计算出了有且仅有一个虚频，虚频振动方向也符合成键规律，对这个过渡态跑IRC时出现了问题。
高级基组参数
(, 下载次数 Times of downloads: 7)
IRC 计算参数
(, 下载次数 Times of downloads: 10)
这个图就是高级基组计算的过渡态跑出来的IRC，红色圆圈是过渡态结构的能量，不是整个IRC能量的最高点，并且发现IRC曲线上还有很多小凸起。所以想问各位老师，为什么低级基组计算的过渡态跑出来的IRC没问题，以低级基组计算的过渡态为初猜结构，采用高级基组进行计算后的过渡态结构，IRC却出现了问题。
(, 下载次数 Times of downloads: 10)

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

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

作者
Author: gjf 时间: 2022-12-2 14:26

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

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

作者
Author: chands 时间: 2022-12-2 14:55

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

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

作者
Author: gjf 时间: 2022-12-2 15:24

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

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

作者
Author: gjf 时间: 2022-12-5 15:55

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

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

作者
Author: chands 时间: 2022-12-5 16:52

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

应该可以。

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