计算化学公社

标题: Orca算频率报错 [打印本页]

作者
Author: 413 时间: 2022-8-13 08:59
标题: Orca算频率报错
请问，下面的错误是咋回事呢？

IBatch 15 (of 20)
   CP-SCF ITERATION 0: 3.3814e-01 ( 138.3 sec 0/ 27 done)
   CP-SCF ITERATION 1: 1.0756e-01 ( 156.2 sec 0/ 27 done)
   CP-SCF ITERATION 2: 4.6905e-02 ( 157.7 sec 0/ 27 done)
   CP-SCF ITERATION 3: 1.7326e-02 ( 158.6 sec 0/ 27 done)
   CP-SCF ITERATION 4: 9.0690e-03 ( 160.5 sec 0/ 27 done)
   CP-SCF ITERATION 5: 5.5159e-03 ( 162.8 sec 0/ 27 done)
   CP-SCF ITERATION 6: 2.0631e-03 ( 164.1 sec  11/ 27 done)
   CP-SCF ITERATION 7: 9.2866e-04 ( 101.7 sec  27/ 27 done)
                  *** THE CP-SCF HAS CONVERGED ***
IBatch 16 (of 20)
   CP-SCF ITERATION 0: [file orca_cpscf/cpscf2.cpp, line 1087, Process 4]: Error (CP-SCF): Failed to read trial vector 15,0

ORCA finished by error termination in SCF Hessian

作者
Author: wzkchem5 时间: 2022-8-13 15:24
感觉像是I/O问题，比如硬盘有坏扇区

作者
Author: 413 时间: 2022-8-13 15:36

wzkchem5 发表于 2022-8-13 15:24
感觉像是I/O问题，比如硬盘有坏扇区

看着很像硬盘的问题

作者
Author: 413 时间: 2022-8-17 10:03

wzkchem5 发表于 2022-8-13 15:24
感觉像是I/O问题，比如硬盘有坏扇区

请问我算频率的时候，又碰到了一个新的错误，第一圈的scf已经收敛了，为啥还会报这种cp-scf的错误？

                  *** THE CP-SCF HAS CONVERGED ***
IBatch 3 (of 24)
   CP-SCF ITERATION 0: 1.0686e+11 ( 490.7 sec 0/ 21 done)
   CP-SCF ITERATION 1: 3.1932e+10 ( 281.3 sec 0/ 21 done)
   CP-SCF ITERATION 2: 1.1976e+10 ( 288.3 sec 0/ 21 done)
   CP-SCF ITERATION 3: 8.6535e+09 ( 295.5 sec 0/ 21 done)
......
   CP-SCF ITERATION 126: 2.6698e-03 ( 996.7 sec  19/ 21 done)
   CP-SCF ITERATION 127: 2.6715e-03 (1041.5 sec  19/ 21 done)

                        ***********
                        *E R R O R*
      ******************       ******************
      *** THE CP-SCF CALCULATION IS UNCONVERGED ***
      *********************************************

Warning : ERROR CODE RETURNED FROM CP-SCF PROGRAM
Codes    : res=-1,cmd=RUN_CPSCF
Consequence: aborting the Z-vector calculation

[file orca_scfhess/scfhess2.cpp, line 2272, Process 0]: Error (ORCA_SCFHESS): Solution of the CP-SCF equations failed.
[file orca_scfhess/scfhess2.cpp, line 2272, Process 14]: Error (ORCA_SCFHESS): Solution of the CP-SCF equations failed.
[file orca_scfhess/scfhess2.cpp, line 2272, Process 3]: Error (ORCA_SCFHESS): Solution of the CP-SCF equations failed.
......

作者
Author: abin 时间: 2022-8-17 10:58

413 发表于 2022-8-17 10:03
请问我算频率的时候，又碰到了一个新的错误，第一圈的scf已经收敛了，为啥还会报这种cp-scf的错误？

...

THE CP-SCF CALCULATION IS UNCONVERGED

作者
Author: wzkchem5 时间: 2022-8-17 18:37

413 发表于 2022-8-17 03:03
请问我算频率的时候，又碰到了一个新的错误，第一圈的scf已经收敛了，为啥还会报这种cp-scf的错误？

...

SCF和CP-SCF是两个完全不一样的东西，SCF收敛不代表CP-SCF也一定收敛。
但是两者的不收敛原因有一定的共通性，像你这种几乎收敛但没有完全收敛的情况，一般是因为数值误差问题，比如积分阈值太松、格点太小。这个对SCF、CP-SCF都是成立的。

作者
Author: 413 时间: 2022-8-17 21:28

wzkchem5 发表于 2022-8-17 18:37
SCF和CP-SCF是两个完全不一样的东西，SCF收敛不代表CP-SCF也一定收敛。
但是两者的不收敛原因有一定的共 ...

DEFGRID3
会有帮助吗

作者
Author: wzkchem5 时间: 2022-8-18 00:11

413 发表于 2022-8-17 14:28
DEFGRID3
会有帮助吗

有可能可以

作者
Author: 413 时间: 2022-8-18 21:22

wzkchem5 发表于 2022-8-18 00:11
有可能可以

我先试下

作者
Author: 413 时间: 2022-8-21 10:27
本帖最后由 413 于 2022-8-21 10:35 编辑

wzkchem5 发表于 2022-8-18 00:11
有可能可以

请问我在算频率时，碰到的CP-SCF的那个错误，跑到128圈没收敛就报错了。
我查了手册，可以将最大循环步数调大。
于是我在输入中加入了
%elprop
MaxIter 512 # max. number of iterations in CPSCF
end

可是为什么CP-SCF还是只能跑到128步呢？

作者
Author: wzkchem5 时间: 2022-8-21 14:49

413 发表于 2022-8-21 03:27
请问我在算频率时，碰到的CP-SCF的那个错误，跑到128圈没收敛就报错了。
我查了手册，可以将最大循环步 ...

因为%elprop里的maxiter只负责electric properties计算里的CP-SCF迭代的最大迭代次数，不负责其他计算里的CP-SCF迭代的最大迭代次数。
CP-SCF是很多计算都会用到的，并不单单是频率计算会用到

作者
Author: 413 时间: 2022-8-21 15:19

wzkchem5 发表于 2022-8-21 14:49
因为%elprop里的maxiter只负责electric properties计算里的CP-SCF迭代的最大迭代次数，不负责其他计算里 ...

那就是这个命令吗？
%Method
Z_MaxIter 512
end

作者
Author: wzkchem5 时间: 2022-8-21 16:24

413 发表于 2022-8-21 08:19
那就是这个命令吗？
%Method
Z_MaxIter 512

这种试一下就知道，不需要问人

作者
Author: 413 时间: 2022-8-21 17:13

wzkchem5 发表于 2022-8-21 16:24
这种试一下就知道，不需要问人

谢谢，我在试了，就是计算量好大

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