Despite these confirmations, we employed the quadratic model developed by Matyushov and Voth 101,102 in order to quantify the potentially missed "non-Marcus" ∆A R→O contributions associated with differences in the solvent's thermal fluctuations upon change of the solute's oxidation state (due to e.g. the use of a non-polarizable force field for water). The calculation of this ∆A R→O term 101,102 (Equation ( ) in reference ) only requires the reorganization free energies of the reduced and oxidized states, which, as mentioned in the Theory Section, can be directly determined from the respective ∆E standard deviations.