求助 含Ag、I体系计算SOC矩阵元基组的选取

是aweia

我这有个500+原子的结构需要计算SOC矩阵元我试过
! B3LYP/G 6-31G*  miniprint tightSCF
%basis
NewGTO Ag "lanl08 " end
NewGTO I "lanl08 " end
end

但是在SCF ITERATIONS部分有warning

1. 15 -42293.1522134138-2094.680332722717218.36087118  0.28978659  0.7628759 0.7000
   
2. 16 -43458.6658805500-1165.513667136176152.66908579  0.19573475  0.9169294 0.7000
   
3. 17 -43648.8117072105-190.145826660584811.57395473  0.65469091  1.0836630 0.7000
   
4. 18 -39575.67178405654073.139923153991313.27466043  0.38185399  1.9248090 0.7000
   
5. 19 -35226.38945507884349.282328977759395.04141083  0.45001905  1.6988564 0.7000
   
6.               *** Restarting incremental Fock matrix formation ***
   
7.                                   *** Resetting DIIS ***
   
8. 
   
9. WARNING: the maximum gradient error descreased on average only by a factor   1.0
   
10.          during the last 20 iterations
    

复制代码

再往下仍然会报错（不知道这部分是在干啥）

1. --------------------------------------------------------------------------------------------
   
2.    Iter.        energy            ||Error||_2        Shift      TRadius  Mac/Mic        Rej.
   
3. --------------------------------------------------------------------------------------------
   
4. WARNING: 4052355 diagonal Hessian elements are negative!
   
5.     0      -7092.573272558230     3.917999e+02                     0.400 (TRAH MAcro)   No
   
6. WARNING : negative HOMO - LUMO gap : -127.243686
   

复制代码

我认为是我基组选择的不对
但不知道选择哪些合适，
我测试过程中用3-21G的时候没有这样的报错，是否是只能用全电子基组，但3-21G应该不行吧，用什么基组比较合适呢，原子数太多也不能随便换，望各位前辈赐教，谢谢。

wzkchem5

这不是报错，只是警告，表示SCF收敛困难，但不说明结果不能用。
SOC一般不建议用赝势计算，算得不准。宁可用QM/MM或者ONIOM，把重要的原子描述得很精确（用全电子相对论方法）、不重要的原子用半经验描述，也不要整体用赝势来算。