优化晶胞报错，请老师帮忙建立输入文件

mars936

我使用Multiwfn建立输入文件，使用cp2k.psmp 1840242.inp |tee 1840242.out运行程序，开始即报错，输入文件，cif文件以及错误信息见附件见，请老师查看谢谢各位老师。

Uus/pMeC6H4-/キ

我怀疑该问题有可能和MPI并行时调用ELPA库处理三斜晶胞有关，实测2025.1版的CP2K用ssmp可以跑，popt即使是-np 1指定单进程跑也会报同样的错，而在&GLOBAL字段添加PREFERRED_DIAG_LIBRARY SCALAPACK后popt又可以跑了。建议拿该输入文件去github的官方仓库跟开发者提issue反映。

另外，该输入文件GEO_OPT任务仅是在盒子尺寸不变的前提下对原子位置做几何优化，而如果要优化晶胞得做CELL_OPT的变胞优化任务。此外有机小分子共晶是用PBEsol合适还是PBE-D3(BJ)合适也可以考虑一下。正确使用psmp的命令应当写作

1. mpirun -np {MPI并行进程数} -x OMP_NUM_THREADS={OpenMP并行线程数} cp2k.psmp file.inp | tee file.out

复制代码

mars936

Uus/pMeC6H4-/キ 发表于 2025-10-7 21:29
我怀疑该问题有可能和MPI并行时调用ELPA库处理三斜晶胞有关，实测2025.1版的CP2K用ssmp可以跑，popt即使是- ...

感谢您的回复。我是固定晶胞仅优化H原子的位置，使用了PBEsol-D3(BJ)，因为尝试解决上述问题，本次上传的inp文件是没加D3校正的，不过都是报同样的错误。

mars936

Uus/pMeC6H4-/キ 发表于 2025-10-7 21:29
我怀疑该问题有可能和MPI并行时调用ELPA库处理三斜晶胞有关，实测2025.1版的CP2K用ssmp可以跑，popt即使是- ...

使用ssmp顺利完成了优化任务，在&GLOBAL字段添加PREFERRED_DIAG_LIBRARY SCALAPACK后popt也能顺利跑完，感谢。