cp2k-5.1 单点能计算中止的问题

kdclark · 发表于 Post on 2020-7-7 09:32:15

大家好。我最近初学cp2k，编译完试了几个tests没有问题，然后就自己建了几个结构试了下单点能计算，大部分能算完，但是个别会在SCF中止。out文件内中止信息和程序运行错误信息如下：

SCF WAVEFUNCTION OPTIMIZATION
Step Update method Time Convergence Total energy Change
------------------------------------------------------------------------------
1 P_Mix/Diag. 0.40E+00 0.6 NaN -319.8326073128 -3.20E+02
*******************************************************************************
* ___ *
* / \ *
* [ABORT] *
* \___/ KS energy is an abnormal value (NaN/Inf). *
* | *
* O/| *
* /| | *
* / \ qs_ks_methods.F:841 *
*******************************************************************************

复制代码

Backtrace for this error:
#0 0x7fea8aa7a2ed in ???
#1 0x7fea8aa79503 in ???
#2 0x7fea8a0f6f1f in ???
#3 0x7fea8a0f6e97 in ???
#4 0x7fea8a0f8800 in ???
#5 0x5606ea2b6648 in __message_passing_MOD_mp_abort
at /root/cp2k-5.1/src/mpiwrap/message_passing.F:1060
#6 0x5606ea3140d4 in __base_hooks_MOD_cp_abort
at /root/cp2k-5.1/obj/local/sopt/base_hooks.F90:68
#7 0x5606ea31429f in __base_hooks_MOD_cp__b
at /root/cp2k-5.1/obj/local/sopt/base_hooks.F90:155
#8 0x5606e87a3599 in __qs_ks_methods_MOD_qs_ks_build_kohn_sham_matrix
at /root/cp2k-5.1/src/qs_ks_methods.F:841
#9 0x5606e87a6c8d in rebuild_ks_matrix
at /root/cp2k-5.1/src/qs_ks_methods.F:1182
#10 0x5606e87a6c8d in __qs_ks_methods_MOD_qs_ks_update_qs_env
at /root/cp2k-5.1/src/qs_ks_methods.F:1091
#11 0x5606e8d29579 in __qs_scf_MOD_scf_env_do_scf
at /root/cp2k-5.1/src/qs_scf.F:438
#12 0x5606e8d384f3 in __qs_scf_MOD_scf
at /root/cp2k-5.1/src/qs_scf.F:231
#13 0x5606e88d7425 in __qs_energy_MOD_qs_energies
at /root/cp2k-5.1/src/qs_energy.F:85
#14 0x5606e926fa89 in __qs_force_MOD_qs_calc_energy_force
at /root/cp2k-5.1/src/qs_force.F:118
#15 0x5606e87d9f07 in __force_env_methods_MOD_force_env_calc_energy_force
at /root/cp2k-5.1/src/force_env_methods.F:239
#16 0x5606e8444124 in cp2k_run
at /root/cp2k-5.1/src/start/cp2k_runs.F:321
#17 0x5606e8445bac in __cp2k_runs_MOD_run_input
at /root/cp2k-5.1/src/start/cp2k_runs.F:1099
#18 0x5606e843afb1 in cp2k
at /root/cp2k-5.1/src/start/cp2k.F:281
#19 0x5606e83b463e in main
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文件。

pc.inp (2.76 KB, 下载次数 Times of downloads: 36)

pc.out (188.66 KB, 下载次数 Times of downloads: 4)

hhtp.inp (4.99 KB, 下载次数 Times of downloads: 12)

py.inp (2.49 KB, 下载次数 Times of downloads: 7)

yoghurt117 · 发表于 Post on 2020-7-9 22:12:49

是不是版本的问题，我的机器上用6.1算就出内存故障，用5.1就没事，你换个试试。

kdclark · 发表于 Post on 2020-7-10 08:45:25

yoghurt117 发表于 2020-7-9 22:12
是不是版本的问题，我的机器上用6.1算就出内存故障，用5.1就没事，你换个试试。

谢谢回复，那我换个版本换个机器编译下试试。

fineren · 发表于 Post on 2020-7-10 11:16:55

7.1都出了大半年了，咋还用5呢

yoghurt117 · 发表于 Post on 2020-7-10 19:15:02

fineren 发表于 2020-7-10 11:16
7.1都出了大半年了，咋还用5呢

很多时候，软件这东西它并不是越新越好.......

kdclark · 发表于 Post on 2020-8-20 11:15:39

本帖最后由 kdclark 于 2020-8-20 11:17 编辑

参考楼上两位回复，最近两天换了个机器，重新编译了cp2k-7.1，make tests全部通过了，然而在算pc.inp时还是出现了上面的情况，另外两个inp依然顺利算出，现在感觉也许是pc.inp文件的问题，但是换了数个输出inp的程序（avogadro，cp2k_Editor等），还是不清楚inp文件的问题出在哪。。。。有哪位大佬可以帮忙运行一下pc.inp，我想确定一下是不是它的问题。。。

sobereva · 发表于 Post on 2021-2-10 04:01:13

lanl2dz赝势基组怎么能和GTH赝势搭配，且不说别的，先把这个致命问题解决了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] cp2k-5.1 单点能计算中止的问题