orca SCF不收敛的问题求助

413 · 发表于 Post on 2021-6-7 13:42:19

请教各位，我算的这个体系SCF不收敛是什么鬼。

我试过拿高斯收敛的波函数，转成mkl，然后转成gbw。
依然会有下面的错误，请问还有啥办法吗？

Starting the XC term evaluation                   ... done ( 1.1 sec)
  promolecular density results
   # of electrons  = 3431.021613477
   EX             = -14510.607127539
   EC             = -200.622333729
   EX+EC          = -14711.229461269
Transforming the Hamiltonian                      ... done ( 0.3 sec)
Diagonalizing the Hamiltonian                   ... done ( 0.3 sec)
Back transforming the eigenvectors                ... done ( 0.1 sec)
Now organizing SCF variables                      ... done
                  ------------------
                  INITIAL GUESS DONE (  10.0 sec)
                  ------------------
--------------
SCF ITERATIONS
--------------
ITER    Energy       Delta-E       Max-DP    RMS-DP    [F,P]    Damp
            ***  Starting incremental Fock matrix formation  ***
--------------------------------------------------------------------------
mpirun noticed that process rank 22 with PID 202286 on node cn086 exited on signal 9 (Killed).
--------------------------------------------------------------------------

ORCA finished by error termination in SCF
Calling Command: mpirun -np 32  /apps/orca/4.0.1/orca_scf_mpi complex_a.gbw b
[file orca_tools/qcmsg.cpp, line 432]:
  .... aborting the run

zjxitcc · 发表于 Post on 2021-6-7 13:48:04

这不能称之为SCF不收敛，毕竟SCF刚开始，SCF表示锅从天上来。这看不出啥，如果没有更多信息，只能理解为内存不够了。如果你的分子很普通，可以晒出输入文件，大家一看便知。

413 · 发表于 Post on 2021-6-7 13:53:57

zjxitcc 发表于 2021-6-7 13:48
这不能称之为SCF不收敛，毕竟SCF刚开始，SCF表示锅从天上来。这看不出啥，如果没有更多信息，只能理解为内 ...

哦哦，我换台机器试下哈。可能是内存的问题。谢谢哈。

413 · 发表于 Post on 2021-6-7 14:04:34

zjxitcc 发表于 2021-6-7 13:48
这不能称之为SCF不收敛，毕竟SCF刚开始，SCF表示锅从天上来。这看不出啥，如果没有更多信息，只能理解为内 ...

谢谢哈。换了台机器。可以跑了。是内存的问题。

413 · 发表于 Post on 2021-6-7 16:08:37

zjxitcc 发表于 2021-6-7 13:48
这不能称之为SCF不收敛，毕竟SCF刚开始，SCF表示锅从天上来。这看不出啥，如果没有更多信息，只能理解为内 ...

再向您请教下哈。
跑倒是能跑了，但是看现在的能量变化，多半会挂掉。我拿了gaussian的输出作为初始波函数猜测。
请问您还有什么建议吗？

SCF ITERATIONS
--------------
ITER    Energy       Delta-E       Max-DP    RMS-DP    [F,P]    Damp
            ***  Starting incremental Fock matrix formation  ***
  0 -14870.5194251119 0.000000000000 7.45894027  0.01518066  0.2697914 0.7000
  1 -14869.5892037229 0.93022138898217.48903319  0.02698390  0.2014162 0.7000
  2 -14865.7531547905 3.836048932382551.96918934  1.08971846  0.3272674 0.7000
  3 -14786.9836098479  78.7695449425661485.71963037  1.86388274  0.4340102 0.7000
  4 -14508.5599429389 278.4236669090397328.27610196  6.97539154  0.8230833 0.7000
  5 -13939.9423390510 568.6176038879204732.12929817  5.14994026  0.9341669 0.7000
  6 -13808.0908455292 131.8514935217241724.24580802  2.78575626  0.7130251 0.7000
  7 -13594.8097015694 213.2811439598663045.05870251  3.68538799  0.8990258 0.7000
  8 -13470.9571112119 123.8525903574922121.35139021  2.87436759  0.9258706 0.7000
  9 -13416.4238013946  54.5333098172341900.34915797  2.55255051  1.0255426 0.7000
10 -13494.4884120933 -78.0646106986901760.77849361  2.63939581  1.0603534 0.7000
11 -13518.7037151593 -24.2153030659471804.74118212  2.56703102  1.0714884 0.7000
                           ***Turning on DIIS***
12 -13527.8565581072  -9.1528429479442253.73786697  3.02119061  1.1170250 0.7000
13 -13344.5715590649 183.2849990423641905.52268554  2.49518868  1.1779902 0.7000
14 -13906.3466220789-561.7750630140291339.06049015  1.71637034  0.7771690 0.7000
15 -14218.5674794504-312.220857371460955.27400838  1.18780957  0.5486193 0.7000
16 -14353.3315785032-134.764099052816671.23131504  0.83378691  0.5160707 0.7000
17 -14402.0019410173 -48.6703625140943836.13052244  3.88093765  0.6736166 0.7000
18 -13482.1802187001 919.8217223171382085.71590445  3.38774817  1.0689628 0.7000
19 -12061.82541332811420.3548053720151976.22272393  2.93709519  1.8874050 0.7000
            *** Restarting incremental Fock matrix formation ***
                                 *** Resetting DIIS ***
20 -12951.3871932088-889.5617798807293969.69736599  4.85857909  1.6198844 0.7000
21 -13893.7168566663-942.3296634574141519.16654839  2.01619719  0.6683754 0.7000
22 -14176.1456953484-282.428838682172952.77325218  1.71754041  0.4258142 0.7000
23 -14419.1127022608-242.967006912404752.71145059  1.32164788  0.3238083 0.7000
24 -14558.0730538593-138.960351598515534.20453476  0.94382225  0.2301301 0.7000
25 -14642.7019278165 -84.628873957205378.73912210  0.65916462  0.3311451 0.7000
26 -14696.4673779420 -53.765450125502207.94176526  0.53020474  0.4361619 0.7000
27 -14535.5029949404 160.964383001628313.12985224  0.40948330  0.7193302 0.7000

zjxitcc · 发表于 Post on 2021-6-7 16:12:57

413 发表于 2021-6-7 16:08
再向您请教下哈。
跑倒是能跑了，但是看现在的能量变化，多半会挂掉。我拿了gaussian的输出作为初始波函 ...

“拿了gaussian的输出作为初始波函数猜测”有点含糊，有的人高斯和ORCA里用的方法、基组都不一样，也自称用作初始猜测，这种操作导致不收敛的概率当然大一些。应当尽可能用一样的方法、基组，且要添加一些控制精度的关键词。

你得交待两个程序里你用的分别是什么方法、基组，关键词分别是什么，传轨道用的哪个程序。

413 · 发表于 Post on 2021-6-7 16:30:43

zjxitcc 发表于 2021-6-7 16:12
“拿了gaussian的输出作为初始波函数猜测”有点含糊，有的人高斯和ORCA里用的方法、基组都不一样，也自称 ...

泛函都是用的pbe
基组都是用的赝势基组lanl2dz
先从chk到fchk，然后从fchk到mkl，然后从mkl到gbw。
高斯的关键词：
#p opt=(maxcycle=512) pbepbe/genecp scf=(conver=6, maxcycle=512)
...
0 2
...

C 0
6-31G(d)
****
O 0
6-31G(d)
****
H 0
6-31G(d)
****
Au Pt 0
lanl2dz
****

Au Pt 0
lanl2dz

orca的关键词：
! opt pbe lanl2dz angs PRINTBASIS autoaux
%maxcore 5000
%scf MaxIter 1024 Convergence loose end
%pal nprocs 32 end
%basis
NewGTO H "6-31G(d)" end #
DelECP H
NewGTO O "6-31G(d)" end #
DelECP O
NewGTO C "6-31G(d)" end #
DelECP C
end
* xyz 0 2

zjxitcc · 发表于 Post on 2021-6-7 16:37:18

本帖最后由 zjxitcc 于 2021-6-7 16:40 编辑

413 发表于 2021-6-7 16:30
泛函都是用的pbe
基组都是用的赝势基组lanl2dz
先从chk到fchk，然后从fchk到mkl，然后从mkl到gbw。

你这离正确传轨道的操作还有不少距离。
（1）还是没讲传轨道用的哪个程序
（2）高斯关键词要加nosymm int=nobasistransform，而maxcycle=512, conver=6这些没用的就不用加了。
（3）应该在高斯中检验波函数稳定性，即加stable=opt，然后再传轨道。

413 · 发表于 Post on 2021-6-7 16:54:14

zjxitcc 发表于 2021-6-7 16:37
你这离正确传轨道的操作还有不少距离。
（1）还是没讲传轨道用的哪个程序
（2）高斯关键词要加nosymm int ...

我是用的Multiwfn_3.8_dev_bin_Win64传的轨道。

zjxitcc · 发表于 Post on 2021-6-7 16:59:22

413 发表于 2021-6-7 16:54
我是用的Multiwfn_3.8_dev_bin_Win64传的轨道。

无论用哪个，至少你得采纳我在8L说的两条建议，不然不是合理的传轨道方式。如果你重算之后还有问题，可以再贴上来。

sobereva · 发表于 Post on 2021-6-7 17:31:09

写scf=(conver=6没意义，既然是要传，何不传充分收敛的

当前计算最好用def2-SVP，对Au Pt本身就是赝势基组，在ORCA里还有标配的辅助基组，比autoaux产生的便宜不少。而且都用def2-SVP在写关键词的时候还省事很多

另外，Gaussian里都opt过了，放ORCA里还opt干嘛

angs、PRINTBASIS写这些都没意义

wzkchem5 · 发表于 Post on 2021-6-7 17:57:38

413 发表于 2021-6-7 09:08
再向您请教下哈。
跑倒是能跑了，但是看现在的能量变化，多半会挂掉。我拿了gaussian的输出作为初始波函 ...

首先，你试过直接用orca默认的初猜做计算吗？只有orca自己的所有收敛方法都不能收敛，才应该考虑从高斯传轨道。绝不能把从高斯传轨道当作默认的计算方法，也不应该认为它天生比用orca直接算要靠谱。

hdhxx123 · 发表于 Post on 2021-6-7 18:07:28

wzkchem5 发表于 2021-6-7 17:57
首先，你试过直接用orca默认的初猜做计算吗？只有orca自己的所有收敛方法都不能收敛，才应该考虑从高斯传 ...

如果闭壳层体系这类波函数稳定性相对不复杂的情况我赞同，若是周期表靠后的或者铁磁类难收敛开壳层体系用高斯传轨道为默认我觉得比较好

413 · 发表于 Post on 2021-6-7 19:06:19

zjxitcc 发表于 2021-6-7 16:59
无论用哪个，至少你得采纳我在8L说的两条建议，不然不是合理的传轨道方式。如果你重算之后还有问题，可以 ...

好的，谢谢哈。我在重新算单点，然后重新传。有问题再向您请教。

413 · 发表于 Post on 2021-6-7 19:08:11

sobereva 发表于 2021-6-7 17:31
写scf=(conver=6没意义，既然是要传，何不传充分收敛的

当前计算最好用def2-SVP，对Au Pt本身就是赝势基 ...

之前所有的结构都是用orca算的。全部重新用高斯算太麻烦了。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] orca SCF不收敛的问题求助

浏览过的版块