计算化学公社

标题: SobMECP查找双自由基的s0/t1交叉点失败 [打印本页]

作者
Author: BinWang 时间: 2026-3-6 10:31
标题: SobMECP查找双自由基的s0/t1交叉点失败
本帖最后由 BinWang 于 2026-3-6 10:33 编辑

各位老师：
   大家好！我有一个双自由基体系。最稳定的三重态结构T0的能量比最稳定的单重态S1的能量低0.5 kcal/mol。最稳定的单重态计算三重态能量，也是三重态能量更低。两者的结构仅仅相差0.38 Anstrong。
    以单重态结构作为初猜，我利用SobMECP查找两者的交叉点，发现第一步最接近，其余跑了100步，跑偏了。这是第一步的数据：
   ProgFile Exists - OK 2  First Input OK
  3
  4  Max Gradient El.: 0.000596 (0.000700)  YES
  5  RMS Gradient El.: 0.000060 (0.000500)  YES
  6  Max Change of X: 0.000417 (0.004000)  YES
  7  RMS Change of X: 0.000042 (0.002500)  YES
  8  Difference in E: 0.000893 (0.000050) NO
我利用过渡态作为初猜，最好的两步是：
  Max Gradient El.: 0.001269 (0.000700) NO
RMS Gradient El.: 0.000278 (0.000500)  YES
Max Change of X: 0.003802 (0.004000)  YES
RMS Change of X: 0.000862 (0.002500)  YES
Difference in E: 0.000039 (0.000050)  YES

Step Number 75 -- MECP not yet converged

Max Gradient El.: 0.000893 (0.000700) NO
RMS Gradient El.: 0.000209 (0.000500)  YES
Max Change of X: 0.006891 (0.004000) NO
RMS Change of X: 0.001551 (0.002500)  YES
Difference in E: 0.000014 (0.000050)  YES
我的高斯输入文件是：
%mem=20GB
%nproc=16
%chk=singlet.chk
#n  uwb97xd/genecp nosym force guess(read) SCF=DIIS

First State

1 1

。。。。。坐标省略

Ir 0
SDD
****
C N O H S 0
6-31G**
****

Ir 0
SDD

这种情况，各位老师有什么建议？SobMECP程序里面的STPMX = 0.1d0，我减小了5倍，10倍，50倍，都不起作用。
有没有可以用多参考方法查找mecp的程序？或者我目前得到的某个结构放松收敛限制其实就可以用了？这三个结构都比较接近收敛限，该选哪一个？sobMECP里面还有别的参数可以动一下吗？

作者
Author: sobereva 时间: 2026-3-7 06:06
两个自旋态之间大概率本来就没有交叉点。T0和S1极小点下都是T0能量最低，而且两个极小点结构又十分接近，很大程度已经说明了这一点。

作者
Author: BinWang 时间: 2026-3-7 09:05

sobereva 发表于 2026-3-7 06:06
两个自旋态之间大概率本来就没有交叉点。T0和S1极小点下都是T0能量最低，而且两个极小点结构又十分接近，很 ...

非常感谢Sob老师的回复！
两个稳定T0/S1的双自由基结构之间没有交叉点，这个是确定的。这个双自由基是三元环开环得到的，我用开环的过渡态作为初猜也没有找到交叉点，这个有什么办法？因为过渡态三重态的能量是要高于单重态的。
我一楼的帖子上面贴出来了利用过渡态找mecp的最接近的两个点，这两个点距离mecp的收敛标准就差一点点。

作者
Author: sobereva 时间: 2026-3-8 09:21

BinWang 发表于 2026-3-7 09:05
非常感谢Sob老师的回复！
两个稳定T0/S1的双自由基结构之间没有交叉点，这个是确定的。这个双自由 ...

在三重态反应的IRC上计算各个点的单重态能量，理应会有三重态和单重态能量正好近乎简并的点，或许是更好的初猜

作者
Author: BinWang 时间: 2026-3-9 09:22

sobereva 发表于 2026-3-8 09:21
在三重态反应的IRC上计算各个点的单重态能量，理应会有三重态和单重态能量正好近乎简并的点，或许是更好 ...

非常感谢Sob老师！

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