Guangchao Liang, Xuan Zhao, Charles Edwin Webster, Nathan J. DeYonker
Accurate computationally derived reduction potentials are important for catalyst design. In this contribution, relatively inexpensive density functional theory methods are evaluated for computing reduction potentials of a wide variety of organic, inorganic, and organometallic complexes. Astonishingly, SCRF single points on B3LYP optimized geometries with a reasonably small basis set/ECP combination works quite well--B3LYP with the BS1 [modified-LANL2DZ basis set/ECP (effective core potential) for metals, LANL2DZ(d,p) basis set/LANL2DZ ECP for heavy nonmetals (Si, P, S, Cl, and Br), and 6-31G(d') for other elements (H, C, N, O, and F)] and implicit PCM solvation models, SMD (solvation model based on density) or IEFPCM (integral equation formalism polarizable continuum model with Bondi atomic radii and α = 1.1 reaction field correction factor). The IEFPCM-Bondi-B3LYP/BS1 methodology was found to be one of the least expensive and most accurate protocols, among six different density functionals tested (BP86, PBEPBE, B3LYP, B3P86, PBE0, and M06) with thirteen different basis sets (Pople split-valence basis sets, correlation consistent basis sets, or Los Alamos National Laboratory ECP/basis sets) and four solvation models (SMD, IEFPCM, IPCM, and CPCM). The MAD (mean absolute deviation) values of SCRF-B3LYP/BS1 of 49 studied species were 0.263 V for SMD and 0.233 V for IEFPCM-Bondi; and the linear correlations had respectable R2 values (R2 = 0.94 for SMD and R2 = 0.93 for IEFPCM-Bondi). These methodologies demonstrate relatively reliable, convenient, and time-saving functional/basis set/solvation model combinations in computing the reduction potentials of transition metal complexes with moderate accuracy. © 2017 Wiley Periodicals, Inc.
As one of the most important electrochemical parameters of transition metal complexes, the reduction potential is widely used to quantitatively understand the ligand and metal effects in various chemical phenomena and catalytic reactions. However, accurate prediction of reduction potential in transition metal system encounters great difficulties and challenges. We present a computational study of a diverse set of organic, inorganic, and organometallic complexes to report on relatively reliable and inexpensive DFT functional/basis set/solvation model combinations. These combination show acceptable accuracy in computing reduction potentials of transition metal complexes.
