ORCA计算单点能报错问题

桃子量子

输入文件
! DLPNO-CCSD(T) normalPNO RIJK aug-cc-pVTZ aug-cc-pVTZ/JK aug-cc-pVTZ/C tightSCF noautostart miniprint nopop
%maxcore  3000
%pal nprocs   10 end
%cpcm
smd true
SMDsolvent "water"
end
* xyz  -2   2
C                 -1.38548700    0.51532300   -1.02463400
C                 -0.43977600   -0.13433600    0.01104200
C                  0.95426500    0.52308200    0.06569800
S                  2.22338200   -0.52951600    0.93321800
O                  3.20837200    0.42079200    1.42303000
O                  1.43856700   -1.16953200    1.98206600
O                  2.68961800   -1.42756800   -0.11227500
F                 -0.32794000   -1.42888400   -0.32171900
F                 -1.02962800    0.18681500   -2.25736300
F                 -1.03300400   -0.03611300    1.20651100
F                  1.33646100    0.77535100   -1.19949800
F                 -1.36932900    1.83726000   -0.90673100
F                  0.82706600    1.68910400    0.72195000
F                 -2.62362500    0.08388000   -0.81486800
O                  4.44596500    3.10050300   -1.77667300
H                  3.82969200    2.70104500   -1.09034900
H                  4.26395500    2.53776100   -2.55360900
*

输出文件后面的报错
WARNING: the maximum gradient error descreased on average only by a factor   1.1
          during the last 20 iterations


                      *** Initiating the TRAH-SCF procedure ***


--------------------------------------------------------------------------------------------
   Iter.        energy            ||Error||_2        Shift      TRadius  Mac/Mic        Rej.
--------------------------------------------------------------------------------------------
    0      -1508.422820536128     2.646621e-03                     0.400 (TRAH MAcro)   No


job aborted:
[ranks] message


[0-9] process exited without calling finalize


---- error analysis -----


[0-9] on PC-202305051251
D:\orca\orca_scf_mpi.exe ended prematurely and may have crashed. exit code 62


---- error analysis -----


ORCA finished by error termination in SCF
Calling Command: mpiexec -np 10  D:\orca\orca_scf_mpi.exe E:\sjn\C3F7SO3-1H2O-AB\A-1\C3F7SO3-1H2O-a-1-E.gbw b E:\sjn\C3F7SO3-1H2O-AB\A-1\C3F7SO3-1H2O-a-1-E
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run
使用ORCA计算单点能报错了，想请教一下报错是因为SCF不收敛吗，还是什么原因，输入文件的坐标是用高斯优化和频率计算得到的结构坐标，

如果想要修改，应该怎么修改呢，

zjxitcc

SCF收敛较难，进入了TRAH步骤。解决这个问题的办法不止一种，我这里提供一种“一次必成”的方法（因为多次尝试各种关键词是需要人工和机时成本的）。我们可以先用Gaussian做个UHF计算，gjf文件如下

1. %chk=test.chk
   
2. %mem=128GB
   
3. %nprocshared=64
   
4. #p UHF/aug-cc-pVTZ nosymm int=nobasistransform scrf=smd
   
5. 
   
6. title
   
7. 
   
8. -2 2
   
9. C                 -1.38548700    0.51532300   -1.02463400
   
10. C                 -0.43977600   -0.13433600    0.01104200
    
11. C                  0.95426500    0.52308200    0.06569800
    
12. S                  2.22338200   -0.52951600    0.93321800
    
13. O                  3.20837200    0.42079200    1.42303000
    
14. O                  1.43856700   -1.16953200    1.98206600
    
15. O                  2.68961800   -1.42756800   -0.11227500
    
16. F                 -0.32794000   -1.42888400   -0.32171900
    
17. F                 -1.02962800    0.18681500   -2.25736300
    
18. F                 -1.03300400   -0.03611300    1.20651100
    
19. F                  1.33646100    0.77535100   -1.19949800
    
20. F                 -1.36932900    1.83726000   -0.90673100
    
21. F                  0.82706600    1.68910400    0.72195000
    
22. F                 -2.62362500    0.08388000   -0.81486800
    
23. O                  4.44596500    3.10050300   -1.77667300
    
24. H                  3.82969200    2.70104500   -1.09034900
    
25. H                  4.26395500    2.53776100   -2.55360900
    

复制代码

完成后依次运行

1. fch2mkl test.chk
   
2. orca_2mkl test_o -gbw

复制代码

即产生ORCA输入文件test_o.inp，波函数文件test_o.gbw。里面写好了基组信息、坐标，收敛的UHF轨道和一些简单关键词。打开test_o.inp文件将前几行修改为

1. %pal nprocs 18 end
   
2. %maxcore 10000
   
3. ! UHF TightSCF DLPNO-CCSD(T) normalPNO RIJK aug-cc-pVTZ/JK aug-cc-pVTZ/C
   
4. %cpcm
   
5.   smd true
   
6.   SMDsolvent "water"
   
7. end
   
8. %scf
   
9.   Thresh 1e-12
   
10.   Tcut 1e-14
    
11.   CNVDamp False
    
12. end

复制代码

提交给ORCA，SCF将在11圈内顺利收敛，进入DLPNO-CCSD(T)计算步骤，不用操心SCF问题。

注1：这里的内存和核数是演示，请根据你的机器实际情况进行修改。
注2：传轨道后没有立即收敛，而是迭代了11圈，是因为高斯的SMD与ORCA的SMD实现细节不一样；高斯没有RIJK，而ORCA开了RIJK。
注3：fch2mkl是开源免费的MOKIT中的一个小程序，用于Gaussian->ORCA传轨道。

桃子量子

好的谢谢老师，我先去计算一下，有什么不懂再来请教，非常感谢您

wzkchem5

另外也可以在orca输入文件里用NoTRAH关键字，好处是省机时和人力（省去了高斯计算以及转换轨道格式的时间）。当然这样收敛的概率不是100%，不过zjxitcc的方法的收敛概率也不是100%，高斯SCF收敛做得再好也不能保证收敛，说“一次必成”是过于夸张了。

桃子量子

wzkchem5 发表于 2023-11-23 16:58
另外也可以在orca输入文件里用NoTRAH关键字，好处是省机时和人力（省去了高斯计算以及转换轨道格式的时间） ...

好的，谢谢您

桃子量子

zjxitcc 发表于 2023-11-23 12:35
SCF收敛较难，进入了TRAH步骤。解决这个问题的办法不止一种，我这里提供一种“一次必成”的方法（因为多次 ...

经过您的指点，结果已经成功算出啦，非常感谢