计算化学公社

标题: cp2k-5.1 单点能计算中止的问题 [打印本页]
作者
Author: kdclark    时间: 2020-7-7 09:32
标题: cp2k-5.1 单点能计算中止的问题
大家好。我最近初学cp2k，编译完试了几个tests没有问题，然后就自己建了几个结构试了下单点能计算，大部分能算完，但是个别会在SCF中止。out文件内中止信息和程序运行错误信息如下：
  1. SCF WAVEFUNCTION OPTIMIZATION

  3.   Step     Update method      Time    Convergence         Total energy    Change
  4.   ------------------------------------------------------------------------------
  5.      1 P_Mix/Diag. 0.40E+00    0.6            NaN      -319.8326073128 -3.20E+02

  7. *******************************************************************************
  8. *   ___                                                                       *
  9. *  /   \                                                                      *
  10. * [ABORT]                                                                     *
  11. *  \___/                KS energy is an abnormal value (NaN/Inf).             *
  12. *    |                                                                        *
  13. *  O/|                                                                        *
  14. * /| |                                                                        *
  15. * / \                                                     qs_ks_methods.F:841 *
  16. *******************************************************************************
复制代码
  1. Backtrace for this error:
  2. #0  0x7fea8aa7a2ed in ???
  3. #1  0x7fea8aa79503 in ???
  4. #2  0x7fea8a0f6f1f in ???
  5. #3  0x7fea8a0f6e97 in ???
  6. #4  0x7fea8a0f8800 in ???
  7. #5  0x5606ea2b6648 in __message_passing_MOD_mp_abort
  8.         at /root/cp2k-5.1/src/mpiwrap/message_passing.F:1060
  9. #6  0x5606ea3140d4 in __base_hooks_MOD_cp_abort
  10.         at /root/cp2k-5.1/obj/local/sopt/base_hooks.F90:68
  11. #7  0x5606ea31429f in __base_hooks_MOD_cp__b
  12.         at /root/cp2k-5.1/obj/local/sopt/base_hooks.F90:155
  13. #8  0x5606e87a3599 in __qs_ks_methods_MOD_qs_ks_build_kohn_sham_matrix
  14.         at /root/cp2k-5.1/src/qs_ks_methods.F:841
  15. #9  0x5606e87a6c8d in rebuild_ks_matrix
  16.         at /root/cp2k-5.1/src/qs_ks_methods.F:1182
  17. #10  0x5606e87a6c8d in __qs_ks_methods_MOD_qs_ks_update_qs_env
  18.         at /root/cp2k-5.1/src/qs_ks_methods.F:1091
  19. #11  0x5606e8d29579 in __qs_scf_MOD_scf_env_do_scf
  20.         at /root/cp2k-5.1/src/qs_scf.F:438
  21. #12  0x5606e8d384f3 in __qs_scf_MOD_scf
  22.         at /root/cp2k-5.1/src/qs_scf.F:231
  23. #13  0x5606e88d7425 in __qs_energy_MOD_qs_energies
  24.         at /root/cp2k-5.1/src/qs_energy.F:85
  25. #14  0x5606e926fa89 in __qs_force_MOD_qs_calc_energy_force
  26.         at /root/cp2k-5.1/src/qs_force.F:118
  27. #15  0x5606e87d9f07 in __force_env_methods_MOD_force_env_calc_energy_force
  28.         at /root/cp2k-5.1/src/force_env_methods.F:239
  29. #16  0x5606e8444124 in cp2k_run
  30.         at /root/cp2k-5.1/src/start/cp2k_runs.F:321
  31. #17  0x5606e8445bac in __cp2k_runs_MOD_run_input
  32.         at /root/cp2k-5.1/src/start/cp2k_runs.F:1099
  33. #18  0x5606e843afb1 in cp2k
  34.         at /root/cp2k-5.1/src/start/cp2k.F:281
  35. #19  0x5606e83b463e in main
  36.         at /root/cp2k-5.1/obj/local/sopt/cp2k.F90:45
复制代码

请问这是什么原因导致的？avogadro的插件和cif2cell都试过，晶胞数据应该没问题。
另外，我是第一次手动编译的，在/usr/bin文件夹会有cp2k.popt的程序，后来又重新用脚本编译了一次，在/cp2k-5.1/exe/local文件夹得到cp2k.sopt的程序，但是其他几个没有出现，这是什么问题呢？不过两个程序都能进行计算，并且得出了相同的结果和错误。。。
下面附的是发生中止的pc.inp/out文件，以及两个可算完的inp文件。

(, 下载次数 Times of downloads: 36) (, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 12) (, 下载次数 Times of downloads: 7)
作者
Author: yoghurt117    时间: 2020-7-9 22:12
是不是版本的问题，我的机器上用6.1算就出内存故障，用5.1就没事，你换个试试。
作者
Author: kdclark    时间: 2020-7-10 08:45
yoghurt117 发表于 2020-7-9 22:12
是不是版本的问题，我的机器上用6.1算就出内存故障，用5.1就没事，你换个试试。

谢谢回复，那我换个版本换个机器编译下试试。
作者
Author: fineren    时间: 2020-7-10 11:16
7.1都出了大半年了，咋还用5呢
作者
Author: yoghurt117    时间: 2020-7-10 19:15
fineren 发表于 2020-7-10 11:16
7.1都出了大半年了，咋还用5呢

很多时候，软件这东西它并不是越新越好.......
作者
Author: kdclark    时间: 2020-8-20 11:15
本帖最后由 kdclark 于 2020-8-20 11:17 编辑

参考楼上两位回复，最近两天换了个机器，重新编译了cp2k-7.1，make tests全部通过了，然而在算pc.inp时还是出现了上面的情况，另外两个inp依然顺利算出，现在感觉也许是pc.inp文件的问题，但是换了数个输出inp的程序（avogadro，cp2k_Editor等），还是不清楚inp文件的问题出在哪。。。。有哪位大佬可以帮忙运行一下pc.inp，我想确定一下是不是它的问题。。。
作者
Author: sobereva    时间: 2021-2-10 04:01
lanl2dz赝势基组怎么能和GTH赝势搭配，且不说别的，先把这个致命问题解决了



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