LBO和mayer键级计算结果和键能结果趋势不一致

loovfnd

本身体系是CNO有机分子，R-离子亲核进攻A-A，A-A中间是碳碳单键。在这个计算中打算研究碳碳键的断裂过程。依据实验有两种情形，实验上并不能确定：
I：当一个的R-亲核进攻A-A时，形成RA(-)-A为一负离子，C-C键断裂，形成RA(-)和A自由基。



II：当二个的R-分别亲核进攻A-A时，形成RA(-)-AR(-)为二负离子，C-C键断裂，形成RA(-)。

分析这两种情形时，考虑了键长，键级（LBO，mayer），键能（整体优化，片段未优化）。方法B3YP/6-311G*
对比时发现：
I：键长1.63，键级（0.513，0.954），键能 51.3Kcal/mol
II：键长1.62，键级（0.559，0.970），键能 8.7Kcal/mol

分析上述结果很疑惑，
从键能上看I过程中碳碳键更稳定，但是从键级，键长上看II中碳碳键更稳定。I和II中所计算键级的差别算大吗？为什么LBO和mayer键级计算结果和键能结果趋势不一致？

还需要从什么方面考虑呢，使用那个数据更有说服力。

loovfnd

sobereva 发表于 2020-5-25 16:47
键能的定义有任意性，结果与片段状态的选取直接相关。这点在《通过柔性力常数考察键的强度》（http://sober ...

谢谢社长，重新计算一下键能，使用6-311+G**

sobereva

键能的定义有任意性，结果与片段状态的选取直接相关。这点在《通过柔性力常数考察键的强度》（http://sobereva.com/364）里专门说了。当前情况下键级更能体现当前体系中本身的键的强度，从键长上也同样反映出后者更强一些。
另外，当前计算键能的级别也太差了，有阴离子，算能量必须得带弥散函数。
另外也注意检查收敛的波函数的稳定性