计算化学公社

标题: 优化晶胞报错，请老师帮忙建立输入文件 [打印本页]
作者
Author: mars936    时间: 2025-10-7 20:37
标题: 优化晶胞报错，请老师帮忙建立输入文件
我使用Multiwfn建立输入文件，使用cp2k.psmp 1840242.inp |tee 1840242.out运行程序，开始即报错，输入文件，cif文件以及错误信息见附件见，请老师查看谢谢各位老师。
作者
Author: Uus/pMeC6H4-/キ    时间: 2025-10-7 21:29
本帖最后由 Uus/pMeC6H4-/キ 于 2025-10-7 21:32 编辑

我怀疑该问题有可能和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
复制代码
作者
Author: mars936    时间: 2025-10-7 21:42
本帖最后由 mars936 于 2025-10-7 21:44 编辑
Uus/pMeC6H4-/キ 发表于 2025-10-7 21:29
我怀疑该问题有可能和MPI并行时调用ELPA库处理三斜晶胞有关，实测2025.1版的CP2K用ssmp可以跑，popt即使是- ...

感谢您的回复。我是固定晶胞仅优化H原子的位置，使用了PBEsol-D3(BJ)，因为尝试解决上述问题，本次上传的inp文件是没加D3校正的，不过都是报同样的错误。
作者
Author: mars936    时间: 2025-10-7 22:42
Uus/pMeC6H4-/キ 发表于 2025-10-7 21:29
我怀疑该问题有可能和MPI并行时调用ELPA库处理三斜晶胞有关，实测2025.1版的CP2K用ssmp可以跑，popt即使是- ...

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



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