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.

Quantum Chemistry Computation on pK_a Values of Organic Compounds in Solution

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Computing Pka Values with a Mixing Hamiltonian Quantum Mechanical/Molecular Mechanical Approach.

Calculations of Solvation Free Energy Through Energy Reweighting from Molecular Mechanics to Quantum Mechanics.

Ab Initio Calculations on Pka Values of Benzo-Quinuclidine Series in Aqueous Solvent

New Advances in Quantum Chemistry Calculation Methods of Biomacromolecular System

Quantum-chemical Predictions of Redox Potentials of Organic Anions in Dimethyl Sulfoxide and Reevaluation of Bond Dissociation Enthalpies Measured by the Electrochemical Methods.

An Approach to Computing Solvent Reorganization Energy.

Theoretical calculations of the p K a values of 1-aryl-4-propylpiperazine drugs in aqueous solution

Computation of standard equilibrium acidity of C–H acids in ionic media: shedding light on predicting changes of chemical behavior by switching solvent system from molecular to ionic

Computational Approaches to Estimating Octanol-Water Partition Coefficient and Water Solubility of Organic Compounds

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

Quantum-Chemical Predictions Of Pk(A)'S Of Thiols In Dmso

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

A New Quantum Method For Electrostatic Solvation Energy Of Protein

pKa Prediction in Non-Aqueous Solvents

DFT Study of Solvent Effects for Some Organic Molecules Using a Polarizable Continuum Model

Theoretical calculation of pKa for ammonium ions in aqueous solution

Calculations of Solvation Free Energies of Nh2-, Nh3 and Nh4+ by Using the Discrete-Continuum Model

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

Comparison of Prediction Methods for Octanol-Air Partition Coefficients of Diverse Organic Compounds.