Influence of explicit hydration waters in calculating the hydrolysis constants for geochemically relevant metals

The effect of including a second, explicit solvation shell of water molecules is examined on the DFT calculation of a selection of aquo-metal ions' pK(a) values (deprotonation constants). Our goal is to gauge the accuracy of how certain cluster approximations and implicit solvation models (PCM or COSMO) affect the results. The ions in this study include: Al(3+)((aq)); Fe(3+)((aq)); Cr(3+)((aq)); Mn(3+)((aq)); Co(3+)((aq)); Ga(3+)((aq)); Fe(2+)((aq)); Mg(2+)((aq)); and Be(2+)((aq)). Overall, we find that experimental pK(a) constants can be calculated successfully to within 1-2 pH units, provided a consistent hydration structure is maintained for products and reactants.