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

嘤嘤嘤 · 发表于 Post on 2024-6-10 18:50:04

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

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

sobereva · 发表于 Post on 2021-7-1 05:48:06

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

这种结论不是普适的

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

wzkchem5 · 发表于 Post on 2021-6-27 21:38:53

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

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

hdhxx123 · 发表于 Post on 2021-6-27 20:49:36

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

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

LJF · 发表于 Post on 2021-6-25 15:14:43

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

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

exity · 发表于 Post on 2021-6-9 16:46:19

orca手册专门有一章说SCF收敛的

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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写这些都没意义

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

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

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

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

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

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

浏览过的版块