计算化学公社

标题: ORCA进行频率计算,PROPERTY INTEGRAL和SCF RESPONSE计算速度问题 [打印本页]

作者
Author:
SDAlsj    时间: 2026-2-7 13:12
标题: ORCA进行频率计算,PROPERTY INTEGRAL和SCF RESPONSE计算速度问题
本帖最后由 SDAlsj 于 2026-2-7 13:17 编辑

用ORCA进行400-420原子(只含C、H、O、N)的频率计算,计算级别是RI-BLYP-D3BJ/def2-SVP,CPCM溶剂模型,gCP考虑BSSE,有几个问题想请教下各位老师:
1.我发现计算耗时主要是在PROPERTY INTEGRAL和SCF RESPONSE上,想问下这两个模块是在计算什么,Hessain矩阵是在SCF RESPONSE里计算吗?
2.PROPERTY INTEGRAL是用batch来进行分解的。核数多但是单核内存少,虽然每个batch算的比较快,但是每个batch里包含的原子数也比较少;相反,核数少但单核内存多,虽然每个batch算的慢,但是batch里包含的原子数比较多。在这一部分上,核数多但单核内存少和核数少单核内存多计算速度反而差不多,想问下怎么调整核数和单核内存是最合理的呢?
3.SCF RESPONSE是否核越多算的越快呢?

4.读取几何优化(几何优化和频率计算严格在同一级别下进行)任务结束后的波函数到频率计算里能否加速PROPERTY INTEGRAL和SCF RESPONSE的速度呢?
5.我用16核并行,每个核14400MB内存(总内存为256GB),在计算到SCF RESPONSE模块上时出现ORCA finished by error termination in SCF RESPONSE [file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run 错误,不知道是每个核内存太少了,还是所用内存总量超过了总内存量的80%问题
谢谢各位老师指点迷津

作者
Author:
wzkchem5    时间: 2026-2-7 14:13
1. 简单来说,property integral计算的输出是SCF response的输入,SCF response的输出经过少许进一步的操作后得到Hessian
2. 在保持nprocs*maxcore基本不变、且小于机器总可用物理内存的条件下,若batch数多,则减少nprocs
3. SCF response也涉及分batch的问题。固定maxcore增加nprocs必然算得更快,但可能因为内存不够用而爆掉;保持nprocs*maxcore不变的情况下增加nprocs的话,很难说算得快还是慢
4. 不能,只能加速SCF迭代的速度
5. 这个报错没有任何信息量,必须给出输出文件最后至少几十行(不含排队系统输出),以及完整的stderr输出
作者
Author:
SDAlsj    时间: 2026-2-7 14:45
wzkchem5 发表于 2026-2-7 14:13
1. 简单来说,property integral计算的输出是SCF response的输入,SCF response的输出经过少许进一步的操作 ...

好的,谢谢老师,明白了
5.的问题应该还是单核爆内存了,它是在算了一天算完property integral后进入SCF Response时报错的,我现在用8核,增大单核内存量,没有报错,一切正常,以后我也会注意这种报错问题的提问方式。
我还有个小问题,解析频率的计算是不支持续算对嘛?比方说我property integral算完了,在SCF Response这个地方报错了,是没办法让它从SCF Response这里续算的对嘛
作者
Author:
wzkchem5    时间: 2026-2-7 20:49
SDAlsj 发表于 2026-2-7 14:45
好的,谢谢老师,明白了
5.的问题应该还是单核爆内存了,它是在算了一天算完property integral后进入SCF ...

对,只支持续算SCF部分(注意不是SCF Response部分)。因为SCF Response涉及的中间量如果要存硬盘的话,要写的内容太多,需要的读写时间太长,为SCF Response增加断点续算功能会拖慢计算
作者
Author:
SDAlsj    时间: 2026-2-7 22:20
wzkchem5 发表于 2026-2-7 20:49
对,只支持续算SCF部分(注意不是SCF Response部分)。因为SCF Response涉及的中间量如果要存硬盘的话, ...

好的,明白了,谢谢老师




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3