求助关于弥散函数和电荷比较的问题

maxminmax

计算软件：高斯
计算目的：希望通过静电势、NPA电荷、福井函数几个参数，比较氯苯和苯酚阴离子（PhO-）中与取代基连接的碳的电正性强弱
问题：
（1）苯酚阴离子基组需要加弥散（计划只对O原子加），为了使结果具有可比性，是否氯苯分子也要加弥散？
（2）如果问题（1）中氯苯分子需要加弥散，那么可以只对氯原子加弥散吗？
（3）如果优化时不加弥散，单点时再加弥散，对上述几个参数的计算结果是否有影响？
（3）对于上述苯酚阴离子和氯苯分子来说，是否可以通过对比两者的静电势、NPA电荷、福井函数数值，说明与取代基连接的碳的电正性强弱问题？（一些资料说不能直接比较数值，只能比较变化趋势，没有想明白！！）

麻烦各位老师，抽空帮忙看一下，谢谢！！

阮文升-河内国大

（1）氯苯不需要加弥散
（2）如果加弥散，为了不加弥散的原子保持统一，应当只对氯原子加弥散。
（3）加弥散可能导致原子电荷结果更烂，因为弥散的基组可能会显著侵入周围原子。对较远的原子的原子电荷则不会有负面影响。
（4）其实我没太明白你这取代基是啥意思……

maxminmax

谢谢帮忙解答！这里取代基指的是氯和氧

sobereva

衡量原子的正电性强弱，用Multiwfn计算Hirshfeld电荷，以及Multiwfn计算局部原子表面上的静电势即可说明。没必要也不建议用NPA电荷，福井函数和这个也没直接关系。看：
Multiwfn支持的预测化学体系反应位点和反应活性方法一览
http://sobereva.com/767（http://bbs.keinsci.com/thread-58790-1-1.html）
谈谈怎么计算“原子的静电势”
http://sobereva.com/641（http://bbs.keinsci.com/thread-29998-1-1.html）
原子电荷、静电势都可以在不同体系间对比

如果你算福井函数是为了考察亲核反应性，福井函数和静电势不是很好选择，见上面的767博文

弥散函数仅对于部分原子电荷会造成不良影响，如Mulliken、SCPA、Lowdin电荷等，对于上面提到的Hirshfeld电荷、静电势分析，都没有丝毫不良影响，只不过是增加耗时而已。看下文了解不同原子电荷的特点和基组的敏感性
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！
http://sobereva.com/714（http://bbs.keinsci.com/thread-46067-1-1.html）
《原子电荷计算方法的对比》（https://doi.org/10.3866/PKU.WHXB2012281）

当前你关心的是碳，没有带显著负电荷，所以当前情况没必要给任何原子加弥散