求助，计算含AlCl3的熔盐体系，Al的原子电荷数极小

KaZikl

研究是含AlCl3的熔盐体系，温度973K，计算使用DZVP-MOLOPT-SR-GTH，PBE，NPT系综，体系总原子数302，已知的是AlCl3在973K熔盐中是以AlCl4-形式存在的，但是计算出来的原子电荷数依然很奇怪，不知道哪里出了错误，请各位老师指点，因为molden文件过大，所以文件以网盘形式给出，计算输入文件如下，不知道是不是我的设置有什么错误
https://pan.baidu.com/s/1Kz_kPpzPlp_aNIaOW-q1rg 提取码: 1234
此外学生还想问一个问题，学生那NPT跑这个体系跑了10ps，控压1bar,但是最后跑完体系约在100bar左右，盒子的边长也在20.35埃左右波动，虽然看了卢天老师那篇关于如何判断是否达到平衡的文章，但是还是不知道是否是正常现象，请各位老师解答

sobereva

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，导致信息零零碎碎，这点在置顶的新社员必读贴里明确说了。

restart文件的信息不便于直接阅读和检查设置，直接提供最初的inp文件

嫌CM5电荷太小的话就尝试其它的，下文说了
使用Multiwfn对周期性体系计算Hirshfeld(-I)、CM5和MBIS原子电荷
http://sobereva.com/712（http://bbs.keinsci.com/thread-46269-1-1.html）

KaZikl

sobereva 发表于 2025-11-9 00:15
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，导致信息零零碎碎，这 ...

谢谢卢老师解答