计算化学公社

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

作者
Author: BinWang 时间: 7 hour ago
标题: 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里面还有别的参数可以动一下吗？

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