Calculation of one-electron redox potentials revisited. Is it possible to calculate accurate potentials with density functional methods?

Density Functional calculations have been performed to calculate the one-electron oxidation potential for ferrocene and the redox couples for a series of small transition metal compounds of the first-, second-, and third-row elements. The solvation effects are incorporated via a self-consistent reaction field (SCRF), using the polarized continuum model (PCM). From our study of seven different density functionals combined with three different basis sets for ferrocene, we find that no density functional method can reproduce the redox trends from experiment when referencing our results to the experimental absolute standard hydrogen electrode (SHE) potential. In addition, including additional necessary assumptions such as solvation effects does not lead to any conclusion regarding the appropriate functional. However, we propose that if one references their transition metal compounds results to the calculated absolute half-cell potential of ferrocene, they can circumvent the additional assumptions necessary to predict a redox couple. Upon employing this method on several organometallic and inorganic complexes, we obtained very good correlation between calculated and experimental values (R(2) = 0.97), making it possible to predict trends with a high level of confidence. The hybrid functional B3LYP systematically underestimates the redox potential; however, the linear correlation between DFT and experiment is good (R(2) = 0.96) when including a baseline shift. This protocol is a powerful tool that allows theoretical chemists to predict the redox potential in solution of several transition metal complexes a priori and aids in the rational design of redox-active catalysts.