新版CP2K2026.2计算出cube文件，数据间无空格，导致Multiwfn处理出现闪退情况解决

JCG_666666

最近在用最新版的CP2K2026.2进行电荷密度差计算，根据社长大佬的博客过程进行处理，发现cub数据导入Multiwfn后，包括最新版，按回车都会出现闪退问题，如下图，后来发现，原始的cube数据含负号之间没有空格，最后通过代码增加了数据间的空格，再导入multiwfn，解决由空格导致不识别数据的闪退问题。

报错图片如下：



引起报错的原因如下：

数据和负号之间没有空空格导致multiwfn不识别，


通过代码增加空格可以解决cub数据导入multiwfn的闪退问题

UW_0728.

可以向上游报告：https://github.com/cp2k/cp2k/issues

sobereva

Gaussian手册里对cube格式原本的定义是数值按照E13.5来输出，我认为这样很不好，会浪费整数位，所以Multiwfn里都是按照1PE14.5E3格式输出数据，有效记录位数能多一位，而且确保会有空格、指数位始终是三位，而不会造成以自由格式读入时的兼容性问题。VESTA、ChimeraX、gview等可以载入cube文件的可视化程序也都能正常载入Multiwfn产生的cube文件。

Gaussian输出的cube文件严格满足E13.5，以如下方式输出，可见相邻数据间总是有个空格

1.   2.91780E-14  5.01961E-14  8.47212E-14 -1.40287E-13  2.27899E-13 -3.63217E-13
   
2.   5.67913E-13  8.71134E-13  1.31090E-12  1.93521E-12  2.80253E-12  3.98134E-12
   
3. ...

复制代码

CP2K当前的记录方式不能算错，也算E13.5，但强行留了三个指数位，就导致数据是负值时相邻数据之间没有空格，会造成无法被自由格式的程序载入的兼容性问题。有兴趣者可以向开发者反馈，建议像Multiwfn一样改用6(1PE14.5E3)格式来输出。

当前CP2K用户可以修改源代码里的输出cube文件的地方改一下格式来解决。