Accurate pKa Prediction in First-Row Hexaaqua Transition Metal Complexes Using the B3LYP-DBLOC Method

Abstract

Acid dissociation constants are computed with density functional theory (DFT) for a series of ten first-row octahedral hexaaqua transition metal complexes at the B3LYP/LACV3P** level of theory. These results are then scaled, primarily to correct for basis set effects (as in previous work on predicting pKa's in organic systems1-5). Finally, localized orbital corrections (LOCs), developed by fitting properties such as ionization potentials, electron affinities, and ligand removal energies in prior publications,3,4,6,7 are applied without any further parameter adjustment. The combination of a single scale factor with the DBLOC (localized orbital corrections for first row transition metals) corrections (and thus a single adjustable parameter in all) improves the mean unsigned error from 5.7 pKa units (with no parameters) to 0.9 pKa units (maximum error 2.2 pKa units), which is close to chemical accuracy for this type of system. These results provide further encouragement with regard to the ability of the B3LYP-DBLOC model to provide accurate and robust results for DFT calculations on transition metal containing species.