关于自旋多重度的一点疑惑

群论是什么

各位老师好，我在利用Gaussian模拟反应计算是遇到一些疑问。计算方法和基组为：b3lyp/6-311g(d,p) em=gd3bj。过渡态命令为opt=(calcfc,ts,noeigen) freq。IRC命令为IRC(calcfc,maxpoint=20,step=10,lqa)。
反应如下图，右侧结构记为R，通过对IRC检验的两侧分别优化得到反应物和产物结构。

之前优化时R采用的默认自旋多重度为1（过渡态也取1），结果记为R1（能量-420.867011 Ha）,优化后结构如下：

我认为此处应该为3（α-β+1），结果记为R3（能量-420.840346 Ha），优化后结构如下：

但之前在学习时，接受了这样一个印象：“选择能量最低的自旋多重度为基态。”
从能量角度来说R3>R1，那这个反应中R的基态还应该是R1么？

Stardust0831

产物有两个自由基仍然可能是单重态的，可以看看《谈谈片段组合波函数与自旋极化单重态》（http://sobereva.com/82）
我的建议是对两种自旋多重度都做一下波函数稳定性测试。如果仍然是单重态能量最低，就可以说明当前产物仍然是单重态。

wzkchem5

此外产物结构不一定是双自由基，也可能是苯环上是碳正，羧基是负离子

naonao5205

在这个case里面应当是一个碳正离子和羰基负离子啊
第一个TS2，产物画的不符合基础化学知识。