求助如何理解阴离子带负电荷但所算得静电势是正值

liangsirskku

大佬们好，
小弟在算表面氯原子对钙钛矿(FAPbI3, FA=formamidium CN2H5 +)表面性质的影响时。将FAPbI3晶胞按（100）面切掉，建立slab模型。再将表面一I-离子替换为Cl-。继而限制slab深层原子的位置，优化表层结构，计算potential。
slab模型优化前后如下所示。Functional：GGA/PBE，Pseudopotential：OTFG Ultrasoft。因为简单尝试，计算质量设置为Fine。

所得FA2PbI5Cl中，Cl-平面静电势却很奇怪，如下图示，I-周围静电势是橙色负值，Pb2+周围是蓝色正值，但是Cl-周围居然也是蓝色负值。请问大佬们该如何理解这样的结果？

可以看到Cl-的形式电荷-1，实际电荷-0.54，这是比I-的实际电荷更负的（I-的结果这里没有展示）。
我也以为可能是赝势函数的影响，但是做了尝试并无差异。
我注意到偶极校正可能会影响（http://bbs.keinsci.com/forum.php ... 21610&fromuid=29697 (出处: 计算化学公社)）但是添加与否并没有改变结果。
我尝试改变Cl-为其他一价负离子X- （X=Br-, F-, SH-）结果居然类似，X离子周围电势为正值。

这里Br-和I-颜色类似，我做了标记。F-, SH-很明显能判断，无需标记。
这个正的电势的计算结果(从定性上）是否合理？如果合理，那么该如何理解？如果不合理，那么问题可能会出在哪里？或者，我对SCF Local Potential的理解是否出了偏差？

最后，恳请大佬指北。谢谢