Keggin polyoxoanions in aqueous solution: Ion pairing and its effect on dynamic properties by molecular dynamics simulations

The dynamics of Keggin polyoxoanions in aqueous solution in the presence of monovalent cations is analyzed through molecular dynamics simulations. Together with structural information yielding the radial distribution functions of Li+, Na+, and K+ with three polyoxometalates (POMs) bearing 3-, 4-, and 5- charges, the diffusion coefficient of these POMs is calculated. We found that the effect of the microscopic molecular details of the solvent is a key aspect to interpreting the structural and dynamic data because a competition between electrostatic interactions between the ions and the stability of the solvation shell is established. Furthermore, we show that solvent-shared structures weakly bound to the POM anion play a crucial role in the determination of the dynamic properties of the anion. The nature of these ion pairs, structurally characterized for the first time, is consistent with experimental data available.