Simulations of the water exchange dynamics of lanthanide ions in 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) and water

The dynamics of ligand exchange on lanthanide ions is important for catalysis and organic reactions. Recent 17O-NMR experiments have shown that water-exchange rates of lanthanide ions in water/1-ethyl-3-methylimidazolium ethyl sulfate (water/[EMIm][EtSO4]) increase as a function of increasing charge density. The trend of water-exchange rates in this solvent is opposite to that observed in water. Since the lanthanide ions and ionic liquids investigated in that work were highly charged, an advanced polarizable potential is desirable for accurate simulations. To this end, we have developed atomic multipole optimized energetics for biomolecular applications (AMOEBA) parameters for all lanthanides and [EMIm][EtSO4], and molecular dynamics simulations with the optimized parameters have been carried out to provide possible explanations for these observed behaviors from the experiments. In water, the association of a water molecule with the first hydration shell can lead to water exchange. Smaller lanthanide ions exhibit slower water-exchange rates than larger ones because they form smaller aqua complexes, preventing the binding of incoming water molecules from the outer hydration shells. By contrast, smaller lanthanide ions undergo faster water exchange in water/[EMIm][EtSO4] because the dissociation of a water molecule is a key step for water-exchange events in this solvent. The first shell [EtSO4]− anions bind closer to the smaller lanthanide ions, resulting in more steric crowding effects on the surrounding water and facilitating the release of water molecules.