Abstract: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.

Theoretical calculation of pKa for ammonium ions in aqueous solution

Proton Solvation Free Energy in Liquid Ammonia at Room Temperature

Predicting pKa Values of Para-Substituted Aniline Radical Cations vs. Stable Anilinium Ions in Aqueous Media

Computation of Pk(A) Values of Substituted Aniline Radical Cations in Dimethylsulfoxide Solution

Accurately Predicting Protein p Ka Values Using Nonequilibrium Alchemy

Determination of Amino Acids' p K a : Importance of Cavity Scaling within Implicit Solvation Models and Choice of DFT Functionals

pKa Prediction in Non-Aqueous Solvents

Ab initio procedure for aqueous-phase pKa calculation: the acidity of nitrous acid

NMR and computational studies of ammonium ion binding to dibenzo-18-crown-6

New insights in the computational pKb determination of primary amines and anilines

X-ray Absorption Spectra for Aqueous Ammonia and Ammonium: Quantum Mechanical versus Molecular Mechanical Embedding Schemes

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model.

Equilibrium Constant Studies for Complexation between Ammonium Ions and 18-Crown-6 in Aqueous Solutions at 298.15 K

Ab Initio Calculations Of Pk(A) Values Of Transition-Metal Hydrides In Acetonitrile

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

Predicting Rate Constants for Nucleophilic Reactions of Amines with Diarylcarbenium Ions Using an ONIOM Method

Jump Rotational Mechanism of Ammonium Ion in Aqueous Solutions

Unravelling the structural and dynamical properties of concentrated aqueous ammonium nitrate solutions: MD simulation studies

A Statistical Approach To The Prediction Of Pk(A) Values In Proteins

Ab initio simulations of α- and β-ammonium carbamate (NH4·NH2CO2), and the thermal expansivity of deuterated α-ammonium carbamate from 4.2 to 180 K by neutron powder diffraction

Toward Prediction of the Chemistry in Ionic Liquids: an Accurate Computation of Absolute Pk(a) Values of Benzoic Acids and Benzenethiols.