计算化学公社

标题: ORCA计算单点能报错问题 [打印本页]
作者
Author: 桃子量子    时间: 2023-11-23 09:51
标题: 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不收敛吗，还是什么原因，输入文件的坐标是用高斯优化和频率计算得到的结构坐标，

如果想要修改，应该怎么修改呢，
作者
Author: zjxitcc    时间: 2023-11-23 12:35
本帖最后由 zjxitcc 于 2023-11-23 12:36 编辑

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

  6. title

  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传轨道。


作者
Author: 桃子量子    时间: 2023-11-23 15:36
好的谢谢老师，我先去计算一下，有什么不懂再来请教，非常感谢您
作者
Author: wzkchem5    时间: 2023-11-23 16:58
另外也可以在orca输入文件里用NoTRAH关键字，好处是省机时和人力（省去了高斯计算以及转换轨道格式的时间）。当然这样收敛的概率不是100%，不过zjxitcc的方法的收敛概率也不是100%，高斯SCF收敛做得再好也不能保证收敛，说“一次必成”是过于夸张了。
作者
Author: 桃子量子    时间: 2023-11-23 19:09
wzkchem5 发表于 2023-11-23 16:58
另外也可以在orca输入文件里用NoTRAH关键字，好处是省机时和人力（省去了高斯计算以及转换轨道格式的时间） ...

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

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



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3