Predicting Stability Constants for Uranyl Complexes Using Density Functional Theory

The ability to predict the equilibrium constants for the formation of 1:1 uranyl/ligand complexes (log K 1 values) provides the essential foundation for the rational design of ligands with enhanced uranyl affinity and selectivity. We use density functional theory (B3LYP) and the integral equation formalism polarizable continuum model (IEF-PCM) to compute aqueous stability constants for UO2 2+ complexes with 18 donor ligands. Theoretical calculations permit reasonably good estimates of relative binding strengths, while the absolute log K 1 values are significantly overestimated. Accurate predictions of the absolute log K 1 values (root-mean-square deviation from experiment