Multiwfn计算hirshfeld-I电荷迭代结果变化幅度巨大以及不明报错

Suchent

各位老师好，最近我在使用Multiwfn计算一批团簇的Hirshfeld-I电荷来预测氮气最可能在其中哪个金属上发生吸附。在计算以下的三个体系（zr4b3-,zr4b4-以及zr4c4-）时，发生了一些迭代计算错误。首先是zr4b3-,zr4b4-体系，迭代最后一步电荷变化幅度明显已经大到不太正常的地步。而zr4c4-的结果同样类似且在没有其他提示信息的情况下直接报错退出了。请问这种情况下该如何计算Hirshfeld-I电荷，或我应该换用其他原子电荷预测发生吸附的原子？（计算中的大部分rad文件是Multiwfn中example下自带的，另外补充的B+4.rad，B-3.rad，C-4.rad文件是自己在PBE0/def2tzvp级别下按照手册3.9.13计算后转换的，Multiwfn版本是今年9月3日发布的3.8版本）

sobereva

到目前为止Multiwfn 3.8版根本没发布过

Suchent

不好意思可能是我下载的时候没太注意，当时是在主页下载的Multiwfn_3.8_dev_bin_Linux_noGUI.zip的压缩包然后在服务器上安装的，这是运行开头输出的版本信息

sobereva

没有你的波函数文件和rad文件我没法判断

对于当前目的，没必要用H-I电荷，如果你想考察的是静电作用主导的物理吸附的位点，用MK或CHELPG电荷最为合适，因为它们是拟合静电势电荷，对分子外围静电势重现性最好，也因此和分子间静电相互作用关系最密切。更好的方法是按下文绘制分子表面静电势图，因为比如若最理想吸附位点是两个原子之间，以原子作为对象靠原子电荷考察就体现不了，而且表面静电势图也非常直观
使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图（含视频演示）
http://sobereva.com/443（http://bbs.keinsci.com/thread-11080-1-1.html）

你用的版本是3.8(dev)

Suchent

非常感谢老师的建议，我会尝试静电势相关电荷和表面图。计算电荷中使用的fchk文件和自己额外生成的rad文件已经在一楼上传，其余rad文件用的是Multiwfn中example文件夹下默认的