求助涉及自由基的单电子转移过程的能量计算

Orlea

各位老师，我在尝试复现文献的反应过程的数据时遇到了一些问题，反应过程如下：

这是一个单电子转移反应，用到了Marcus理论计算反应过程的能垒

根据文献中给出的公式，需要计算反应过程中重组能和吉布斯自由能的变化，就能算出反应过程的能垒。


我想复现文献中的数据，就尝试用四点法计算重组能，之前看到分子间单电子转移过程，有把A、B各取两个点计算的，也有A、B分别进行四点法计算再加和的，考虑到文献并没有交代，我就都试了一下

benzene（A-）和SO4-·(B)
对A-进行结构优化，关键词是opt freq 6-31g(d) 0 1
在A-的结构上计算A的能量E(A [A-]) 6-311++g(3df,2p) m062x scrf=(SMD, solvent=water) 1 2
对A进行结构优化，关键词是1 2
在A的结构上计算A的能量E(A [A])  1 2


对B进行结构优化，关键词是 -1 2
在B的结构上计算B-的能量E(B- [B]) -2 1
对B-进行结构优化，关键词是 -2 1
在B-的结构上计算B-的能量E(B- [B-]) -2 1
最后加和得到E(A [A-]) - E(A [A]) + E(B- [B]) - E(B- [B-])也就是λ

对于吉布斯自由能的变化，就是将Thermal correction to Gibbs Free Energy的校正量和6-311级别下的单点能相加得到自由能，然后将产物的自由能和反应物的自由能相减得到变化值
想请教一下各位老师，这种涉及自由基的过程，计算过程是不是和(A-) + B → A + (B-)这种阴离子过程一致，调整电荷数和自旋多重度就可以了呢？