ORCA做CASSCF优化问题求助

尚艳磊

各位前辈们好，最近在学习用ORCA做CASSCF计算，遇到一些问题。
我有已经用DFT方法优化好的TS结构，但在用单参考态做单点能计算的时候这个TS的T1有点大，想做一个多参考态的验证计算。
选择（5,5）对应1个自由基轨道、C-N和O-O键的成键与反键轨道，高斯结构已附在结尾。
我首先在高斯里面用uhf/3-21g做了个uno初猜，然后在Multiwfn确认了轨道是我想要的顺序，接着在ORCA里面转换为了gbw文件。
为了保证在做CASSCF优化的时候轨道不会被旋转，我采用了如下的输入参数：
%pal nprocs 60 end
%maxcore 5000
%scf MaxIter 512 end
! 3-21g moread optTS NumFreq VeryTightSCF NoIter defgrid3 KDIIS
%moinp "TS-casscf-321g.gbw"
%scf Thresh 1e-12 Tcut 1e-14 end
%casscf nel 5 norb 5 end


但是总是优化不出来，有次终于5个YES达到结构优化收敛了，但是频率计算好像报错了
ORCA finished by error termination in CASSCF
Calling Command: /work/home/ac5qzwfvnn/orca503/orca_casscf TS-casscf_D00052.casinp.tmp >TS-casscf_D00052.lastscf
Check for CASSCF-logfiles
ORCA_NUMCALC: ORCA finished with an error in the energy calculation

ORCA finished by error termination in CASSCF
Calling Command: /work/home/ac5qzwfvnn/orca503/orca_casscf TS-casscf_D00026.casinp.tmp >TS-casscf_D00026.lastscf
Check for CASSCF-logfiles
ORCA_NUMCALC: ORCA finished with an error in the energy calculation
Please check TS-casscf.lastint.dpl52 for details
Please check TS-casscf.lastscf.dpl52 for details
Please check TS-casscf.lastint.dpl26 for details
Please check TS-casscf.lastscf.dpl26 for details
[file orca_numfreq/numfreq_utils.cpp, line 2106]: Skipping calculation . . .

[file orca_numfreq/numfreq_utils.cpp, line 2106]: Skipping calculation . .


查阅提示的文件，发现[file orca_casscf/orca_casscf.cpp, line 8170]: ERROR (ORCA_CASSCF): Convergence Failure. Retry with modified convergence settings!
似乎是SCF不收敛。然后在ORCA论坛发帖询问，得到答复是 it makes no sense to try a CASSCF optimization with noiter. You need to converge the wavefunction or you cannot get a proper gradient.

有两个疑问：1. FrankNeese教授说加noiter做casscf优化没有意义，那如果不加的话，轨道在优化过程会进行交换，这就不是我想要考察的轨道了呀。我也试了不加noiter，结果就是最终结构不是DFT下的TS结构，相差很多。2. 收敛的波函数，我试着用单点能读取uno.gbw得到了收敛的gbw，然后再优化还是不行呢？

感谢各位前辈们答疑解惑

hebrewsnabla

如果你的活性空间选得是合理的，那么在轨道优化之后是可以保持的。不能保持通常说明不合理。

当然不能用noiter来解决这个问题。这就好比dft的scf不收敛用(5/13=1)一样。

3-21g太小了。

wjc404

你这个体系(TS-casscf.inp)的成键方式有些怪异，我用B3LYP-D3/def2-TZVP算频率没有发现虚频。

zjxitcc

同意楼上所说，这个过渡态很奇怪，你检验过DFT波函数稳定性吗？如果有，并且稳定，可否说说你的DFT具体计算级别，比如晒出所有关键词，还有稳定性检测结果？