sobereva 发表于 2025-12-7 13:56
当前体系用Mulliken电荷说事就可以，知名度极高，对当前基组也能兼容，并且结合Multiwfn绘制的密度差图还 ...

谢谢老师！

electrochemyjy 发表于 2025-12-7 13:16
老师，我用MBIS跟Mulliken算出来的电荷转移方向结果与前面电子密度差的结果对上了，感谢老师的指点！不过 ...

当前体系用Mulliken电荷说事就可以，知名度极高，对当前基组也能兼容，并且结合Multiwfn绘制的密度差图还可以相互印证

可以一定程度说明

sobereva 发表于 2025-12-5 20:27
从密度差和电荷位移曲线来说，Pt算是稍微得了一些电子。
你这个体系，Pt和C96H24接触面积非常大，这种情况 ...

老师，我用MBIS跟Mulliken算出来的电荷转移方向结果与前面电子密度差的结果对上了，感谢老师的指点！不过我在用MBIS计算时有提示：Warning: Negative electron density has been found on some grids! In this case MBIS result will be inaccurate. Note that negative electron density at core region can be caused by PAW calculation of VASP, if it is the present case, please consider not to use PAW or use other method to calculate atomic charges。我的体系计算时用的方法是PBE/DZVP-MOLOPT-SR-GTH，请问在这种情况下是不是用Mulliken算出来的结果好点呢？后面我对这个体系进行电荷分解分析，结果显示Pt(111)向C96H24转移了2.2个电子，C96H24向Pt(111)转移了2.55个电子，总结果是Pt(111)净得到0.35个电子。这个电子转移方向与前面用MBIS、Mulliken的方法相符合。

另外，还想请教老师一个问题，我想评估C96H24在跟Pt(111)后具体是C96H24上哪些原子向Pt(111)贡献的电子更多，这个信息可以通过直接比较相同坐标状态下，C96H24在C96H24与C96H24/Pt(111)上净电荷变化来获得吗？

最后再次感谢老师的回答！

附上C96H24/Pt(111)原子电荷计算情况(将片段中原子净电荷加和起来)
Hirshfeld方法：片段Pt(111)：0.0467；片段C96H24：-0.0529
Hirshfeld-I方法：片段Pt(111)：1.472；片段C96H24：-1.478
CM5方法：片段Pt(111)：0.965；片段C96H24：-0.971
MBIS方法：片段Pt(111)：-0.985；片段C96H24：0.979
Mulliken方法：片段Pt(111)：-1.564；片段C96H24：1.564

sobereva 发表于 2025-12-5 20:27
从密度差和电荷位移曲线来说，Pt算是稍微得了一些电子。
你这个体系，Pt和C96H24接触面积非常大，这种情况 ...

谢谢老师！我这就去试下

从密度差和电荷位移曲线来说，Pt算是稍微得了一些电子。
你这个体系，Pt和C96H24接触面积非常大，这种情况下片段净电荷对原子空间划分极为敏感，Hirshfeld-I的结果未必合理，建议同时用Multiwfn算MBIS和Mulliken电荷进行对比。

-0.006是数值积分误差，直接无视