The dynamical processes taking place in the first coordination shells of the gadolinium (III) ion are important for improving the contrast agent efficiency in magnetic-resonance imaging. An extensive study of the gadolinium (III) ion solvated by a water cluster is reported, based on molecular dynamics simulations. The AMOEBA force field [P. Y. Ren and J. W. Ponder, J. Phys. Chem. B 107, 5933 (2003)] that includes many-body polarization effects is used to describe the interactions among water molecules, and is extended here to treat the interactions between them and the gadolinium ion. In this purpose accurate ab initio calculations have been performed on Gd3+-H2O for extracting the relevant parameters. Structural data of the first two coordination shells and some dynamical properties such as the water exchange rate between the first and second coordination shells are compared to available experimental results. We also investigate the charge transfer processes between the ion and its solvent, using a fluctuating charges model fitted to reproduce electronic structure calculations on [Gd(H2O)n]3+ complexes, with n ranging from 1 to 8. Charge transfer is seen to be significant (about one electron) and correlated with the instantaneous coordination of the ion.
