sobereva 发表于 2015-7-11 22:12
我说的分裂就是指LUMO和LUMO+1的能量差。单独算一个苯的时候LUMO二重简并，能量差为0。

明白了，谢谢老师。

我说的分裂就是指LUMO和LUMO+1的能量差。单独算一个苯的时候LUMO二重简并，能量差为0。

sobereva 发表于 2015-7-11 21:04
UB3LYP对这个阴离子自由基体系恰好不收敛，UBLYP、UM062X试了一下都很容易就收敛了，ROB3LYP也能顺利收敛。 ...

感谢sobereva老师的详细讲解.

001.png和MO.png中的结构不是同一个分子的，所以会有差异。不好意思，是我没有说清楚。

我当时只改变了基组，把6-31g*替换为6-31g**发现了同样的问题，没有考虑方法的问题。在这里又跟您学到了！！

我还有一个问题，您说简并被破坏，但依然很接近，才分裂了0.035eV，您是怎么算出分裂的具体数值的呢？

谢谢老师！

UB3LYP对这个阴离子自由基体系恰好不收敛，UBLYP、UM062X试了一下都很容易就收敛了，ROB3LYP也能顺利收敛。

我在ROB3LYP下算的阴离子的SOMO，就是你001.png中的中性分子LUMO。

我不知道为什么你在MO.png把氯的朝向扭转了。

估计你知道，苯的LUMO其实是二重简并的，轨道基本特征就是你001.png的那两个图，虽然加入了基团使得简并被破坏，但依然很接近，我算了下，才分裂了0.035eV。所以此时说哪个是LUMO，其实并不很明确。加和不加SCF=QC的时候你可以仔细看看轨道能量和总体能量差异，原理上当然是能量最低的那个更合理。能量相对高的有可能因为收敛得不彻底或者收敛到了其它波函数，导致两个能量本来就很相近的轨道顺序也产生了差异。