BOND-VALENCE METHODS FOR PKA PREDICTION: CRITICAL REANALYSIS AND A NEW APPROACH

Abstract

Bond-valence methods for the prediction of (hydr)oxide solution monomer and surface functional group acidity constants are examined in light of molecular structures calculated using ab initio methods. A new method is presented that is based on these calculated structures, and it is shown that previously published methods have neglected one or more of four essential features of a generalized model. First, if the apparent pKa values of solution monomers are to be used to predict intrinsic pKa values of surface functional groups, similar electrostatic corrections must be applied in both cases. In surface complexation models, electrostatic corrections are applied by representing a charged surface as a uniform plane of charge density, and an analogous correction can be made to solution monomers by treating them as charged spheres. Second, it must be remembered that real surfaces and real monomers are not homogeneous planes or spheres. Rather, charge density is distributed rather unevenly, and a further electrostatic correction (which is often quite large) must be made to account for the proximity of electron density to the point of proton attachment. Third, the unsaturated valence of oxygen atoms in oxyacids, hexaquo cations, and oxide surfaces is strongly correlated with acidity after electrostatic corrections are made. However, calculation of unsaturated valence for oxyacids and oxide surfaces must be based on realistic MeO bond lengths (taking into account bond relaxation), which can be obtained from ab initio structure optimizations. Finally, unsaturated valence must be divided between possible bonds (four for oxygen atoms) to reflect the fact that O-H bonds are localized to particular regions of the O atoms.