计算化学公社

标题: 关于ORCA的DLPNO-CCSD(T)方法占用巨量磁盘和存储不足报错 [打印本页]

作者
Author:
DoorBell    时间: 2019-3-27 15:01
标题: 关于ORCA的DLPNO-CCSD(T)方法占用巨量磁盘和存储不足报错
老师们好!最近在算一个钴的团簇,已经用DFT完成了粗优化,之后又用ORCA 4.1.1的DLPNO-CCSD(T)计算单点能。

但是无论在PC上还是超算上都会在这一步产生几百G的.tmp文件:
  1. ...
  2. --------------------------------
  3. INITIAL GUESS AND PNO GENERATION
  4. --------------------------------

  5. PNO truncation parameters                   ....
  6.     PAOOverlapThresh   =        1.000e-08

  7.     TCutPairs          =        1.000e-04
  8.     TCutPNO            =        3.330e-07
  9.     TCutPNOSingles     =        9.990e-09
  10.     TCutMP2Pairs       =        1.000e-05
  11.     TCutMKN            =        1.000e-03
  12.     TCutDO             =        1.000e-02

  13. Pair selection                          .... not used
  14. Type of local MP2 treatment                 .... semi-local MP2
  15. Strategy for PNO selection                  .... occupation number selection
  16. Pair density normalization                  .... MP2 norm
  17. Spin component scaling                      .... not used

  18.    .... Finished loop over pairs
  19. ...
复制代码
之后进行到此处就会报错OUT OF MEMORY!:
  1. ...
  2.    .... Finished loop over pairs

  3. PNO Occupation Number Statistics:
  4.   | Av. % of trace(Dij) retained        ...  95.447125430886
  5.   | sigma^2 in % of trace(Dij) retained ...   2.27e+01
  6.   | Av. % of trace(Di) retained         ...  99.997798967622
  7.   | sigma^2 in % of trace(Di) retained  ...   2.96e-06
  8. Distributions of % trace(Dij) recovered:
  9.   |       >= 99.9  ...   293 (  9.3 % of all pairs)
  10.   |   [60.0, 70.0) ...     1 (  0.0 % of all pairs)
  11.   |   [70.0, 80.0) ...    42 (  1.3 % of all pairs)
  12.   |   [80.0, 90.0) ...   304 (  9.7 % of all pairs)
  13.   |   [90.0, 99.0) ...  1905 ( 60.6 % of all pairs)
  14.   |   [99.0, 99.9) ...   600 ( 19.1 % of all pairs)
  15. Distributions of % trace(Di) recovered :
  16.   |       >= 99.9  ...    84 (100.0 % of all I-pairs )

  17. Making pair pair interaction lists          ... done

  18.                        ===========================
  19.                        1740 OF 3145 PAIRS ARE KEPT
  20.                        ===========================

  21. Total time spent in the initial guess                 ...  8761.445 sec
  22. SL-MP2 correlation energy (all non-screened pairs)    ...      -7.092094379993 Eh

  23. Initial PNO correlation energy                        ...      -7.003322263097 Eh
  24. sum of pair energies estimated for screened out pairs ...      -0.002482968226 Eh
  25. sum of MP2 pair energies for pairs that were not kept ...      -0.054621149365 Eh
  26. sum of PNO error estimates for the kept pairs         ...      -0.034236781988 Eh
  27.                                                           --------------------
  28. sum of all corrections                                         -0.091340899579
  29. Initial total correlation energy                               -7.094663162676
  30. Thresholds for map construction and integral transformation for fine MP2 and CCSD(T) calculation:
  31.    TCutMKN   ...     1.0e-03
  32.    TCutDO    ...     1.0e-02
  33.    TCutPairs ...     1.0e-04
  34.    TCutCMO   ...     1.0e-03
  35.    TCutCPAO  ...     1.0e-03

  36. Time for aux screen maps:            0.098
  37. Time for maps after fine MP2:            0.360
  38. -----------------------------
  39. LOCAL RI TRANSFORMATION (IJV)
  40. -----------------------------

  41. Orbital window:      40 to 123
  42. Basis functions:    920 (360 shells)
  43. Aux. functions:    5940 (1660 shells)

  44. Processing maps (0.0 sec)
  45. Average map sizes:
  46.    Aux shells -> MOs(i)           82.1
  47.    Aux shells -> MOs(j)           84.0
  48.    MOs        -> AO shells       349.5

  49. Calculating integrals (20.8 sec, 316.439 MB)
  50. Sorting integrals (1.0 sec, 316.372 MB)
  51. Total time for the integral transformation: 21.9 sec
  52. --------------------------------
  53. LOCAL RI TRANSFORMATION (VABPAO)
  54. --------------------------------

  55. Number of PAOs:     920
  56. Basis functions:    920 (360 shells)
  57. Aux. functions:    5940 (1660 shells)

  58. Processing maps (0.0 sec)
  59. Average map sizes:
  60.    Aux shells -> PAOs            920.0
  61.    PAOs       -> AO shells       360.0

  62. Calculating integrals (1139.8 sec, 69876.943 MB)
  63. Finished


  64. -------------------------------------
  65. Pair Pair Term precalculation with
  66. RI-(ij|mn) and (im|jn) transformation
  67. ON THE FLY
  68. -------------------------------------

  69.    IBatch   1 (of   1)              ...
  70. [file orca_tools/qcmem.cpp, line 918, Process 3]:  OUT OF MEMORY ERROR!

  71. [file orca_tools/qcmem.cpp, line 918, Process 3]:  OUT OF MEMORY ERROR!

  72. [file orca_tools/qcmem.cpp, line 918, Process 1]:  OUT OF M

  73. ...乱码...

  74. [file orca_tools/qcmem.cpp, line 918, Process 1]:  OUT OF MEMORY ERROR!

  75. cmem.cpp, line 918, Process 2]:  OUT OF MEMORY ERROR!

  76. [file orca_tools/qcmem.cpp, line 918, Process 2]:  OUT OF MEMORY ERROR!

  77. [file orca_tools/q

  78. ...乱码...

  79. [file orca_tools/qcmem.cpp, line 918, Process 5]:  OUT OF MEMORY ERROR!

  80. ]:  OUT OF MEMORY ERROR!

  81. [file orca_tools/qcmem.cpp, line 918, Process 4]:  OUT OF MEMORY ERROR!

  82. [file orca_tools/qcmem.cpp, line 918, Process 0]:  OUT OF MEMORY ERROR!

  83. [file orca_tools/qcmem.cpp, line 918, Process 0]:  OUT OF MEMORY ERROR!


  84. job aborted:
  85. [ranks] message

  86. [0-5] process exited without calling finalize

  87. ---- error analysis -----

  88. [0-5] on DESKTOP-KEETL36
  89. D:\ProgramFiles\ORCA\orca_mdci_mpi.exe ended prematurely and may have crashed. exit code 12

  90. ---- error analysis -----

  91. ORCA finished by error termination in MDCI
  92. Calling Command: mpiexec -np 6  D:\ProgramFiles\ORCA\orca_mdci_mpi.exe F:\ORCA\test.mdciinp.tmp F:\ORCA\test
  93. [file orca_tools/qcmsg.cpp, line 458]:
  94.   .... aborting the run

  95. [file orca_tools/qcmsg.cpp, line 458]:
  96.   .... aborting the run
复制代码
输入文件:
  1. ! RKS SlowConv NODIIS Grid6 DLPNO-CCSD(T) RI def2-TZVP AutoAux

  2. %scf
  3.     maxIter 1500
  4.     MaxIntMem 40960
  5.     MaxDisk 40960
  6. end

  7. %maxCore 1024

  8. %pal
  9.     nProcs 40
  10. end

  11. %output
  12.     XYZFile true
  13. end

  14. * xyz 0 1
  15. ...
  16. *
复制代码
原子簇如图所示。

由于不知道哪些信息是有用的,所以就多复制了一些文字,还请各位老师原谅。

非常感谢!




作者
Author:
DoorBell    时间: 2019-3-27 15:04
产生的.tmp文件是这样的:
作者
Author:
plus    时间: 2019-3-27 16:46
%maxCore 10000
并行核数*10G<总内存
作者
Author:
量化小菜鸡    时间: 2019-3-28 01:19
这种计算特别费内存,10个G也未必够。你试着增大%MaxCore减小nprocs看看。
作者
Author:
DoorBell    时间: 2019-3-28 13:21
量化小菜鸡 发表于 2019-3-28 01:19
这种计算特别费内存,10个G也未必够。你试着增大%MaxCore减小nprocs看看。

已经降低了并行核数并提高了内存,但是仍然产生了几百G的文件...一共就只有500G的存储......





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