计算化学公社

标题: 求助：优化过渡态用不同基组单点能会发生很大变化 [打印本页]

作者
Author: zym01 时间: 2020-10-27 09:25
标题: 求助：优化过渡态用不同基组单点能会发生很大变化
本帖最后由 zym01 于 2020-10-27 09:53 编辑

用的方法和基组分别是m062x 6-311++g(d,p) 和 m062x 6-311++g(3df,2p) ,都添加了溶剂模型scrf=(smd,solvent=water)。同一个过渡态，首先我在低基组进行过渡态优化，然后进行高基组的过渡态优化，两个吉布斯自由能差值达到了21kcal，经过排查我发现是单点能出现了较大的变化，请问这样正常吗？高基组的单点能和预期不符，过渡态的能量比反应物还低，但是低基组的单点能比较正常，是否可以直接用低基组算出的单点能？ (, 下载次数 Times of downloads: 7) （另一个文件没法上传，被限制了）

作者
Author: 喵星大佬 时间: 2020-10-27 10:08
哪有比不同级别计算的过渡态的绝对能量的。。。。。

作者
Author: zym01 时间: 2020-10-27 10:24
本帖最后由 zym01 于 2020-10-27 10:28 编辑

喵星大佬发表于 2020-10-27 10:08
哪有比不同级别计算的过渡态的绝对能量的。。。。。

因为我要用高基组算单点，然后发现过渡态的吉布斯自由能比反应物的更低，所以就比较了一下，做过IRC确定没有其他的中间体了，想知道该怎么办，谢谢。

作者
Author: 喵星大佬 时间: 2020-10-27 10:28
本帖最后由喵星大佬于 2020-10-27 10:29 编辑

zym01 发表于 2020-10-27 10:24
因为我要用高基组算单点，然后发现过渡态的吉布斯自由能比反应物的更低，所以就比较了一下，做过IRC确定 ...

反应物和过渡态都在相同级别优化取热力学校正量然后再在相同的更高级别算得电子能量？

还有，6-311++g(3df,2p)这种基组就别用了，def2-TZVPD或者ma-TZVP就好了
基组变大了绝对能量降低不是应该的么

作者
Author: zym01 时间: 2020-10-27 10:35
本帖最后由 zym01 于 2020-10-27 10:36 编辑

喵星大佬发表于 2020-10-27 10:28
反应物和过渡态都在相同级别优化取热力学校正量然后再在相同的更高级别算得电子能量？

还有，6-311++g ...

1：是的
2：经过高基组优化出来反应物的单点能升高，不知道为什么会这样

作者
Author: 喵星大佬 时间: 2020-10-27 10:45

zym01 发表于 2020-10-27 10:35
1：是的
2：经过高基组优化出来反应物的单点能升高，不知道为什么会这样

检查一下波函数稳定性，可能没有收敛到稳定的波函数

作者
Author: zym01 时间: 2020-10-27 10:49

喵星大佬发表于 2020-10-27 10:45
检查一下波函数稳定性，可能没有收敛到稳定的波函数

好的，谢谢

作者
Author: 火星哥 时间: 2023-8-2 15:37

喵星大佬发表于 2020-10-27 10:28
反应物和过渡态都在相同级别优化取热力学校正量然后再在相同的更高级别算得电子能量？

还有，6-311++g ...

我发现XAT与SN2的单点在低基组和高基组下能量也是有差异，大小关系反转。目前只能靠势阱来说明XAT是优势，要不计算结果与实验结论不符

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