Efficient computational methods for accurately predicting reduction potentials of organic molecules

A simple computational approach for predicting ground-state reduction potentials based upon gas phase geometry optimizations at a moderate level of density functional theory followed by single-point energy calculations at higher levels of theory in the gas phase or with polarizable continuum solvent models is described. Energies of the gas phase optimized geometries of the S-0 and one-electron-reduced D-0 states of 35 planar aromatic organic molecules spanning three distinct families of organic photooxidants are computed in the gas phase as well as well in implicit solvent with IPCM and CPCM solvent models. Correlation of the D-0 - S-0 energy difference (essentially an electron affinity) with experimental reduction potentials from the literature (in acetonitrile vs SCE) within a single family, or across families when solvent models are used, yield correlations with r(2) values in excess of 0.97 and residuals of about 100 mV or less, without resorting to computationally expensive vibrational calculations or thermodynamic cycles.