Prediction of acidity constants of some important selenium oxoacids in aqueous solution by computational techniques

The acidity constants (pK_a) of two important series of selenium oxoacids including H₂SeO_n (n = 1-5) and H₂Se₂O_n (n = 1, 3, 4, 6 and 7) were predicted by using two different computational methods for the first time. The calculations were performed using MP2 level of theory employing a very large basis set, 6-311++G(3df,3pd), and a new density functional method, M06-2X, with the same basis set, separately. A new continuum solvation model, SMD, based on the quantum mechanical charge density of a solute molecule interacting with a continuum description of the solvent, was used to account the solvent effects. The calculated values of pK_a were corrected using the different correlation equations, reported in the literature, to improve the accuracy of results. Also, the calculated results in this work were compared with the corresponding results related to sulfur oxoacids obtained in our previous work. The same trend was observed for the variation of the acidity constants of selenium oxoacids and corresponding sulfur oxoacids (with the number of oxygen atoms). Comparison of the pK_a1 of selenium and sulfur oxoacids showed that the acidity strength of sulfur is higher than selenium while for pK_a2, selenium oxoacids are stronger acids than sulfur oxoacids. The predicted values of pK_a of selenium oxoacids are important and useful in the different chemical aspects of these compounds in chemistry and biochemistry.