SobMECP查找双自由基的s0/t1交叉点失败

BinWang

各位老师：
       大家好！我有一个双自由基体系。最稳定的三重态结构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里面还有别的参数可以动一下吗？

sobereva

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

BinWang

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

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