cp2k导出的电子密度和静电势cube计算极值点格点不匹配的问题

自留地

各位好，我利用cp2k做DFT计算导出的cube文件（周期性体系，计算需要）可以用卢老师的方法（使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图http://sobereva.com/443）在VMD中成功画出静电势图，由此可以判断cube文件并没有问题。但是，我在利用Multiwfn计算极值点时一直遇到如下问题：
Error: The grid setting of this cube file is inconsistent with that of the grid data stored in memory!

我的具体做法如下：
Multiwfn density.cub
12    //表面定量分析
5      //导入外部cube函数文件
3      //确认选择外部文件进行插值
0      //开始分析，此时会产生cube template file并提示输入静电势cube文件，我在输入esp.cub之后就出现了上述问题

我分析的原因：esp.cub和density.cub 格点是完全一致的，而template.cub是电子密度0.001 a.u.的isosurface cube，所以跟cp2k产生的不可能一致。

那么，如何产生跟template.cub一致的esp.cub呢？(实在是抱歉，cube文件过大，无法上传)

还请大家指点。拜谢

sobereva

目前主功能12不支持周期性体系，以后会支持