A theoretical study on the p K a values of selenium compounds in aqueous solution

The importance of selenium compounds in chemistry and biochemistry has become apparent. The p K a values are important for understanding the reactivity and applicability of selenium compounds. To this end, we calculated the the p K a values of 16 selenols in aqueous solution by using two kinds of DFT methods and PCM-Bondi, SMD and SMD sSAS models. It is observed that the highest precision with a lowest RMSE value of 0.38 p K units was provided by using the ωB97XD method with the SMD model. Therefore, this theoretical protocol was employed to calculate the p K a values of selenol compounds and selenoxo forms of selenocarboxylic acid compounds including R = alkyl, cycloalkyl, benzyl, allyl and aryl as well as the substituent effects. The results suggest that selenols can give higher p K a values than selenocarboxylic acids. Generally, for both substituted selenols and selenocarboxylic acids, EWGs can decrease the p K a values and EDGs can increase the p K a values, and the p K a of the EWG of –COOH is larger and the p K a values of the EDGs of –NH 2 and –N(CH 3 ) 2 are smaller. In addition, linear relationships between the p K a values with the natural charges of Se atoms, E HOMO , etc. were obtained in different types of selenols and selenocarboxylic acids.