Ab initio benchmark calculations on Ca(II) complexes and assessment of density functional theory methodologies

A set of benchmark results for the geometries, binding energies, and protonation affinities of 24 complexes of small organic ligands with Ca(II) is provided. The chosen level of theory is CCSD(T)/CBS obtained by means of a composite procedure. The performance of four density functionals, namely, PW91, PBE, B3LYP, and TPSS and several Pople-type basis sets, namely, 6-31G(d), 6-31+G(d), 6-31+G(2d,p) and 6-311+G(d) have been assessed. Additionally, the nature of the metal ligand bonding has been analyzed by means of the Symmetry Adapted Perturbation Theory (SAPT). We have found that the B3LYP hybrid functional, in conjunction with either the polarized double-ζ 6-31+G(2d,p) or the triple-ζ 6-311+G(d) basis sets, yields the closest results compared to the benchmark data. The SAPT analysis stresses the importance of induction effects in the binding of these complexes and suggests that consideration of classical electrostatic contributions alone may not be reliable enough for the prediction of relative binding energies for Ca(II) complexes.