计算化学公社

标题: orca SCF不收敛的问题求助 [打印本页]

作者
Author: 413 时间: 2021-6-7 13:42
标题: orca SCF不收敛的问题求助
请教各位，我算的这个体系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

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

作者
Author: 413 时间: 2021-6-7 13:53

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

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

作者
Author: 413 时间: 2021-6-7 14:04

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

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

作者
Author: 413 时间: 2021-6-7 16:08

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

作者
Author: zjxitcc 时间: 2021-6-7 16:12

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

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

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

作者
Author: 413 时间: 2021-6-7 16:30

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

作者
Author: zjxitcc 时间: 2021-6-7 16:37
本帖最后由 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，然后再传轨道。

作者
Author: 413 时间: 2021-6-7 16:54

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

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

作者
Author: zjxitcc 时间: 2021-6-7 16:59

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

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

作者
Author: sobereva 时间: 2021-6-7 17:31
写scf=(conver=6没意义，既然是要传，何不传充分收敛的

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

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

angs、PRINTBASIS写这些都没意义

作者
Author: wzkchem5 时间: 2021-6-7 17:57

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

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

作者
Author: hdhxx123 时间: 2021-6-7 18:07

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

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

作者
Author: 413 时间: 2021-6-7 19:06

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

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

作者
Author: 413 时间: 2021-6-7 19:08

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

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

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

作者
Author: 413 时间: 2021-6-7 19:15

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

我试过orca默认的初猜，也是存在不收敛的情况。

作者
Author: exity 时间: 2021-6-9 16:46
orca手册专门有一章说SCF收敛的

作者
Author: LJF 时间: 2021-6-25 15:14

hdhxx123 发表于 2021-6-7 18:07
如果闭壳层体系这类波函数稳定性相对不复杂的情况我赞同，若是周期表靠后的或者铁磁类难收敛开壳层体系用 ...

我同意。本人的经验，对相对论效应比较明显的元素来说，用ORCA的全电子基组和Gaussian的赝势都能收敛的情况下，往往也是用赝势会准确些。

作者
Author: hdhxx123 时间: 2021-6-27 20:49

LJF 发表于 2021-6-25 15:14
我同意。本人的经验，对相对论效应比较明显的元素来说，用ORCA的全电子基组和Gaussian的赝势都能收敛的情 ...

你可能误会了我的意思，我强调的是Gaussian的波函数稳定性包括收敛性方面做得比ORCA好，并不涉及基组方面的讨论

作者
Author: wzkchem5 时间: 2021-6-27 21:38

LJF 发表于 2021-6-25 08:14
我同意。本人的经验，对相对论效应比较明显的元素来说，用ORCA的全电子基组和Gaussian的赝势都能收敛的情 ...

我的经验是，一般赝势的准确度还是不如ZORA+相对论性全电子基组。我遇到过有的时候赝势和ZORA/全电子基组算出来的键长的差别高达0.05埃，对键能影响达到20kcal/mol，而且都是ZORA/全电子基组和实验对得更好。而且从原理上，ZORA/全电子基组也比赝势要严格得多。
你确定你发现赝势更准确，不是因为误差抵消的缘故吗

作者
Author: sobereva 时间: 2021-7-1 05:48

LJF 发表于 2021-6-25 15:14
我同意。本人的经验，对相对论效应比较明显的元素来说，用ORCA的全电子基组和Gaussian的赝势都能收敛的情 ...

这种结论不是普适的

什么相对论哈密顿、什么赝势、什么全电子基组、什么赝势基组，细节都得十分明确才有可比性

作者
Author: 嘤嘤嘤 时间: 2024-6-10 18:50

413 发表于 2021-6-7 14:04
谢谢哈。换了台机器。可以跑了。是内存的问题。

请问您内存给了多少啊？多少核数？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)