利用RI-HFX计算Ni的素晶胞失败

guyingbing

Sob老师，各位老师：

    晚上好。我想利用杂化泛函计算含Ni氧化物的形成能。这需要得到单质镍的能量。选用杂化泛函PBE0，利用RI-HFX计算计算Ni的能量。结果在Intel 64核与Intel 96核的超算上，运行失败。具体如下，我总共测试了两个不同基组，第一个是pob-TZVP-rev2，采用全电子计算，输入文件、输出文件与报错文件分别为cp2k1.inp，cp2k1.out， cp2k1.err。第二个是采RI-HFX基础上结合ADMM技术，以DZVP-MOLOPT-SR-GTH-q18为基组进行测试，，输入文件、输出文件与报错文件分别为cp2k2.inp，cp2k2.out， cp2k2.err。两次任务都采用Ni的素晶胞，边长2.4918A，仅1个Ni原子，K点采用12，12，12。两个测试的结果非常类似都是卡在了SCF计算的第一步，随后任务就被杀掉，并在报错文件中留下响应的错误信息。以前分别是我两次测试输入文件、输出文件与报错文件。请sob老师审阅，批评，与改正。

此致
                 敬礼

sobereva

有warning，要读的波函数文件不存在，先把这个解决了
*** WARNING in qs_initial_guess.F:306 :: User requested to restart the    ***
*** wavefunction from the file named: cp2k-RESTART.kp. This file does not ***
*** exist. Please check the existence of the file or change properly the  ***
*** value of the keyword WFN_RESTART_FILE_NAME. Calculation continues     ***
*** using ATOMIC GUESS.      

如果解决了这个问题还不行，8成是内存不够。可再尝试pob-DZVP-rev2

guyingbing

sobereva 发表于 2025-6-9 23:39
有warning，要读的波函数文件不存在，先把这个解决了
*** WARNING in qs_initial_guess.F:306 :: User re ...

读取波函数之后，还是在相同的地方停止了

北大-陶豫

那你用CP2K算其它体系也有相同的问题吗？我之前也遇到过，自己装的CP2K，只要读之前算好的波函数就会在SCF第一步卡死。但是最后也没解决，自己重装不顶用，还是请工程师重装才弄好的。

guyingbing

北大-陶豫 发表于 2025-6-10 10:04
那你用CP2K算其它体系也有相同的问题吗？我之前也遇到过，自己装的CP2K，只要读之前算好的波函数就会在SCF ...

我其实也怀疑是不是超算提供的cp2k编译有问题。