求助，ADCH和NPA电荷分析与EDA片段选取矛盾

GallonT

我在对Y3B57用EDA分的时候，在选取不同片段计算时，得到三重态Y3+(T)与三重态B57-(T)的轨道能最小，根据参考文献中的 (轨道相互作用能 Eorb 的大小是评估不同碎裂方案的适当性的有用标准。较小的 Eorb 表示所选片段与它们在最终复合物中的电子状态非常相似，因此更合适)，从而最后使用三重态Y3+(T)与三重态B57-(T)两个片段具体分析。而我在使用ADCH或者NPA电荷分析的方法对Y3B57分析电荷时，发现过渡金属Y3上都是分布负电荷，老师说我这个是矛盾的，所以我疑惑的地方是：是使用其他电荷分析来分析来符合EDA的结果，还是直接在EDA分析的时候选用例外一种轨道相互作用能不是最小的一种片段选取Y3-与B57+来分析

GallonT

还有个问题，老师，那我在用sobEDA分析的时候，对于片段的自选多重度该如何选取呢

GallonT

嗷嗷，我这里通过ADCH分析得到Y3上-0.48 ，NPA分析得到Y3上-0.33。

sobereva

“轨道相互作用能 Eorb 的大小是评估不同碎裂方案的适当性的有用标准。较小的 Eorb 表示所选片段与它们在最终复合物中的电子状态非常相似，因此更合适” 并非绝对合理。片段状态取什么来自于你的实际目的，诸如对于NaCl做分析，取中性的Na和Cl，还是取Na+和Cl-，各有各的意义，前者的Eorb可以涵盖Na->Cl的电子转移能，后者的Eorb可以体现Na+对Cl-的极化，且后者的E_ele可以更直接体现Na+和Cl-的经典静电相互作用。

另外，本来就没必然性Y3非得带正电荷。更何况两种常用的原子电荷计算方法都体现了其带负电荷。注意但凡说原子/片段电荷的时候，一定要说清楚具体数值，而不要光说正负，诸如-0.001和-0.3有天壤之别。

另外，还可以考察Y3部分的氧化态。
使用Multiwfn通过LOBA方法计算氧化态
http://sobereva.com/362（http://bbs.keinsci.com/thread-5118-1-1.html）