小中大【回复】
实际上你的分析本质是基于静电势来分析反应位点。然而NBO电荷(严格来说应叫做NPA电荷)对分子静电势重现得并不好,Mulliken电荷更糟(且基组依赖性还很大),建议用拟合静电势电荷或ADCH电荷来分析。更严格的方法,则是不依赖原子电荷,而是直接分析分子范德华表面上静电势的极值点位置。
而基于静电势分析反应活性本身也并不很可靠,应当再考虑使用Fukui函数、平均局部离子化势(研究它在分子表面的极值点位置)、双描述符等等指标去分析。
PS: 静电势和平均局部离子化势在分子表面的极值点位置都可以靠Multiwfn程序的主功能12实现,实例见手册4.12。Fukui函数、双描述符在Multiwfn里也都可以很容易地计算。