计算化学公社

标题: Ag团簇吸附COOH的频率计算一直失败 [打印本页]

作者
Author: XxxXOO 时间: 2025-5-1 18:51
标题: Ag团簇吸附COOH的频率计算一直失败
本帖最后由 XxxXOO 于 2025-5-11 13:24 编辑

我的体系是58个Ag原子，吸附COOH，使用的基组和泛函是 PBE0 D4  def2-SV(P)
我试了8 ，16，32核内存从10000MB增加到100000MB都成功不了
-------------------------------------------------------------------------------
                           ORCA SCF HESSIAN
-------------------------------------------------------------------------------

               --------- SHARK INITIALIZATION DONE ---------

Hessian of the Kohn-Sham DFT energy:
Kohn-Sham wavefunction type                   ... UKS
Hartree-Fock exchange scaling                   ... 0.250
Number of operators                            ... 2
Number of atoms                               ... 62
Basis set dimensions                            ... 1436
Integral neglect threshold                      ... 2.5e-11
Integral primitive cutoff                      ... 2.5e-12

Setting up DFT Hessian calculations             ...

HESS GRID
---------

General Integration Accuracy    IntAcc    ... 4.629
Radial Grid Type                RadialGrid  ... OptM3 with GC (2021)
Angular Grid (max. ang.)       AngularGrid ... 5 (Lebedev-434)
Angular grid pruning method    GridPruning ... 4 (adaptive)
Weight generation scheme       WeightScheme... Becke
Basis function cutoff          BFCut    ... 1.0000e-11
Integration weight cutoff       WCut       ... 1.0000e-14
Angular grids for H and He will be reduced by one unit
Partially contracted basis set             ... off
Rotationally invariant grid construction    ... off

Total number of grid points                ... 933432
Total number of batches                   ... 14614
Average number of points per batch          ...    63
Average number of grid points per atom    ... 15055

Building densities                            ... done (    2.5 sec)
Calculating rho(r) on the grid                ... done (    14.1 sec)
Building xc-kernel on the grid                ... done (    0.6 sec)
                                                   done (    36.0 sec)

Nuclear repulsion Hessian  (SHARK)             ... done (    0.0 sec)

----------------------------------------------
Forming right-hand sides of CP-SCF equations    ...
----------------------------------------------
One electron integral derivatives  (SHARK)    ... done (    6.2 sec)
Transforming the overlap derivative matrices    ... done (    17.1 sec)
Making the Q(x) pseudodensities                ... done (    15.9 sec)
Adding the E*S(x)*S(y) terms to the Hessian    ... done (    17.6 sec)
Calculating energy weighted overlap derivatives  ... done (    2.3 sec)
F(y): ECP integral derivatives                ... done ( 2647.1 sec)
Two electron integral derivatives (SHARK-RI)    ...
ORCA finished by error termination in SCF Hessian
Calling Command: mpirun -np 16  /work/home/acr8708ots/WHT/orca/orca_scfhess_mpi FreqCOOH.scfhess.inp FreqCOOH
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run
都会出现这个问题，怎么解决呢？

作者
Author: wzkchem5 时间: 2025-5-1 19:36
内存设的不是越大越好。以下两种情况都会导致无法正常计算：
（1）计算所需内存大于设定的内存；
（2）设定的内存大于机器可用物理内存。
不要只考虑到（1）的情况而忽视（2）的可能。
此外，还需要留意会不会是硬盘空间不够

作者
Author: XxxXOO 时间: 2025-5-7 10:22

wzkchem5 发表于 2025-5-1 19:36
内存设的不是越大越好。以下两种情况都会导致无法正常计算：
（1）计算所需内存大于设定的内存；
（2）设 ...

您好，我上传了输出文件，您有时间看一下，万分感谢了

作者
Author: wzkchem5 时间: 2025-5-7 10:58

XxxXOO 发表于 2025-5-7 10:22
您好，我上传了输出文件，您有时间看一下，万分感谢了

我在私信里已经跟你说了，orca不支持在加smeartemp的情况下算频率。
你应该去掉smeartemp，再跑，然后再看是否报错，如果报错，再上传输出文件。

作者
Author: XxxXOO 时间: 2025-5-7 11:11

wzkchem5 发表于 2025-5-7 10:58
我在私信里已经跟你说了，orca不支持在加smeartemp的情况下算频率。
你应该去掉smeartemp，再跑，然后再 ...

好的好的，谢谢谢谢

作者
Author: XxxXOO 时间: 2025-5-11 13:25

wzkchem5 发表于 2025-5-1 19:36
内存设的不是越大越好。以下两种情况都会导致无法正常计算：
（1）计算所需内存大于设定的内存；
（2）设 ...

我试了一下那个方法，还是有报错，输出文件我上传了，有时间可以帮我看一下吗？谢谢

作者
Author: wzkchem5 时间: 2025-5-12 09:25

XxxXOO 发表于 2025-5-11 13:25
我试了一下那个方法，还是有报错，输出文件我上传了，有时间可以帮我看一下吗？谢谢

检查maxcore是不是太小了

作者
Author: XxxXOO 时间: 2025-5-12 09:35

wzkchem5 发表于 2025-5-12 09:25
检查maxcore是不是太小了

我设置的1000，应该不小吧

作者
Author: XxxXOO 时间: 2025-5-12 09:56

wzkchem5 发表于 2025-5-12 09:25
检查maxcore是不是太小了

在ORCA中计算频率时，可以固定原子吗？

作者
Author: wzkchem5 时间: 2025-5-12 10:24

XxxXOO 发表于 2025-5-12 09:35
我设置的1000，应该不小吧

1000很小，默认值可能就在1000左右。需要特意增大maxcore的场合，一般maxcore至少设4000

作者
Author: wzkchem5 时间: 2025-5-12 10:24

XxxXOO 发表于 2025-5-12 09:56
在ORCA中计算频率时，可以固定原子吗？

算数值频率时可以，具体看手册数值频率相关章节

作者
Author: XxxXOO 时间: 2025-5-13 13:06

wzkchem5 发表于 2025-5-12 10:24
1000很小，默认值可能就在1000左右。需要特意增大maxcore的场合，一般maxcore至少设4000

好的好的，我再试一下

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