Predicting the acidity constant of a goethite hydroxyl group from first principles

Accurate predictions of the acid-base behavior of hydroxyl groups at mineral surfaces are critical for understanding the trapping of toxic and radioactive ions in soil samples. In this work, we apply ab initio molecular dynamics (AIMD) simulations and potential-of-mean-force techniques to calculate the pK(a) of a doubly protonated oxygen atom bonded to a single Fe atom (FeIOH2) on the goethite (101) surface. Using formic acid as a reference system, pK(a) = 7.0 is predicted, suggesting that isolated, positively charged groups of this type are marginally stable at neutral pH. Similarities and differences between AIMD and the more empirical multi-site complexation methodology are highlighted, particularly with respect to the treatment of hydrogen bonding with water and proton sharing among surface hydroxyl groups. We also highlight the importance of an electronic structure method that can accurately predict transition metal ion properties for goethite pK(a) calculations.