计算化学公社

标题: 请问计算激发态T1、T2的能量相等该怎么解决呢？ [打印本页]

作者
Author: DFT.. 时间: 2026-3-10 10:49
标题: 请问计算激发态T1、T2的能量相等该怎么解决呢？
本帖最后由 DFT.. 于 2026-3-10 10:49 编辑

各位老师好，最近在复现文献数据(计算激发态)时遇到下列问题:
1.计算T1、T2态的能量相等，怀疑出现了势能面交叉，尝试缩小优化步长上限，还是没能解决问题，请问该怎么解决呢？
2.后续引入CB7做优化，应该会有比较大的耗时，请问有比较省时的优化策略吗？
3.Supporting Information中提到，引入CB7的体系采用了QM/MM方法，但Gaussian好像不支持QM/MM吧，可以使用哪个软件的呢？

附件包含S0、S1、T1、T2的输入输出文件及参考文献，麻烦老师帮忙检查一下是否存在问题，感谢各位老师解答。
(, 下载次数 Times of downloads: 11)

作者
Author: wal 时间: 2026-3-10 13:10
确实交叉了 T2漏到了T1
高斯只支持ONIOM，设置力场会比较费劲，可以考虑换软件或者用chemshell接口接别的程序做QM/MM
(, 下载次数 Times of downloads: 15)

作者
Author: DFT.. 时间: 2026-3-10 20:06

wal 发表于 2026-3-10 13:10
确实交叉了 T2漏到了T1
高斯只支持ONIOM，设置力场会比较费劲，可以考虑换软件或者用chemshell接口接别的 ...

感谢您的回复，请问对于这种“T2漏到T1”的情况，怎么处理才能得到正确的T2结构呢？

作者
Author: wal 时间: 2026-3-10 20:57
本帖最后由 wal 于 2026-3-10 21:07 编辑

DFT.. 发表于 2026-3-10 20:06
感谢您的回复，请问对于这种“T2漏到T1”的情况，怎么处理才能得到正确的T2结构呢？

一直监视，发现交叉了就停掉重新交，或者让脚本帮你做这件事
实话讲可以考虑放弃了，这种情况T2不一定真的有极小点。我今天看见你的帖子就用我自己写的的脚本帮你跑了一下，发现优化在T1/T2交叉点卡住了
(, 下载次数 Times of downloads: 10)
拿中间那个近似收敛的点做起始结构，调调步长优化或许能拿到一个root=2，几何收敛的点，但实际上这跟一般讲的T2极小点可能没啥关系
(, 下载次数 Times of downloads: 2)

（暴论，这种我倾向从S0激发到T2之后绕不开圆锥交叉点，总会掉到T1）

作者
Author: sobereva 时间: 2026-3-11 07:01
Gaussian可以做ONIOM(QM:MM)，一个麻烦之处是需要提供QM层原子的力场参数，但可以借助这个解决
借助sobtop给Gaussian便利地补ONIOM(QM:MM)计算牵扯的力场参数的重要程序BaneTop：http://bbs.keinsci.com/thread-53665-1-1.html

作者
Author: DFT.. 时间: 2026-3-11 10:58

wal 发表于 2026-3-10 20:57
一直监视，发现交叉了就停掉重新交，或者让脚本帮你做这件事
实话讲可以考虑放弃了，这种情况T2不一定真 ...

非常感谢您的回复@sobereva

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