计算化学公社

标题: 求助：84原子内盐HOF晶胞计算静电势绘图无正常结果 [打印本页]

作者
Author: Clunos 时间: 2026-3-19 19:09
标题: 求助：84原子内盐HOF晶胞计算静电势绘图无正常结果
各位老师，新手求助下：我在做一个84原子的HOF体系晶胞计算，由内盐形式的对甲基间胺基苯磺酸构成，变胞优化后的晶胞参数如图1，晶胞如图2，绘制静电势图时遇到一些问题。我用的是CP2K做单点计算，PBE-D3(BJ)/TZV2P-MOLOPT-GTH，k点是(3,2,1)，CUTOFF=500，REL_CUTOFF=60，用Multiwfn生成了输入文件，并使计算输出两个.cube文件，运行环境Rocky Linux 10.1，CP2K 2026.1 psmp版本，虚拟机分配16核32G内存，完全使用MPI并行。画静电势时使用的是VMD 1.9.3配合sob老师的CP2K_ESP.tcl脚本，但是在isovalue=0.001时图中原子周围没有等值面，只有氢键区域有少量蓝色等值面（图2）。isovalue=0.06左右时才能看到红色以及原子周边的等值面，但图看起来也不太对劲（图3）。在VESTA中手动按CP2K培训班这部分内容设置时，图糊成了一团，并且显示只有Positive，无法选择Positive and Negative（图4）。之后也换成TZVP-GTH重算了一次单点，得到的结果没有变化。想问下我这可能是出了哪些错误，应该怎么更正？还是说这种内盐形式在输入文件处需要什么特殊处理么？
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0)

作者
Author: sobereva 时间: 2026-3-20 00:05
图1本来就是正确的，搞清楚电子密度等值面是怎么分布的。晶胞里绝大部分都是电子密度>0.001 a.u.的，自然0.001 a.u.的等值面只在远离任何原子的局部出现

作者
Author: Clunos 时间: 2026-3-20 02:27

sobereva 发表于 2026-3-20 00:05
图1本来就是正确的，搞清楚电子密度等值面是怎么分布的。晶胞里绝大部分都是电子密度>0.001 a.u.的，自然0. ...

受教了，感谢sob老师解答

欢迎光临计算化学公社 (http://bbs.keinsci.com/)