Computational evidence for a variable first shell coordination of the cadmium(II) ion in aqueous solution

In this paper, we present a state-of-the-art 100 ns molecular dynamics simulation of a cadmium(II) aqueous solution that highlights a very flexible ion first coordination shell which transits between hexa- and heptahydrated complexes. From this investigation, a dynamical picture of the water exchange process emerges that takes place through an associative mechanism for the solvent substitution reaction. Our procedure starts from the generation of an effective two-body potential from quantum mechanical ab initio calculations in which the many-body ion-water terms are accounted for by the polarizable continuum method (PCM). This approach is computationally very efficient and has allowed us to carry out extremely long molecular dynamics simulations, indispensable to reproduce the dynamic properties of the cadmium(II) aqueous solution. Quantum mechanical ab initio calculations of the hexa- and heptahydrated complexes extracted from MD configurations have revealed stable minima for both clusters with the water molecules arranged in T-h and C-2 symmetries in the hexa- and heptahydrated complexes, respectively, with a slight energetic preference for the heptahydrated one. Finally, a comparison of the calculated hexa- and heptahydrated cluster IR and Raman spectra with the experimental data in the literature, has demonstrated that the IR spectroscopy is not able to distinguish between the two species, whereas the Raman spectrum of the Cd2+-(H2O)(7) cluster provides a better agreement with the experimental data.