Thermodynamic cycle for the calculation of ab initio pKa values for hydroxamic acids

A thermodynamic cycle to calculate pK_a values (Minus log of acid dissociation constants) of hydroxamic acids is presented. Hydroxamic acids exist mainly as amide isomers in the aqueous medium. The amide form of hydroxamic acids has two deprotonation sites and may yield either an N-ion or an O-ion upon deprotonation. The thermodynamic cycle proposed includes the gas-phase N–H deprotonation of the hydroxamic acid, the solvent phase transformation of the N-ion to the O-ion and the solvation of the hydroxamic acid molecule and the O-ion in water. The CBS-QB3 method was employed to obtain gas-phase free energy differences between 12 hydroxamic acids and their respective anions. The aqueous solvation Gibbs free energy changes were calculated at the HF/6-31G(d)/CPCM and HF/6-31+G(d)/CPCM levels of theory using HF/6-31+G(d)/CPCM geometries. For the proton, literature values of the gas-phase free energy of formation and the solvation free energy change were used. The free energy change for the transformation of the N-ion to O-ion in the aqueous medium was calculated by employing CBS-QB3/CPCM in the aqueous medium. For this, the hydroxamic acids were divided in two classes according to the substituent at the carbonyl carbon. A common transformation free energy difference for aliphatic substituted hydroxamic acids and a separate common transformation free energy difference for aromatic substituted hydroxamic acids were obtained. The pK_a calculation yielded a root mean square error of 0.32 pK_a units.