Quantum Chemical Benchmarking, Validation, and Prediction of Acidity Constants for Substituted Pyridinium Ions and Pyridinyl Radicals.

Sensibly modeling (photo)electrocatalytic reactions involving proton and electron transfer with computational quantum chemistry requires accurate descriptions of protonated, deprotonated, and radical species in solution. Procedures to do this are generally nontrivial, especially in cases that involve radical anions that are unstable in the gas phase. Recently, pyridinium and the corresponding reduced neutral radical have been postulated as key catalysts in the reduction of CO2 to methanol. To assess practical methodologies to describe the acid/base chemistry of these species, we employed density functional theory (DFT) in tandem with implicit solvation models to calculate acidity constants for 22 substituted pyridinium cations and their corresponding pyridinyl radicals in water solvent. We first benchmarked our calculations against experimental pyridinium deprotonation energies in both gas and aqueous phases. DFT with hybrid exchange-correlation functionals provide chemical accuracy for gas-phase data and allow absolute prediction of experimental pKas with unsigned errors under 1 pKa unit. The accuracy of this economical pKa calculation approach was further verified by benchmarking against highly accurate (but very expensive) CCSD(T)-F12 calculations. We compare the relative importance and sensitivity of these energies to selection of solvation model, solvation energy definitions, implicit solvation cavity definition, basis sets, electron densities, model geometries, and mixed implicit/explicit models. After determining the most accurate model to reproduce experimentally-known pKas from first principles, we apply the same approach to predict pKas for radical pyridinyl species that have been proposed relevant under electrochemical conditions. This work provides considerable insight into the pitfalls using continuum solvation models, particularly when used for radical species.