Theoretical descriptors response to the calculations of the relative pKa values of some boronic acids in aqueous solution: A DFT study

With the specific target of calculating the pKa values of boronic acids (RB(OH)2) in aqueous solution we inquired the solute–solvent interactions of these acids and their corresponding conjugate base. Relative pKa values were computed for each boronic acid using methylboronic acid (CH3B(OH)2) as a reference. All gas phase computations were performed at (MP2/6-311++G(d,p)//B3LYP/6-31+G(d)) level of theory. Solvation was included in the calculations using the polarized continuum model (PCM) at the HF/6-31G(d,p) level. The geometry optimization of studied structures was performed with DFT computation and the optimized structures were used to carry out Natural Bond Orbital (NBO) analysis. Natural Bond Orbital (NBO) analysis was examined as indicator for the variations observed in the calculated pKa. Natural resonance theory (NRT) was used to calculate the natural bond order and predict bond polarity. The NRT results indicate that the O–H bonds in all studied compounds have an electrovalency character. The results revealed that difference in acid strength is more relevant to the type of substituent. The electron withdrawing substituent will increase the acidic strength. However, the electron releasing substituent will have a reverse effect and decrease the acidic strength. There are good agreement between theoretical values and experimental results for this series of compounds, the average error was found to be less than 1.2 pKa unit.