求助，有得电子和价态变化的反应的吉布斯自由能变数据过大

HydroKagami

如题，本人正在使用Gaussian对双核的Co配合物计算反应的吉布斯自由能变，有这么一步反应： [Co(II)Co(II)L(H2O)2]2+ +e- → [Co(II)Co(I)L]+ +2H2O
所以我将上述两个中间体以及水分子在TPSS-D3(BJ)/def2-SVP下做opt和freq计算，再在TPSS-D3(BJ)/def2-TZVP下算高精度单点，设置[Co(II)Co(II)L(H2O)2]2+的自旋多重度为3，[Co(II)Co(I)L]+的自旋多重度为2。
电子的吉布斯自由能取J.Phys.Chem.A,118,11090(2014)中的数据，即-0.001377a.u.，按照上述方程式求算ΔG，最后得出结果为-0.26305912679941a.u.，换算为kcal/mol高达-165.072232657898。（热力学校正项均为298K，1atm下的数据）
这显然不合理，请问我是在哪一步计算中出现了问题或者犯了什么错误呢？多谢大家指点。